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Lamps and Wavelengths of Mercury 198 
By William F. Meggers and F. Oliver Westfall 


rhe production and extraction of artificial He from ‘24Au bombarded by neutrons 
is described, and the procedure in making electrodeless lamps with He is given in detail 
hese lamps are excited by high-frequency radio waves; they emit the mercury spectrum 
entirely free from isotope shifts and hyperfine structures that characterize natural mercury 


Fabry-Perot interferometers and a prism spectrograph were employed for measuring, relative 


to Cd 6438.4696 angstroms, the wavelengths of the stronger lines in the visible spectrum of 


sole Preliminary values for '{oHg radiations in normal air are reported as 5790.6627, 


5769.5984, 5460.7532, 4358.3377, 4077.8383, and 4046.5717 angstroms, with probable errors 


not exceeding +0.0001 angstrom 


I. Introduction 


paper entitled “On the feasibility of estab- 
a light wave as the ultimate standard of 
Michelson and Morley [1]! in 1889 

bed the method later employed in measur- 
meter and reported the first interferometer 

s of wavelengths for one sodium line, one 
line, and three mercury lines. They also 
that “the brilliant green (mercury) line 
beautifully clear circles even with a differ- 
of path of half a million waves, so that in all 
itv this will be the wave to be used as the 
standard of le ngth”’ 1892 
son [2] discovered that many spectral lines 


However, in 


“the green mercury line is 
This 


tion disqualified radiations from natural 


plex and that 
the most complex yet examined”’. 


as precision standards for more than half 
vy, during which time the causes and cures 
ral line complexities were sought and 
The mass spectrometer disclosed [3] that 
mercury is a mixture of seven isotopes 
ss numbers 19€, 198, 199, 200, 201, 202, 
the interferometer resolved the compo- 
complex lines so that they could be 
tively explained [4] as resulting in part 
tope shifts and in part from hyperfine 


due to nuclear spins of ': and *; units for 
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isotopes 199 and 201, respectively. [It was obvious 
in 1931 that single sharp lines could be obtained 
from any even-mass isotope of mercury if such an 
isotope could be produced in sufficient quantity 
and purity for spectroscopic light sources. The 
latter requirements ruled out, at least temporarily, 
any physical methods of separating the isotopes of 
natural mercury. Fortunately another solution 
was found in transmutation 

Although nuclear reactions with certain elemen- 
tary particles were known to transmute some 
elements into others, such experiments yielded 
only infinitesimal amounts that could not be seen 
Thus, in 1934 Fermi, et al. 


a mixture of beryllium powder and radon as a 


or weighed. [5] used 
neutron source to transmute gold into mercury, 
but it was wholly impracticable to make tangible 
amounts of mercury this way. In 1940 Alvarez 
boldly suggested that neutrons obtained from cy- 
clotron operations might transmute sufficient gold 
into mercury to be detected spectroscopically. 
This experiment was successful and demonstrated 
that the green line of artificial mereury was devoid 
of complex structure [6, 7]. Attempts to produce 
ponderable amounts of 'SHg, required for many 
durable lamps, were frustrated until a more pow- 
erful source of neutrons became available. The 
chain-reacting uranium piles proved ideal for this 
purpose, and with the cooperation of the United 
States Atomic Energy Commission, sufficient arti- 
ficial “SHg has been made since 1945 to provide 
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several hundred satisfactory spectroscopic light 
sources Now that the isotopic curse has been 
lifted from mercury and ‘SHg has become avail- 
able in abundance, there are good reasons for 
believing “that in all probability the brilliant 


green (mercury) line will be the wave to be used 
as the ultimate standard of length’ [1, 8, 9, 10] 
It is the purpose of this paper to describe the 
preparation and operation of ‘SHg lamps and to 
report some preliminary results of measurements 


on wavelengths emitted by these lamps 


II. Production and Extraction of Artificial 


S0Lg 


Gold of high purity (999.9 fine) in the form of 
foil or chips was first placed in a fused quartz well 
attached by means of a graded seal to Pyrex tub- 
ing leading to a vacuum pump through a trap 
immersed in liquid air. The enclosed gold was 
electrically 


for about 20 hr to remove 


surrounded by an heated furnace, 
maintained at 800° C 
possible contamination by natural mercury. No 
mereury contamimation could be detected in cold 
foil purchased from the New York Assay Office 
This mercury-free gold was then submitted to 
the U.S. Atomie Energy Commission for irradia- 
tion with neutrons in a chain-reacting pile The 
production of artificial He from ot occurs 
when neutrons that attach themselves to nuclei of 
gold atoms produce a highly radioactive isotope 
half-life 2.7 days) by 8 


sly: Sy mbolically 


of Au, which decays 


emission to become stable 
i Au’+ In°—'SAu’—'SHo’+ B 


After a period for neutron bombardment and 
another for radioactive decay, the gold was re- 
turned to our laboratory, where it was treated for 
the removal of artificial “SHg in the same manner 
as formerly for extraction of natural Hg contami- 
nation. Mercury that distilled from the heated 
gold first condensed in tiny droplets on cooler 
glass walls above the furnace from whence it was 
driven to the liquid-air trap by flaming with a gas- 
air torch. To be certain that all the artificial 
mereury had been extracted, this procedure was 
continued several hours after the last visible 
traces of condensation could be detected above the 
furnace 


The yield and purity of ‘SHg are both under 
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control; the former depends on the qua v of 
gold and the time of exposure to a given ny \trop. 
flux density, the latter is inversely proportions 
to the flux density. Mass-spectrometer ayalyses 
[11] of artificial “SHg invariably show so: = 
tamination by ‘GHg, suggesting that a fraction of 
the ‘Au atoms, before decaying to 'SHe, aceres 
another neutron to become ee. whi ther 
decays to "ona 
that such "SHg contamination in impure “He 


Kessler [12] has demonstrate 


can be determined from the intensity of a forbid. 
den line. A purity of 99.6 percent, or greate; 
can be guaranteed by choosing a low heutron- 
flux density, as lower yield can be compensate 
either by more gold or by longer exposur: 

For many vears it seemed hopeless to effect 
large-scale production of pure isotopes by electro. 
magnetic separation, but this outlook was chang 
in 1945 when the colossal calutrons constructed fi 
separating uranium isotopes were applied to tl 
enrichment of stable isotopes. Keim 3] re. 
ported recently that “jHg has thus been enriche 
from 11 to 98 percent, *“3Hg from 7 to 
Whether or not 


rated natural isotopes will be produced in sufficier 


-92 percent 
these electromagnetically sepe- 
quantity and purity to compete with artificia 
sllg remains to be seen. The production an 
extraction of high-purity ‘SHg from pure ‘SAu is 
such a simple and sure procedure that it can hard 
be surpassed. 


98 


Ill. Manufacture of '35Hg Lamps 
The first SHg lamp made at this Bureau 
1946 was an electrodeless Pyrex-glass tube (5-n 
bore, 10-cm length) containing an invisible quan 
tity of ‘SHg and some argon gas; it was exeiti 
with a 100-Me oscillator constructed by Wiens {7 
and it emitted intense structureless lines of t! 
50 hr before fail 
Attempts to reyu\ 


mercury spectrum for about 


on account of clean-up [9]. 


nate this lamp by heating were frustrated by ! 


liberation of adsorbed gases whose spect 
pressed that of mercury. This experiet 
gested the use of fused quartz or Vy 
956% silica), which withstand higher t 
tures than Pyrex and are therefore mor 
out-gassed before ‘“H¢g is inserted. 

In 1947 about 60 mg of 'SHg of 99.9 
purity [11] were obtained from neutron-i" 
gold, and four experimental lamps wet! 
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with 5 mg of “SHg and 5-mm (Hg) pressure of 


pure argon in each. Two of these lamps were 
Geissler type of Vycor glass with internal elec- 
trodes of molybdenum, figure 1, A, and two were 
electrodeless type, one of Vycor and one of fused 
quartz, figure Ll. B We are indebted to N. C. 
Beese of the Research Westinghouse 
Lamp Division, Bloomfield, \. J., for making the 


Division, 


Geissler lamps and fused-quartz lamp, and for 
advice on their proper preparation 
The first experiments designed to measure wave- 


lengths emitted by these lamps soon led to the 


abandonment of the Geissler-type on account of 


lower intensity and inferior sharpness of lines 
when compared with electrodeless lamps. The 
latter are the last word in simplicity and service- 
ability; they simplify manufacture by eliminating 
sealed-in wires, graded seals, and internal elec- 
trodes, thus facilitating exhausting and outgassing. 
voltage, sparking, 
defects of 
Because these electrodeless lamps have 


In service they avoid high 


sputtering, and all other Geissler 
tubes 
been found to be highly satisfactory, their prepara- 
tion will be described in some detail 

The essentials for making ‘“{Hg lamps are, in 
addition to 'SHg, clean glass tubing, vacuum 
pumps, a bottle of pure argon gas, a low-range gas 
pressure indicator, and a gas-oxygen hand torch 
Either clear fused quartz or Vycor glass is recom- 
mended because of high melting point and trans- 
parency to ultraviolet. Unless the extreme ultra- 
violet is desired, fused quartz should be avoided 
because its high transparency to Hg 1849 A radia- 
tion permits the formation of obnoxious ozone in 
the laboratory Vycor is preferred because it 
transmits practically all Hg radiations except that 
forming ozone. The glass is first cleaned with 
hot chromic acid and rinsed with distilled water. 

To eliminate the hazard of contamination with 
natural mercury, the use of mercury pumps and 
manometers should be avoided. Evacuation was 
effected with an oil forepump and oil-condensa- 
tion pump; and for measuring the pressure of 
argon gas, either a calibrated dark space in an 
attached Geissler tube or a mechanical pressure 
gage was employed 

The procedure in making an electrodeless Hg 
lamp may be outlined with the aid of figure 2. 
Sections A and B are cleaned Vycor glass (5-mm 
bore Section A is the lamp blank with con- 


striction, a, and sealed upper end blown bulbous 
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TO PUMP 


Figure 2 Vanufacture of 
i=lamy 


so that the lamp may be viewed end-on if 

Section B is the mercury well with a sid 
provided with a V veor-Pyrex joint, bh, at 
striction,c. A broken thermometer tube n 
on the principle of the medicine droppet 
for withdrawing from the supply am 
amount of ‘SHg and depositing it as a dr 


at the bottom of section B This sectiot 
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ed supply of ‘Hg, 
nercury 


and argon flask), and section 


section B. In 


the pump 
led 
d in subsequent heating for out-gassing, 


onto making these 
il overheating of the droplet, d, and con- 
loss of ‘SOHg, is avoided by wrapping the 
f of section B with wet asbestos paper. 
vacuation section A and the uncovered 
if section Bare heated with a gas-oxygen 
ntil a nonconducting vacuum has been 

After a 
is has been admitted, the assembly AB is 
ff at 
he droplet of mercury into the bulbous end 
on A, A is sealed off at 


tion, ad, thus completing an electrodeless 


measured pressure of pure 


constriction, ¢, and inverted to roll 


Finally, section 


ontaining pure mercury and argon gas 


sence of argon gas retards the clean-up of 
increases the intensity of the mercury 
im, and facilitates operation of the lamp 


mmersed in cold, circulating water. 


Mercury Clean-Up in Electrodeless 
Lamps 


rience with these lamps since 1947, and the 
persuaded us to reduce 
content of future lamps to 3 mg or 
1946 that gradual 
limit the life of a 
but it was not known then that such clean-up 


It was recognized in 


ip of mercury would 


be retarded or corrected. Experiments 


conceived and conducted to determine the 
im mereury required for a long-lived lamp. 
lamps, each containing 1 mg of natural 
Hy 


rv and 5 mm pressure of argon, were 


red and life-tested at equal brightness on 


different high-frequency exciter units emit- 
60, and 90 Me, respectively. 
onstructed by Albert Weiss of this Bureau 
e surprised to find that, with 30 Me, 1 mg 


These units 


vas cleaned up in about 10 hr, with 60 Me 
r, and with 90 Me in 2,630 hr. This experi- 
ompted us to experiment with still higher 


370 Me 


months exhibited no trace of darkening or 


A lamp excited with for 


C1lesS 


p despite operation at a higher level of 


ess. Jacobsen and Harrison [14] have 


reported high intensity and undetectable 


p when a ‘SHg lamp was excited with 
le generated by a magnetron. 


lean-up is manifested by gradual darken- 
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ing of the glass near the terminals of the high- 
frequency exciter until black bands are formed 
the If the lamp 
is held in spring clips no darkening occurs under 


on inner wall of the lamp. 


the clips, but only outside their boundaries. Iden- 
tical effects 
Vycor, and fused silica; blackening continues until 


have been observed with Pyrex, 
all the enclosed mercury has penetrated the glass, 
whereupon the mercury spectrum vanishes and 
pure argon remains. Preliminary experiments 
with various pressures of argon indicate a marked 
dependence of Hg clean-up upon gas pressure. 
Two lamps, each containing | mg of Hg, were filled 
with 5 and 15 mm of argon, respectively, and op- 
erated at constant intensity with 30 Me until the 
Hg spectrum vanished; the first lamp glowed for 
$5 hr and the second for 365 hr. These experi- 
ments, indicating that the rate of Hg clean-up in 
ultrahigh-frequency discharges decreases with fre- 
quency and with gas pressure, are being continued ; 
The 
mechanism of Hg clean-up is imperfectly under- 


the He 


clean-up in electrodeless discharge tubes is ex- 


they will be reported in more detail later. 


stood, but following facts are fixed. 
tremely slow if the oscillatory excitation frequency 
exceeds 100 Me and if the foreign gas pressure 
Hy 
localized near the edges of the high-frequency 
black the 

This leads finally to extinction of the mer- 
light 
heat the blackened portions of the lamp in a gas- 


This 


treatment removes the blackening, recovers visible 


exceeds 5 mm In any case the clean-up is 


terminals, resulting in bands inside 
lamp. 
eury In that event it is only necessary to 
oxygen flame until they glow a dull red. 


droplets of mercury, restores the Hg spectrum, 
and rejuvenates the lamp without detriment to 
the Vycor glass. In tests of Hg clean-up the same 
lamp has been resurrected several times, thus in- 
that 
entirely eliminated, successive rejuvenations guar- 


dicating even if the clean-up cannot be 


antee practically infinite life. 
V. Wavelengths Radiated by '35Hg 


Superior spectroscopic light sources are of no 
value for length measurements until the lengths 
of several radiated light waves have been meas- 
Such measurements have been made by 
light 
so fig with the primary standard of wavelength. 
6438.4696 A. adopted [10] as the 


ured. 


comparing the waves characteristic of 


198 


prov isionally 
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Michel- 


For this purpose 


the red radiation from 


son's type of cadmium lamp 


metric value of 
we employed the method first demonstrated by 
Buisson and Fabry [15] and used for the past 40 
vears in making the most precise determinations 
of relative wavelengths. This method, as is well 
known, is based on the principle that wavelengths 
of monochromatic radiations are inversely propor- 


tional to the corresponding orders of interference 


determined from linear measurements of circular 


fringes produced by the Fabry-Perot interferom- 
eter. The only improvements upon the method 
as originally described consist of simultaneous 
illumination of the optical system by both light 
sources, and application of least squares in the 
reduction of measurements 

In our experiments the mercury lamp was 
imaged inside the cadmium lamp, and both sources 
were then focused upon the slit of a prism spectro- 
graph. Both sources and lenses were initially 
placed on the optical axis of the spectrograph by 
bisecting the beam emerging from the slit when a 
bright light was placed at the spectrograph focus. 
Then the the 


ring-projection lens were inserted in the beam 


interferometer and achromatic 
and adjusted so that the interference fringes were 
accurately centered and focused on the spectro- 
slit. both 


avoids all resulting from. slight 


graph Simultaneous observation of 


sources errors 
changes in air density or optical adjustments be- 
tween successive exposures to separate sources. 
The Michelson-type cadmium lamp was enclosed 
in a furnace maintained at 300 to 320° C., and the 
capillary was viewed end-on. The electrodeless 
mercury lamp was maintained at a temperature 
between 10° and 20° C, and was always viewed 
side-on 

The Fabry-Perot interferometer consisted of ad- 
justable etalons holding a pair of crystal quartz 
plates separated by three invar pins. The plates 
were coated with evaporated aluminum having 
reflectance of 80 percent at 6000 A. Invar separa- 
tors of 3, 5, 25, 40, 50, 67, 90, or 105 mm were em- 
ployed. In order to maintain an approximately 
constant scale of interference fringes, achromatic 
lengths were used with 


lenses of different focal 


different plate spacers; 25 em with 3 or 5 mm, 50 


em with 25, 40, or 50 mm, 75 em with 67 or 90 
mm, and 100 em with 90 or 105 mm. 
The observations reported in this paper were all 


made photographically with the aid of a Hilger 
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medium spectrograph with glass 
spectrograms were recorded on each 4 by 
plate, but no two plates were made with | 
interferometer adjustment—the orders o 
ference were changed slightly from plate 

so that different configurations of fringes 1 
Beginning with low orders, associated wit! 
separation of plates, the wavelengths pu 
for natural mercury were accurate enough 
the correct integral orders of interference 
first approximation to wavelengths chara: 
of “SHg. These 


successively refined by observing with larg 


approximate values wer 
larger orders. 

In general, for any given wavelength th 
sponding order of interference at the cent 
Fabry-Perot pattern consists of an intege! 
fraction. Because the squared values of tl 
diameters are a linear function of the ring 
ber, the fractional order is simply obtau 


prisn Fo 


dividing the square of the first-ring diameter by t! 


constant difference between squares. Fi 


more, because the entire error in the compal 


of wavelengths resides in this determination of fr: 


tional orders, it is desirable to reduce it by 


tributing it among several fringes. — Finally 


derive the most probably correct values, the ob 
tho 


vations should be reduced according to the n 


of least squares. In all our determinatio: 
‘sollg 
radiation the diameters of either five o1 
These 


with an interference comparator, designed 


wavelengths we have measured fo! 


cessive rings. measurements were 
Burns in 1913, constructed in the Bureau’s 

ment shop, and illustrated in an earlier pap: 
At first the diameters of six rings were meas 


and fractional orders were calculated by th 


cedure described and recommended by Meis 


[17]. Because this procedure involves a 
ordinate amount of arithmetic, it was aban 
ID favor of the procedure earlier deseribs d by 
and Barrell [18]. 


ally simple form when five fringes 


The latter reduces to an ¢ 
are met 
and all operations are simple mental arit! 
the findi 
final quotient that represents the fractiona 


except squaring diameters and 


sought. The squares were taken from Gau 
and the quotient from a calculating machin 

Starting with micrometered values of the 
integral orders of the 


spacers, the correct 


and yellow waves of mercury were uniqu 
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m their coincidence rates and observed 
il orders as explained elsewhere {10]. The 
air is happily heuristic for the order of 
nee because coincidences can recur only 
als of 275 waves. After the correct orders 
mercury lines had been established, the 
orresponding to other known lines, includ- 
for red radiation from cadmium, were 
ed. Finally, the wavelengths of several 
ines were determined relative to Cd 
696 A by taking this times its order and 

‘ that product by the orders of the several 

‘He lines. These are the correct relative values 
the observations have been made in stand- 


pressure 760 mm Hg and temperature 15 


provided that the optical distance between 


terferometer surfaces is the same for all 


Because, in general, neither of 


ngths. 
conditions is present, Fabry-Perot wave- 
gih determinations require slight corrections 
normal air and for dispersion of phase change 
eflection 
temperature of the air near the interferom- 
und the barometric pressure were recorded 
each spectrogram, and from these data the 
propriate air-density corrections to wavelengths 
derived directly from tables published by 
The largest corrections 


0.0014 A for 4046 A. 


tions required on account of dispersion of 


rs and Peters [20] 


type amounted to 


change were deduced from systematic 


neces obtained for wavelengths measured 
mall and with large orders of interference 

Such systematic differences observed 
»>mm and with 90 mm separations of plates 
wh in figure 3. All wavelengths (except 
\) measured with these particular plates 


& positive correction that is a function of 





5500 6000 650 
ANGSTROMS 


change on eflectance n 


phase 


fiims ¢ fed 18 d flerences n vavelenaths 


mm and S5-mm separations of plates 
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no 


interference. For 
-(.0005 A 


wavelength and of order of 
4046 A this correction 


when these interferometer plates are separated 40 


amounts to 


mm. It is obvious that these required corrections 
must be evaluated with some care if one desires 
to express wavelengths correctly within +0.0001 
A 

measurements were 


these wav elength 


When 
begun in 1947 a cadmium metal vapor lamp pur- 
chased from the Westinghouse Lamp Division 
[22] was used in lieu of the (lacking) Michelson 
lamp. The red radiation from this lamp was 
measured relative to neon standards by Meggers 
and Humphreys, who reported an increase of 
wavelength with current [8]. Such a shift in 
wavelength of the cadmium line emitted by West- 
inghouse lamps was confirmed by the present 
experiments after “GSHg wavelengths were com- 
pared with the primary standard emitted by 
Michelson-type Cd lamps, which we constructed 
and used during 1948. By measuring the wave- 
leneths of certain “He lines relative to the red radi- 
ation from both types of Cd lamps we arrived at 
the ratio of 1.000 000 23 for Cd (Westinghouse) to 
Cd (Michelson) when the former lamp is operated 
with 3 to 4 amp 

In 1949 we 


cadmium vapor lamp made by the General Elec- 


experimented with an “Osram” 
tric Company of England and were satisfied that 
the wavelength of its red radiation is practically 
the same as that from the standard Michelson 
lamp After 


ences with Michelson lamps on account of thei 


innumerable discouraging experi- 


smallness, faintness, short life, ete., and after dis- 
appointment that the wavelength of cadmium 
light from Westinghouse lamps is not constant, 
it was a pleasure to work with the improved 
cadmium lamp supplied by the General Electric 
Company of England. Unfortunately, the latter 
does not conform to construction and operation 
specifications adopted by the International Com- 
mittee on Weights and Measures [23] and by the 
International Astronomical Union [24]. 

‘sollg 


visible radiations are presented in tables 1 and 2 


Preliminary results for wavelengths of 
Table 1 contains our first results based on West- 
inghouse cadmium lamps but now corrected to the 
scale of wavelengths fixed by the Michelson type 
of cadmium lamp. The results given in table 2 
are based on Michelson cadmium lamps and are 


therefore preferable to those given in table 1 
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Only results from plate separations of 40 mm or is set by the failure of the cadmium red ra 
greater are presented; those from smaller separa- to show sharp interference fringes beyond 4 


tions are naturally less accurate. Because the waves. Even at retardations of 326,000 
error of a wavelength measurement by this method the cadmium fringes appear somewhat wir 
lies entirely in the fractional order, it would appear indistinct, as shown in figure 4, and the fra 
desirable to hold this to a minimum and increase order cannot be found with the same accu 
the integral order indefinitely. However, a limit for "soHg lines. 


TABLE 1! Prelimin , values of ‘SHe wavelengths relative to Cd 6438.4696 A in normal air 


nal angstroms with interferometer-plate separations 


Ob 


Pasir 2 lima al f "wHg wavelengths relative Preliminary measurements for the wavelengths 


vormal air of green and yellow lines of ‘SHg have been re- 
ported [25] from the National Physical Laboratory 
and from the International Bureau of Weights and 
Measures, as follows: 
NPL: 5460.7531, 5769.5985, 5790.6628 A 
IBWM: 5460.7533, 5769.5986, 5790.6630 A 
In addition to the “{QHg wavelengths presented 
in tables 1 and 2, we have measured the wave- 
length of the only argon line recorded on our inter- 
ference spectrograms; the average of 35 values 
obtained from plate separations of 40, 50, 67, 9, 
or 105 mm ts 6965.4307 A. 
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Heat of Fluorination of Cobaltous Fluoride, and the 
Heats of Reaction of Cobaltic Fluoride with Hydrogen 
and with Bis(trifluoromethyl)benzene 
(c,2,0,2’,2’,a’~Fexafluoroxylene) 


By Ralph S. Jessup, F. G. Brickwedde, and Martin T. Wechsler’ 


Calorimetric measurements of the heats of the following reactions have been made, with 


the results indicated 


IC ol c), AI1( 250°C 


L4CoF,(« CF yy (CF 
109 + 20(”) keal 


2H (g), A47(335°C 74 


i for the heats of reactions 1 and 3 are consistent, within the assigned 
with available data on the heat of formation of hvdrogen fluoride Reac 
companied by unknown side reactions, and the value reported represents the 
reactions that actually took place in the calorimeter per 14 moles of CoF, reduced 
} 


value given for the heat of this reaction, as writte! may therefore be cor 


iderably greater than that indicated by the assigned probable error 


I. Introduction The heat of the first reaction was determined at 
about 250° C and that of the second at about 335 
work deseribed in this paper was done dur- C. According to information available to us, the 
sand 1944 at the request of the Manhattan reactions proceed too slowly at lower temperatures 
ct of the U.S. Army Corps of Engineers. At to permit accurate determination of the heat 
ne the work was done the need for data of evolved, and at temperatures above 335° C, prod- 
ite accuracy Was urgent For this reason, ucts other than C,F))(CF are formed in con- 
ttle time could be devoted to development siderable amounts in the reaction of C,H,(CF 
iratus, and the purity of the materials sup- with Cok We were informed that, at 335° C, 
for the measurements was not investigated reaction 2 proceeded as indicated and was reason- 
sults obtained are, therefore, not of high ably free from side reactions 
It is believed, however, that they are The value obtained for the heat of reaction 1, 
clent interest to justify their publication namely, AJ7= —56.1 keal at 250° C. differed so 
vere requested to measure the heats of the greatly from the value, A// 159.2 keal, pre- 
ng two reactions viously reported by Jellinek and Koop? for this 
— ; ee reaction, that it was thought desirable to check 
2C oF. (ce I, (g 2CoF; (¢ 
. the value obtained in the present investigation. 
This was done by determining the heat of the 
‘LD _—- reaction, 
*s)o(g) + 14Cok o(gz) 4 
L4CoF, ( : 2CoF;(c) + Ha(g 2CoF;(c) +2HF (g 


ldre Universit { Mic nt rbor kK. Jellinek and R. Koop, Z. phys hem. A145, 30 Pat) 
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and comparing the experimental result with a 
value calculated from our value for the heat of 
reaction 1, and a value reported in the literature 
for the heat of the reaction, 


H,(g) +5 F.(g) = HF (g). (4) 
The heat of reaction 3 was measured at about 
335° C with the same apparatus used for reac- 


tions 1 and 2. 


II. Method, Apparatus, and Materials 


The method consisted in comparing the temper- 
ature changes produced in a calorimetric system 
by measured amounts of chemical reaction and of 
electric energy. 

Figure 1 is a schematic diagram of the appa- 
ratus. The calorimeter consisted of the Dewar 
flask, Q, containing liquid Aroclor No. 5460 in 
which was immersed the reaction vessel, NV. A 
heating coil of bare nichrome wire for supplying 
electric energy was wound on the insulated sup- 
ports, L, and was provided with potential leads 
for measuring the impressed voltage. By means 
of the screw propeller, A, the bath liquid was driven 
downward past the heater; through the tube, J; 
through openings in the steel support, P; upward 
between the tube, J, and the wall of the Dewar 
flask; and back into the tube, J, through holes 
located just above the propeller. 

In order to reduce heat loss from the Dewar 
flask, it was immersed in an’® unstirred liquid 
Aroclor bath in the iron cylinder, R. This outer 
bath was heated by a 660-watt contact heater, 
U’, at the bottom and by two 2,000-watt immersion 
heaters (not shown in fig. 1), which were coiled 
inside and against the iron cylinder, R. During 
a calorimetric measurement the outer bath was 
maintained manually at an approximately constant 
temperature somewhat higher than that of the 
calorimeter. Loss of heat from the top of the 
calorimeter was reduced by means of the jacket 
cover, G, of the type described by White * consist- 
ing of two sheets of copper separated by the 
asbestos board, H, and bent down at the edges so 


as to dip into the liquid in the outer bath, 2. 


In spite of the precautions taken, the heat loss 


+W. P. White, J. Am. Chem. Soc, 48, 1149 (1926); The Modern Calorim- 
eter p. 152 (Reinhold Publishing Corporation, New York, N. Y., 1928). 
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Figure l. D 


1, Steel support for stirring mo 


INCHES 
iagram of calorimeter. 


tor and shield for propeller shaft 


num resistance thermometer; C, inlet tube leading through outer 


reaction vessel in calorimeter; D, 


transite board; F, propeller shaft 
iH, transite board separating shie 


exit tube (continuation of tube, ( 
of monel metal tubing; G, copper 
lds; J, transite ring between top 


flask and copper shield (has slots to permit passage of inlet and out! 


J, copper baffle tube; K, propeller 
heater supports; N, reaction vessel 


L, heater supports; M, steel rir 
O, coiled outlet tube from react 


P, steel support for reaction vessel (permits passage of bath liq 


baffle tube, J); R, steel container for outer bath; S, galvanized iror 
the space between FR and S was filled with rock wool for insu 


support for Dewar flask; U, ring h 
fluoride powder; XY, nickel gauze 


from the calorimeter 


eater; \, transite base; W, trays 


600-mesh 


was relatively large 


is believed to have been due principally to evapor 


tion of the Aroclor. 
The outer bath was 


insulated by a 1-inch 


of rock wool around the cylindrical wall, a 


the asbestos boards, 


bottom, respectively. 


E and V, at the to; 
The jacket cover, G 
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Tray for holding cobalt fluoride powder in 


reaction vessel. 


the tube, J, were supported from the asbestos 
ward, E. 

The reaction vessel, NV, was a steel cylinder 
ontaining 29 steel trays, W’, of the form illustrated 
n figure 2. A total of about 100 grams of CoF, 
was placed in the trays, the depth of the powder 
being about 5 mm in each tray. The trays were 
arranged with their flat sides staggered, as in- 
licated in figure 1, so that the reactant gases 
flowed over each tray.‘ Above the trays was a 
i00-mesh sereen, Y, to prevent blowing of the 

fluoride powder from the reaction vessel. 
reactant gas entered the system through 
be, C, and was brought to the temperature 
outer bath before entering the calorimeter. 
actant gas was admitted to the reaction 
at the bottom, and the product gases, if 
owed from the top of the reaction vessel 
1 the coiled tube, QO, in which their tem- 
re was reduced to that of the calorimeter 


zn of the reaction vessel with trays to contain the cobalt fluoride 
i by R. D. Fowler, Johns Hopkins University. 
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before entering the tube, ), which was immersed 
in the outer bath. 

The electrical energy supplied to the calorimeter 
in the determination of its energy equivalent was 
determined from measurements of current and 
voltage with a potentiometer, in conjunction with 
a volt box and standard resistor, and measurement 
of the time of heating, using second signals from 
a Riefler clock. The temperature of the calori- 
metric liquid was measured with a 2.5-ohm 
platinum resistance thermometer, 8B, and a 
Mueller bridge. 

The fluorine, the cobaltous fluoride, and the 
C,H,(CF;). used were obtained by R. Rosen, then 
of the Kellex Corporation, from E. I. Dupont 
de Nemours & Co. No quantitative information 
concerning the purity of these materials was 
furnished, but it was stated that the purity 
was relatively high. The C,H,(CF;). was prob- 
ably a mixture of isomers. The hydrogen used 
for the measurement of reaction 3 was prepared 
electrolytically at this Bureau, and was passed 
over hot copper and through a liquid air trap 
before entering the reaction vessel. 

The drying of the CoF, powder proved to be 
difficult. Heating at 350° C in vacuum for 28 
hours did not remove all the water. For this 
investigation drying was completed by fluorinat- 
ing the CoF, at 250° C. The same filling of 
cobalt fluoride was used in all experiments. 
For use in reaction 1 it was reduced to CoF, 
by passing hydrogen through the reaction vessel 
at 335° C, and for use in reactions 2 and 3 it was 
converted to CoF; by passing fluorine through 
the reaction vessel at 250° C. 


III. Details of Procedure 


The energy equivalent of the calorimeter was 
determined before and after each determination 
of the heat of a reaction. The correction for 
thermal leakage, heat of stirring, ete. was cal- 
culated for each calorimetric experiment from 
measurements of rates of change of temperature 
of the calorimeter in “rating periods’ preceding 
and following the “reaction period.”’ 

The amount of reaction was determined from 
the gain or loss in weight of the reaction vessel. 
The hot Aroclor corroded the steel vessel, the 
rate of corrosion increasing with the length of 
time the Aroclor had been held at the elevated 
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temperatures. To reduce the error resulting from densing chambers. The C,H,(CF;). was va; 
corrosion of the reaction vessel, the following steps in a metal container immersed in an oil | 
were taken: (a) The vessel and its connecting 150° C. The tube through which the 
tubes were plated with nickel and chromium; flowed to the calorimeter was electrically 
(b) the calorimeter was filled with fresh Aroclor to prevent condensation. The rate of | 
for each measurement of the heat of a reaction; vapor into the reaction vessel during t} 


(c) the rate of corrosion of the vessel was measured, measurements of the heat of reaction 2 w; 
and a correction for the change in weight due to and 0.03 gram per minute. In order to 
corrosion Was applied to the results. The thermal complete fluorination of C,H4(CFs3)>, reac 
effect of corrosion at a constant rate would be was stopped before half of the cobalt salt 
taken into account in the correction for thermal reaction vessel had been reduced to the diff 


leakage, ete., mentioned above. The amount of 
reaction 3 was determined also from the increase IV. Results 


in weight of a tube containing Ascarite (NaOH ; 
on asbestos) in which the HF generated by the There are given in table 1 the experi 


reaction was absorbed, followed by a drying values of the heats of reactions 1, 2, and 3 rt 


agent to absorb water evaporated from the reduction of the observed value for the heat 
Ascarite. The increase in weight of the Ascarite reaction 3 from 335° to 250° C, made using dat 
tube was consistent with the corrected loss in on the heat capacities of hydrogen * and hydrog 


weight of the reaction vessel within 1 percent, fluoride,” together with the estimated values 


_ Saat > Ll « . . . far 4 
which is within the uncertainty of the measure- cal mole~' deg™' and 18 cal mole! deg rt 
ments heat capacities of CoF, and CoF,, respectively 
In the measurements of the heat of reaction 2, yielded the value A/7/(250° C) 3.6 Keal. A 
value A//(250° C 72.4 keal. was caleulat 
¥ 


combining ou: 


‘ 


the product gases were passed first through a tube 
containing Ascarite plus drying agent, and thence for the heat of reaction 3 b; 
perimental value for the heat of reaction 

AH(250° C) 56.1 keal, with the val 
AH(250° C —64.3 keal for the heat of reaction 4 
The latter value was derived from the val 


to a trap immersed in liquid air, and the increases 
in weight of the trap and of the tube containing 
Ascarite were determined. It was found that the 
sum of the losses in weight of the C, H,(CF;). con- 


tainer and the reaction vessel was equal (within mn ' ; 
PaBLe 1. Values of heat of reaction 


the accuracy of the measurements) to the gain in - 
weight of the Ascarite plus the weight of the 
fluorocarbon product condensed in the liquid air 
trap 

A small flow of helium through the reaction 
vessel was maintained before and after each 
“reaction period ~ While reactions.1 and 3 were 
proceeding, the flow of helium was stopped, the 
absolute pressure of the gases in the reaction 
vessel was held at approximately 3 atmospheres, 
and the outlet valve was opened slightly to permit 
a flow of gas from the system. This flow of gas 
was necessary during reaction 3 for the removal 
of the product gas, HF, and was found to be 
necessary during reaction 1 to prevent the im- 
purities in the fluorine from accumulating in the 


* The estimated probable errors were derived from the agreeme 


f the 


“PAC "esse 
reaction Vess¢ l. results of duplicate experiments, together with estimates o 


During reaction 2 the flow of helium through systematic errors 

° . . of reaction 2 
the reaction vessel was maintained in order to $$$ 
'F. G. Brickwedde, M. Moskow, and J. G. Aston, J. Researct 


exclude air from the system, and to carry the aaa 


products of the reaction to the absorbing and con- ¢M. Murphy and J. E. Vance, J. Chem. Phys. 7, 806 (1939 


See text for a discussion of the probable error 
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64.2 keal’ for the heat of reaction 
', together with data on the heat capaci- 
ydrogen (see footnote 5), fluorine, and 
n fluoride (see footnote 6). Our experi- 
values for the heats of reactions 1 and 3 
fore consistent with the value (see foot- 
or the heat of reaction 4 to well within the 


probable errors 


V. Discussion 
1. Reaction 1 


stated in section 1, the value A// 56.1 
btained in the present work for the heat of 
on 1 at 250° C, is in violent disagreement 
he value AH 159.2 keal at 260° C, 
d by Jellinek and Koop (see footnote 2 
tter value was derived from measurements 
equilibrium concentrations of HF and H, 
tact with CoF, and CoF, over a tempera- 
together with a 
for the heat of 


range from 470° to 595° K, 
reported in the literature 
tion of HF. In view of the consistency 
values obtained in the present work for the 
its of reactions 1 and 3 with the value, AH 
we Cc 64.3 keal, for the heat of reaction 4, 
s believed that the value reported by Jellinek 
Koop must be seriously in error. It may be 
d that the value for the heat of reaction 3, 
ited using the value reported by Jellinek 
Koop for reaction 1 and the value (see foot- 
7) for the heat of reaction 4, is of opposite 
to the value for reaction 3 found experiment- 
the present work 
AH] 52+3  keal, 


heat of reaction 1 has been reported by 


experimental value, 
et al..5 who also state that the value re- 
by Jellinek and Koop leads to values for 
fluorination of hydrocarbons with CoF; 
inconsistent as to sign with their labo- 
yperience., 
f interest to calculate the equilibrium pres- 
*, over CoF, and CoF;. An accurate cal- 
of this pressure is not possible, because 
temperature heat capacities of CoF, and 
not available for a calculation of their 
Values of Chemical Thermodynamic Properties, National 
.dards, Washington, D.C. Table 91 
W. D. Burford, III, J. M. Hamilton, Jr., R. 
iI. Litant, Ind, Eng. Chem. 38, 292 (1917 


{ Reaction 


absolute entropies. An approximate calculation, 
taking 
AF° AH 


In Pes RT RT’ 


”) 


indicates that pp, is of the order of 10~” atmos- 


phere at 250° C. In view of this extremely low 
value for the equilibrium pressure of F,, the reduc- 
tion of CoF, in reactions 2 and 3 must involve 
immediate contact of the gaseous reactant mole- 


cules with the CoF, crystals. 
2. Reaction 2 


The value given in table 1 for the heat of re- 
action 2 is subject to considerable uncertainty, 
because of the fact that the exact nature of the 
reaction between CoF, and C,H,(CF;). in the 
There Is evidence, 


calorimeter is not known. 


which indicates that, under the conditions of our 
experiments, reaction 2 is accompanied by un- 
known side reactions. Thus, Fowler et al (see 
footnote 8) report that the reaction of CoF, with 
C,H,(CF 


those in our experiments 


under optimum conditions (similar to 
gave a crude yield of 
97 percent, and a yield of pure CyFio(CFs), of 83 
percent. For the same reaction, Benner et al.,’ 
report a crude yield of 91 percent, and a recovery 
of HF of 80 percent of theory. 

In our measurements of the heat of reaction 2, 
the amount of C,H,(CF;), passed into the reaction 
vessel was found to exceed that chemically equiva- 
lent to the loss in weight of the cobalt fluoride In 
the first experiment, in which the cobalt fluoride 
was initially completed fluorinated, the excess of 
C,H,(CF;), was 10 percent, while in the second 
experiment, after partial reduction of the cobalt 
fluoride, the excess was 20 percent, although the 
rate of flow of C,H,(CF;). was much lower in the 
second experiment than in the first. 

It is believed that the value given in table 1 for 
the heat of reaction 2 represents the heat produced 
per 14 moles of CoF; reduced in the reaction, or 
reactions, which actually took place in the calori- 
metric experiments, within the limits indicated by 
the assigned probable error. This quantity would 
be useful in connection with the design of equip- 


ment for fluorinating C,H,(CF;),, and data for this 


purpose were of primary interest when the work 


R. G. Benner, A. F. Benning, F. B. Downing, C. F 
. Linch, H. M. Parmelee, and W. V. Wirth, Ind 





was undertaken. Because of the uncertainty indi- ‘during fluorination of CoF, with gaseous { 

cated above as to what reactions actually took and during fluorination of C,H,(CF;), by 

place in the calorimetric experiments, the true it over CoF;. They used a preliminary 

value for the heat of reaction 2 may differ from AH 58 keal, from the present work for | 

the value reported in table 1 by considerably more of reaction 1. Using our final value for ¢ 

than might be inferred from the assigned probable of reaction 1, their value for reaction 2 become 

error. AH 445 keal. This value is subject to the san 
Benner et al (see footnote 9) reported the value —_s uncertainty as ours with respect to the natur 

SH 460 keal for the heat of reaction 2. This the reaction to which it applies. 

value was estimated from the observed heat- 

transfer rates, in their fluorination apparatus, WaAsHINGTON, January 17, 1950. 
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Permanence of Glass Standards of Spectral 
Transmittance 


By Kasson S. Gibson and Marion A. Belknap 


Shortly after the introduction of commercially available photoelectric spectrophotom- 


eters, the National Bureau of Standards began issuing calibrated disks for checking the 


photometric scale of spectrophotometers. 


precise information was available regarding their permanence. 


the time these standards were selected no 


The present paper presents 


data showing that the glass standards kept in the files for 16 years with only occasional use 


have not changed. Data are also given s 


viasses similar to the standards when exposed under giass to south skvlight 


light) for periods up to 3 years 
I. Introduction 


One of the essential parts of a spectrophotom- 
r is the photometric equipment, and the ac- 
racy of the photometric scale is one of the major 
ms with which both the manufacturer and the 
ser of the spectrophotometer should be concerned. 
This photometric scale may be based on optical 
neiples and devices, as in all visual spectro- 
hotometers and the Hardy or General Electric 
rding spectrophotometer, or it may be based 
electrical principles and devices, as in many of the 


W hat- 


the method, the user of the spectrophotometer 


wr photoelectric spectrophotometers. 


should have some means available to insure that 
instrument is giving reliable results or to warn 
hat something is wrong. 

1 visual spectrophotometer this can be done 

ins of rapidly rotating sectored disks, whose 

es can be mechanically calibrated with high 

vy. On photoelectric spectrophotometers, 

for one reason or another, such a pro- 

is not generally feasible. An acceptable 

ve, and one that is in general much more 

ent, is the use of solutions or glass disks 

spectral transmissive properties are ac- 
known over the desired spectral range. 

st thought one might wish for a series of 

hose respective transmittances do not vary 

which will not 


ivelength, displace the 


hich do not reflect strongly, and for which 


Standards of Spectral Transmittance 


1owing the nature of the changes occurring in 


including sun- 


the transmittances can be independently deter- 
mined with high accuracy. Such filters would 
enable a test of the photometric scale to be made 
that would be independent of all other sources of 
error, and if available they would probably be 
useful. They are not available, however, and in 
fact would probably not be as useful as the series 
of moderately selective filters that are discussed 
in the present paper. This is because the non- 
selective filters would not enable the detection of 
wavelength errors no matter how great; they 
would give no information regarding the effects 
of wide slits; nor would they detect the presence 
of stray spectral energy. The selective filters, on 
the other hand, will often yield valuable informa- 
tion about these things in addition to affording 
the desired check of the photometric scale over 
most of its range. 

Shortly after the appearance of commercially 
available photoelectric spectrophotometers in 1930, 
this Bureau began the preparation and issuance of 
These 
were described at a meeting of the Optical Society 
of America in 1933 {1}. 
in other 


glass standards of spectral transmittance 


They have been noted 


Bureau publications since that date 


(2, 3, 4], and they have been listed in past and 
present Bureau test fee schedules To date about 


170 of these filters have been issued? 


Figures in brackets indicate the lite 
paper 
? Solution standards of spectral transmittancy are described 


in NBS Letter Circular L929 [3] and in NBS Circular C484 





At the time that decision was made on the types 
of filters to use, little quantitative information 
was available as to the permanence of colored 
glasses. Manganese glasses were known to be 
unstable, and the spectral transmittances of many 
glasses were known to vary in the ultraviolei 
upon exposure to strong ultraviolet radiant flux 
On the other hand it can be assumed that many 
kinds of colored glass have high stability in the 
visible spectrum, else cathedral windows would 
have paled or otherwise deteriorated to a notable 
degree through the centuries of exposure. 

Since no evidence was available that the com- 
mon types of colored glass other than manganese 
would change, the selection of the filters was made 
wholly from the standpoint of suitable spectral 


transmittances in suitable thicknesses. Of neces- 


sity it was assumed that these glasses would re- 


main constant until sufficient time hed elapsed 
that changes could be ascertained even though 
small The kinds of glass selected are described 
in the next section. For each type two glasses 
were selected for initial measurement, and these 
have been kept in the Bureau files, unused except 
for occasional measurement These remeasure- 
ments were not systematic but were sufficient to 
show that no important changes were taking place 

Naturally there were often differences in the 
third decimal between the original and the repeat 
values of transmittance, but these were usually 
slight and always within the combined experi- 
mental uncertainties. Nevertheless, considera- 
tion of these small differences made it seem desir- 
able to carry out some kind of accelerated exposure 
test in order to find out the nature of any real 
changes that might be expected This was accord- 
ingly carried out over a 3-year period 

The purpose of the present paper is to show the 
nature and extent of the changes produced under 
this accelerated exposure and to indicate the 
degree of certainty within which the Bureau 
standards (not so exposed have or have not 
changed It is believed that this information 
will be of interest to most of the users of these 


standards 
II. The Standards 


1. Description 


Four types of glass were finally chosen for these 


standards: (1) A “cobalt blue’’ glass from the 
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Corning Glass Works, designated by tl! 
merly as G55A and at present as 5551. 
been used mostly in 3-mm_ thickness 
“earbon yellow” glass from the Cornir 
Works, designated by them formerly as H 
transmission) Yellow, and at present . 
their number 3307. This has been used 
in 2-mm thickness. (3) A Jena “copper 
glass, supplied by the Fish-Schurman Corp 
and carrying the Jena number BG14. ‘I 
been used mostly in 2-mm_ thickness 
Jena “selenium orange’’ glass, supplied 
Fish-Schurman Corporation, and carryi: 
Jena number OG3. This has also be: 
mostly in 2-mm thickness 

When it is desired not to be specify 
glasses will be referred to by the color 
given above, but without quotes—namely 
blue, carbon yellow, copper green, and s 
orange. 

The usual procedure in preparing thi 
ards has been to secure a supply of each kin 
glass (all pieces from the same melt) and prey 
as many disks (29.8 mm in diameter) as poss 


from the supply at hand. From these disks 


or two would be selected and retained as Bur 
standards. The others would be kept as sta 
ards to be issued. In this paper referenc: 
be made to the “Bureau standards” or the “‘is 
standards” where it is desired to differentiat 
tween them, or simply to the “standards 
the “‘olasses”’ if such differentiation is 
portant 

The general types of transmittance cu 
these glasses are shown in figure 1. Fo 
type are shown curves for the pair of 
originally selected as Bureau standards it 
and a third curve for the Bureau standard sele« 
from the glasses now being used. In gene 
spectral transmittances of the issued sta 
are similar to but not identical with thi 
illustrated in figure 1. The broken-line 
below 390 and above 750 mu, show the s 
transmittances of these glasses in the ult: 
and near infrared, when these have been m« 

The wavelengths used in reporting va 
transmittance for these standards (sh. 
table 1) are often the wavelengths of lines 
mercury or helium spectra. These were 
originally with the intention of making the 
metric measurements at these wavelengt! 
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so 


f the four types of glass standards: 1, selenium orange 


2, copper green 3. cobalt blue 


, Car bon /é llou 


s prepared in 193 Heavy lines 


violet and near infrared spectral tr 


sources, and the Hg 404.7, 435.8, 546.1, and 
577.0 and 579.1 doublet) lines and the 


01.6 and 587.6 lines have been extensively 


such purpose. The advantage of this, 

n pointed out, is that values are obtained 
independent of any possible wavelength 
twidth errors. After extensive trial it was 
however, that the intensities at the other 
He wavelengths were so low that accuracy 
st rather than gained by their use. Al- 
ncandescent sources are now used at most 
other wavelengths in the standardization 
ments, there is still some advantage in 
¢ these wavelengths instead of going to 
went even number (for example, 471.3 
of 470) because of slightly greater certainty 
n wavelength calibrations. Furthermore, 
f the value of earlier work at these Hg 
wavelengths would have been lost by 


anew at an adjacent wavelength 
2. Temperature Effects 


probable that the spectral transmittance 
sses changes at all wavelengths if the tem- 


standards of Spectral Transmittance 


standards no 


nh 


smittance 


perature change is sufficiently large (for example, 
from — 180° C to +450° C). Illustrations of such 
changes for a considerable number of types of 
glasses have been collected and published in NBS 
Circular 484 [4]. In many instances such changes 
are important even for the usual room temperature 
variations. For use as a standard of spectral 
transmittance, accordingly, it is necessary that 
the values be obtained and certified for some se- 
lected temperature and that a formula be given 
enabling the user of the glass to compute the 


spectral transmittance of the glass at room ot 


instrument temperatures other than that specified 

Room temperatures in laboratories, if not air- 
conditioned, often vary from around 20° C on 
winter mornings to above 30° C in summet 
weather A temperature of 25° C was therefore 
selected as the most useful temperature for the 
standardization measurements, and the spectral 
transmittances of all issued standards are certified 
at this temperature 

In addition measurements were made at 40° C 
on the Bureau standards in 1933 and at 55° C 
on three of the four newly selected Bureau stand. 
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ards in 1946. These more recent measurements 
were not only made at about twice the tempera- 
ture interval used in the former measurements, 
but were also made with greater instrumental pre- 
cision. They were thus considered more reliable 
than the earlier ones and have superseded them for 
certificate use. The two sets of values were w holly 
consistent within the uncertainties of measure- 
ment. 

As already noted, at the start of the work stand- 
A suitable 


‘coefficients’’, 


ards of different thicknesses were used 


way of reporting these temperature 
d, was therefore taken simply from the relations: 


log, (i ba, +- ft 
log le ba, 
logy T's = bag 


where 7 is the transmittance at the subscript 
temperature, a is the absorbance index [5], 6 is the 
thickness in millimeters, and & is a constant 
log,o0.92, approximately) 
We define d as equal to 


logo 7’; logy Tos) 


(ly al» 


and assume that 


From this we get 


log» T's 


logio Tr, +-6d (° 


as given in the certificates sent out with the issued 
standards. These values of d are given in table | 
for the four types of standards. The assumption 
of linearity in the variation of a» with @ is justified 
by the small temperature changes involved and 
by the uncertainties of measurement of d; the 
values of d are mostly small and are generally re- 


ported to but one significant figure. 


3. Use of Standards 


In addition to the values of transmittance at 
oe" 


values of uncertainty at 


and the values of d, the Bureau reports 
the respective wave- 
lengths. These uncertainties are the “huge er- 


rors” [6], taken as 4.9 times the probable errors. 
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TABLE 1. Values of d. the te mperature coe fic en 


four types of glass used as standards of spectra 
tance. 


[For definition of d, see text 


remperature coefficient, d, of glass« 
designated as 


Carbon Copper 


vellow green 


These huge errors have been computed twice, or 
in 1933 from the values obtained on the six orig 
glasses of each type when computed to comm 
thickness, and again in 1945-46 from values o 
tained from repeated determinations on the sai 
glasses with three different types of instrume! 
[2]. 

While the number of determinations ent 
into the computations of the huge error wa 
small in either case to give a highly significa 
result, there was satisfactory agreement 
two sets of data, and there has been no e 
that the values reported in the respective c 
tificates do not cover the over-all uncertainty) 
measurement. These reported values of 
tainty are usually close to those given in 
2 to 4. 

The method of obtaining the values of 
mittance for the issued standards has varied 
what over the years, changing with instr 
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Beckman photoelectric 
and the 
boundaries designated as 2U represent the 


spectrophotometers and values adopted in 1933 from measurements made on the Kénig-Martens 
were obtained by use of stray energy filters in the The 
It is concluded that the standard has not char | 


attached 
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FIGURE 


values of transmittance obtained in 1949 on the Beckman photoelectric (large circles 


slues adopted in 1933 from me 


wdopted 
ted values. It 


and personnel. At present the values for the 
issued standards are measured on the Beckman 
spectrophotometer (2, 7] by direct comparison 
with the carefully calibrated Bureau standards of 
the same type. The uncertainty of this compar- 
ison is very small, and the uncertainties as reported 
in the certificates are always large enough to cover 
this final step. In addition, an absolute deter- 
mination of transmittance is usually made on each 
issued standard on the Kdénig-Martens visual 
spectrophotometer [8] at a few wavelengths, this 
serving continually also as an indirect check on 
the stability of the Bureau standards. 

The intent in the use of the issued standards 
by those receiving them should be not to correct 
the photometric scale of the spectrophotometer 
but to check it. If the instrument gives values 
on the standard that are in agreement with those 
reported, within the combined uncertainties, the 
instrument presumably is functioning correctly 
and can be expected to give reliable results. If 
the values obtained are outside of the combined 
uncertainty, the user should attempt to eliminate 
the cause of the error. The issued standards are 


hardly adequate to use as a means of correcting 
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values having undergone slight revision in 1940 The boundaries 


carbon yellow glass standard kept in the Pureau files with only occasional use 


and the Kk 


surements made on the Kénig-Martens and the Gibson spectroy 


designated as 2[ 


s concluded that the standard has not changed 


the photometric scale, a possible exception to t! 
occurring for low values of transmittance. | 
example, the error found by use of the standa 
might be caused by stray radiant energy, w 
slits, or wavelength errors. If important, « 
should attempt to eliminate the cause, not atten 
to correct the photometric scale, which may 
wholly reliable 

Regarding selection of these standards by thos 
who wish to check the photometric scale of t! 
spectrophotometers, it is generally found that | 
are sufficient, (1) the carbon yellow and (2 
the cobalt blue or the copper green, depending 
the desires of the user. The cobalt blue is 
for checking the photometric scale near th 
T=0.01 or —logyy7=2.00. The selenium « 
has not proved useful for this purpose 
purchase is not recommended. Only a fe 


been issued. 


III. Permanence of Unexposed Samples 


The pairs of glasses originally sele 
retention as Bureau standards have bee! 
mostly in the files, with only occasional | 
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excessive exposure to radiant flux. In 
se glasses had been carefully measured 
,onig-Martens visual and Gibson photo- 
9] spectrophotometers at the various 
wavelengths from 390 to 750 my, and 
f spectral transmittance adopted at that 
\fter the intervening 16 years, one of each 
airs of glasses was remeasured in 1949, 
comparison of these results with those 
in 1933 is given in this section. 
glasses selected for this comparison are 
ited as Corning G55A a" (cobalt blue), 
r HT Yellow a’ 


copper green 


carbon yellow), and Jena 
The spectral transmit- 
of these glasses are shown in figure 1 (the 
of the two light-weight curves in each case). 
ited above, the Jena OG3 glasses have 
to have little use as spectrophotometric 
uds, and no investigation of their perma- 
has been made The most certain test of 
nence would have been to remeasure these 
on the same instruments and under the 
nstrumental conditions as were used in the 


ginal measurements. The older methods are 
however, particularly the visual, and in 
hopes of expediting the work the three sam- 
es were measured on the more rapid Beckman 
ctrophotometer. These measurements were on 
absolute basis, not relative to any other stand- 
and were made with the samples at 25° C as 
They did, in fact, 


w that the samples had in all probability not 


the previous Measurements. 


ged, but it proved necessary also to make 
additional measurements on the Kdénig- 
us for greater certainty, particularly at the 
ransmittances. 
results could be illustrated in various ways, 
most suitable appears to be that used in 
2, 3, and 4. In these figures there are 
the differences between the values of spec- 
nsmittance measured on the Beckman or 
Martens spectrophotometers in 1949, and 
ted in 1933 from the measurements made 
Kénig-Martens and Gibson spectropho- 
These differences are designated as 
and are plotted as circles. In the same 
here are shown the over-all uncertainties, 
at are usually reported on the issued 
s and that would apply closely to the 
values. If the 1949 values fall within 
certainties there would be no certainty of 
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change in the transmittance having occurred. 
Furthermore, at least an equal over-all uncertainty 
should be assumed for the 1949 values, so that 
deviations up to twice the plotted uncertainties 
would not give certain indications of change. 

As may be seen, the values of 7y— 73; in most 
cases fall within the over-all uncertainties plotted. 
In those instances where the values fall outside, 
the evidence, such as it is, is nullified by the other 
measurements shown. 

With special reference to figure 4, Corning HT 
Yellow a’, 


vised slightly in 1940 at the longer wavelengths 


the values adopted in 1933 were re- 


on the basis of some measurements made at that 
time, and in figure 4 the values of 74,— 7% are 
also shown. It is not believed that the sample 
changed between 1933 and 1940, but it was felt 
in 1940 that the revised values were slightly more 
accurate than the 1933 values. Note, however, 
that the values of 7,.— 733 exceed the uncertainties 
only at 750 mu. 

It seems fair to conclude that the spectral trans- 
mittances of these standards have not changed 


over this 16-year period. 
IV. Accelerated Exposure Tests 


It was desired that the accelerated exposures 
should be made under conditions not too widely 
different from ordinary indoor conditions. The 
samples were accordingly placed in a small box 
with a plate glass cover, which was kept outside 
of a south floor. The 
samples were thus exposed to daylight, including 
sunlight, from the south half of the sky, in all 
However, the plate glass and 
The plate 
glass further served to keep the solar ultraviolet 
radiant flux of the shortest wavelengths from the 


window on the fourth 


kinds of weather. 
box kept the samples dry at all times. 


glasses, since conditions would otherwise not have 
been comparable with those indoors. Small holes 
in the box permitted circulation of air in and out 
of the box and helped to prevent condensation of 
moisture on the plate glass cover. 

These exposure tests were made with carbon 
vellow, copper green, and cobalt’ blue glasses 


similar to the respective standards. No corre- 


sponding tests have been made of the selenium 


orange glass. Measurements on the three glasses 
were made prior to the start of the exposures and 


further measurements have been made after the 
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ncluding sunlight 


ttance before exposure 
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500 


© 


~f) 


700 


WAVELENGTH, my 


ange of 


$ pec 
skylight 


pet iods indi 
C 


“ated All measurements were made 
o* 


at £0 


The results are illustrated in figures 5, 6, and 7 


Differing from figures 2 


to 4, the ordinates of 
are plotted as ratios rather than dif- 
These ratios are the values of transmit- 
measured after the 


these figures 
ferences. 
tance 


indicated exposures 
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ttance of a cobalt blue glass 


incl 


to the 
siding 8 enlight 


similiar standards 


(pon exposure under 


divided by the respective values of transn 
before the exposures began. Since all o 
after exposur 
made on the Beckman spectrophotomet: 


measurements, and 


initial 
tive to a closely similar reference glass, t! 
cision of measurement was higher than obt 


in figures 2 to 4, in which the measureme! 
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600 


WAVELENGTH, 


hon yellou 


skylight ne 


und the trat 


on an absolute basis, no reference standard 
nvolved 

all in 
the 
and in fact the consistency of the data in 
that 


percent 


changes indicated 
100+ 0.01 


s probable that 


es 5 to 7 outside value are 


cases indicates the changes are real 
though less than 1 
the 
from figures 5 to 7 


that 


wavelength 


to be 
is that, in view of the 
390 


aps most important conclusion 
and 
the 


one 


changes occurred between 


the range for which 


standards are customarily certified, 
the to show 
In this 


to 7 may 


not issued standards 


expect 


with ordinary laboratory use 
the data of figures 2 to 4 and 5 
to be consistent, and thus the assumption 
or perhaps we should say the result hoped 
933 has been corroborated. 

to 7 be 


v details relating to figures 5 may 


the ultraviolet between 334.2 and 400 mu 
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glass 


luding 8 


the shows a large 


solarization effect 


copper green glass, figure 5, 
of the usual type, the carbon 
a small one, and the 
all It is possible, 


of course, that the cobalt blue glass might solarize 


vellow glass, figure 7, cobalt 


blue glass, figure 6, none at 
at shorter wavelengths. Those interested in this 
effect should read the papers by Coblentz and 
Stair, who made an extensive study of the ultra- 
violet solarization of window glasses and othe 
materials [10, 11]. 

2. The changes in the cobalt-blue glass, figure 6, 
are small throughout, with extremes barely ex- 
ceeding 1 percent. While some of the changes 
shown are probably real, further work would be 
to indicate whether the rather erratic 
behavior the 


The changes, if real, are certainly not important 


necessary 


shown by curves are significant 


for this purpose, and no further measurements are 
contemplated. 
3. The changes in the copper green glass shown 


by the curves of figure 5 are undoubtedly real, but 


again are of little importance (above 390 mu) in 


view of the extent of the exposure necessary to 
bring them about 

+. The major changes shown by the carbon yel- 
low glass of figure 7 are also undoubtedly real 
Again these would not be judged of much im- 
in this connection for the fact 


portance except 
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WAVELENGTH, Mu 


kiGuRE 8 é tion of pure oO ange from 


I 
shorter 


ind 53 days 


1 definite change 


that all of the change occurred during the first 
summer of exposure, there being in fact a slight 
recovery shown by the later measurements. 
Another sample was accordingly selected and 
measurements made at the intervals shown in 
These 


change between 400 and 500 my starts immediately 


figure 8 results show that the large 


upon exposure and is half completed within 8 


days efter start of the exposure. One would 


therefore recommend that this glass not be left 
around where sunlight or other strong radiant flux 


might irradiate it. At that, however, only at 


$35.8 my do the changes shown in figures 7 and 8 


greatly exceed the over-all uncertainties shown 


in figure 4.' 


‘ This is shown by three exampk The transmittance of this glass at 
435.5 = 0.040 1.07 X0.040=0.0428. 0.0428—0.040=0.0028, which is four times 
the tolerance shown in figure 4 at 435.8 mg 2) At 420 mse T=0.085 and the 
resulting change in 7'=0.0025 compared with an uncertainty of 0.002 in figure 
‘ }) At 301.6 mp T=0.2s and the resulting change in T'=0.004, compared 


with an uncertainty of 0.008 in figure 4. 
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€T pos res 


nall to | 


Acknowledgment is made to various indiy 
Most 
original measurements were made by Gerald 
Walker Haupt and Mabel E. Brown, as indicat 


Recordings of the sp 


who have cooperated in this work 


by reference [1] 
transmittances of the standards were ma 
Rece 


brations of new standards of the same typ: 


various times by Harry J. Keegan. 


glass were made by Margaret M. Balcor 
Lois Peterson Gloster, as indicated in re 
{2}. Miss Peterson also made new and 


precise determinations of the change of tra 


tance with temperature than had been p 
with the previous instruments, and ma 
initial measurements on the exposure tests 
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§ 600-Ohm Multiple-Wire Delta Antenna for Ionosphere 
Studies 


By H. N. Cones, H. V. Cottony, and J. M. Watts 


This paper describes the design and the performance of a multiple-wire delta antenna 
developed for use with high output impedance, vertical-incidence ionosphere sounding equip- 
ment. Graphs are presented showing the terminal impedance of this antenna over the sep- 
arating frequency range (1 to 25 megacycles), using various lengths of open-wire transmission 
line The results of pattern measurements using model techniques are given, showing the 
expected radiation characteristics of the full-scale antenna. A practical test of the com- 
parative merits of the antenna is described in which ionosphere records obtained by the 
use of this antenna are compared with those obtained by the use of a larger antenna developed 


previously 


I. Introduction ing requirements above, the physical dimensions 


J, ; of the antenna should be moderate, at least in the 

luring the past 20 years the technique of 
vertical dimension. 

eal-incidence ionosphere soundings by the : 
’ Numerous experimenters have worked on the 

ol pulse transmitters scanning over a wide . . . 
problem of the design of antennas suitable for 

ieney range has been adopted internationally ' 
' ionospheric recorders. Systems of several an- 
recording the virtual height of ionosphere 
; ane © tennas, each operating over a fairly narrow fre- 
rs and their degree of ionization. These are A 
a quency range and switched automatically by the 
ised in predicting the conditions for com- “tity ' it 
' equipment, have been tried. Some systems ob- 

ation in the high-frequency band. One , ee 
ey tained satisfactory operation over a fairly wide 

latest developments in the equipment em- 
, we frequency band by employing low-(Q cage an- 

iv this technique is exemplified by the model “- ae 
Aare ine age tennas. Even antennas that changed their di- 
automatic ionosphere recorder. This in- 
: mensions in the course of operation are known 

ent was developed at the Central Radio ; os 
: to have been tried. The most generally satisfac- 

agation Laboratory of this Bureau. It ; ‘ 
Cdl nine torv results have been obtained by the use of non- 

through the frequency range of 1 to 25 ; ; 


y - resonant antennas, such as rhombic and vee an- 
an interval of time as short as 7's secs. ps ’ 
. tennas. These are particularly suitable for use 
performance requirements imposed by the , a : 
, “hi with the model C-2 ionosphere recorder because 

ent on the antenna system, including the “ght ; : 
; of their inherently high input impedances. Ex- 

ited transmission lines, which are to hold ; 
ta an perimental work has been carried out at this 
the entire frequency range of 1 to 25 Me, , : 
is Bureau looking toward the development of an 
follows: (1) The impedance of the antenna eee 
, effective antenna system for ionosphere sounding. 
must be relatively uniform and, for pur- ; hi ; 
. Some oO e results 0 Ss work ec e r bee 

of maximum power transfer, as nearly me of the results of thi rk have already been 
“ call 8 a —— : i. aint 

0 the output impedance of the transmitter published.’ Here, the results of the work per 


order of 1,000 ohms) as possible; (2) the formed since the date of that paper are described. 
n must be substantially in the vertical 

(3) the antenna must be as efficient H. N. Cones, Impedance characteristics of some experimental broad 
° o,° band antennas for vertical-incidence ionosphere sounding, J. Research N BS 
or as possible. In addition to the operat- 43, 71 (1949) RP 2006 


n; 


’ 
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II. Instrumentation 


As mentioned previously, the important aspects 
of the antenna to be considered are its impedance, 
its radiation pattern, and its radiation efficiency. 
The last factor may be defined as the ratio of the 
power radiated in space by the antenna to the 
This 


non- 


power delivered to it by the transmitter. 


factor is of particular significance in a 


resonant, terminated antenna where actually 
more power is likely to be dissipated in the ter- 
minating resistor than radiated in space. 

For measurements of antenna impedance, a 
balanced recording impedance meter developed at 
this Bureau was used. A description of this in- 
strument was presented at the Winter 1946-47 
IRE Convention in New York City.’ 


is the block diagram of the instrument. 


Figure | 
Briefly, 
it is a beat frequency generator including a 57-Me 
fixed-frequency oscillator beating with a variable 
frequency oscillator covering the band from 57 to 
82 Me. The difference frequency of 0 to 25 Me 
is passed through a wide-band amplifier designed 
to have a constant current output. The output 
voltage is, therefore, directly proportional to the 
absolute magnitude of impedance connected across 
the output terminals. A voltage stabilizing circuit 
is incorporated for the purpose of maintaining 
reasonably constant output current. This voltage 
is rectified, amplified, and applied to a recording 
milliammeter. The recording milliammeter and 
the variable frequency oscillator are mechanically 
coupled, so that the impedance is plotted as a 
function of frequency. No direct indication of the 
phase angle of the impedance has thus far been 
provided. The instrument is calibrated by mak- 
ing several frequency sweeps using fixed nonin- 
ductive resistors across the measuring terminals. 

For measurement of the radiation pattern of the 
antenna, a model antenna technique has been 
employed. The principles of this technique con- 
sist briefly of constructing a model antenna having 
the configuration of the full-sized antenna under 
study, but with the linear dimensions reduced by 
some scaling factor n. The radiation pattern of 
this antenna is then measured in a conventional 
manner at a frequency n times the operating 
frequency of the full-sized antenna. It thus be- 


comes possible to bring the target transmitter 
?H. V. Cottony, A method of rapid continuous measurement of antenna 
mpedance over a wide fre paper presented before IRE Con 


vention (New York, N. ¥ 


juency range 
March 6, 1947 
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much closer physically to the antenna un 
and still maintain a spacing of severa! 
lengths at the scaled-up frequency. The 
sional analysis of this technique shows 
addition to scaling the linear dimensions «: 

a factor n, and the frequency up by n, the : 
tivity should be increased by the same fa 
This scaling factor should be applied to all e! 

of the radiating system. The radiating systen 
not merely the antenna itself, but also the er 
and, insofar as the conductivity is concerned 
scaling factor should, properly, be applied to | 
the material of the antenna (copper) and the so 
It is patently not possible to select a mater 
having conductivity exceeding that of copper 
any scaling factor n. However, at least in 
present case where copper losses may reasona! 
be expected to be small compared with th 
ated power and the terminal resistor losses 
significant error should be introduced = b 
failure to scale the conductivity. 

In the case of the ground, the conductiy 
should be sealed if strictly accurate results ar 
be obtained. However, no completely Satish 
tory method for scaling the ground constants 
available. Furthermore, since the ground co 
stants change from one location to another, a 
in general, will also change to some extent i 
given location with time, no attempt was mad 
in course of these tests to seale the ground « 
stants. Instead, the ground was covered by 


wire netting that effectively provided pert 
reflector conditions at all frequencies. This 
course, is an approximation of the actual operat 
conditions. The principles of model techniq 
are now well known and are described in the t 
nical literature.* * 

The actual model range employed in thes 
consists of a level gravelled surface, covered 


-in. wire mesh. In practice the model o 
antenna was used as the receiver, while a targ 
transmitter, powered by self-contained batte: 
was swung about the antenna ia an are appro) 
mately 13 ft in radius. To simulate the frequ 
cies of 9 to 25 Me, a sealing factor of n equal 
30 was employed, the target transmitter operat 


between 270 and 750 Me. 


To simulate operat 


on investigation of model te« 


Ernest A. Jones, Final report 


determination of the characteristics of low-frequency antennas 
University Research Foundation (Columbus 10, Ohio 

‘George Sinclair, Theory of models of electromagnetic 
IRE 36, No. 11, 1344 to 1370 (Nov. 1948 
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5 and 9 Me, a sealing factor of 85 was 


2 illustrates the terminology employed 
erence to the radiation patterns presented 
paper. It is convenient, for reference pur- 
o locate the antenna under test in the 
ane. That plane and the Y-Z plane at 
ngles to it are the principal planes of the 
a, and all radiation patterns were taken in 
the other of the two planes, and are so 
fied. In making the patterns of the antenna 

\-Z plane (path Y-Z7-X’ 


ation of the target transmitter was always 


the plane of 


dent with the -Z plane, and its equatorial 
tic plane always passed through the Y-Y’ 
For patterns in the )-Z plane, (path Y-Z- 
he .Y-’ axis was always in the plane of 
jarization, whereas the equatorial magnetic 
ine was coincident with the Y-Z plane 
No direct measurements of radiation efficiency 
imade. However, a measure of over-all per- 
nance of the antenna (including its radiation 
ency) was obtained in a series of tests, in 
ch the antenna under discussion was compared 
th an earlier type of antenna and with a num- 
ol simple doublets In these tests an iono- 
recorder was alternately switched from the 
-wire delta to a double-W antenna (see 
wtnoete | In each test a number of records 
re obtained and the results compared with each 
her. This method does not permit measure- 
ent of small differences but is a useful check of 


e effectiveness of the antennas. 


III. Design of the Antenna 


Nonresonant antennas, which are sometimes 


terminated or traveling-wave 
three types: (1) 
Of these. the 


we antenna is suitable primarily for low- 


ed to as 
are, in general, of 


) 


(2) vee, and (3) rhombic. 


radiation. Of the vee and the rhombic 
as, the vee antenna can be shown to be 
iniformly effective as a radiator over the 
ney range than the rhombic antenna. 
er, the conventional-type vee antenna with 
it terminals at the apex and the terminat- 
stors at the open end of the vee presents 
mechanical difficulties when designed to 
in the vertical direction. It requires two 


tht poles for its installation, and additional 


wires must be furnished for termination to ground. 


For these reasons rhombic antennas have fre- 
quently been used for ionosphere sounding. 

A variation of a vee antennna commonly known 
as a delta was devised several years ago by J. W. 
Cox of the Baddow Research Laboratories of the 
British Marconi Co. 


that the termination is made at the vertex of the 


Its principal features are 


vee, and the extremities of the vee are brought 
together by two horizontal wires extending a few 
feet above ground. The radiation due to these 
horizontal wires is minimized by their proximity 
to the ground, and the radiation pattern should be 
approximately that of a vee antenna 

Practical considerations require a delta antenna 
to be designed for a much lower height than one 
having theoretically optimum size and proportions 
To be very efficient a nonresonant antenna should 
be several wavelengths long. Since the equip- 
ment operates at frequencies down to 1 Me, and 
since the desired direction of radiation is straight 
up, the antenna must necessarily be short (mea- 
sured in the desired direction of propagation) as 
compared to the wavelength except at the upper 
This height 
was limited by the lengths of the readily available 


end of the operating frequency band. 
telephone poles to 70 ft. In addition, it was highly 
desirable that the antenna be capable of being 
connected to the equipment by transmission lines 
that might differ in length for each ionosphere 
station according to local conditions. In order 
that the transmission line need not accentuate the 
ratio of maximum to minimum impedances of the 
antenna proper, it is necessary that the geometric 
mean of these impedances be very close to the 
characteristic impedance of the line. The geo- 
metric mean of the impedance of a simple single- 
wire delta antenna is of the order of 800 ohms 
An open-wire transmission line having this value 
of surge impedance has awkwardly large wire 
spacing. A 600-ohm transmission line is, on the 
other hand, a very common one for which standard 
commercial spacers are available. Furthermore, 
600 ohms is a suitable load impedance for the type 
C-—2 ionosphere recorder. It was, therefore, de- 
cided to reduce the mean impedance of the antenna 
to 600 ohms. Such reduction of impedance is 
readily accomplished without affecting either the 
radiation pattern or the mechanical construction 
The use of 


multiple-wire construction offers the additional 


by the use of multiple-wire elements. 
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advantage of being able to reduce the impedance 
variations over the frequency range by varying 
the spacing of the wires that compose the individ- 
ual elements 

The height of the apex having been set by 
practical considerations at 70 ft, it is possible to 
compute the base width of the antenna so that it is 
wide enough to give useful radiated power at 1 
Me and yet not so wide as to have the first null or 
minimum within the operating frequency range of 
the equipment. Guided by such computations, 
a number of scaled-down antennas having various 
base widths were constructed and tested on the 
model range. The final design selected was one 
with the largest base width, which had no nulls 
or minima in the vertical direction at any point in 
the frequency range. The width so selected was 
130 ft measured corner-to-corner. 

Following the determination of the major dimen- 
sions of the antenna, the task remaining consisted 
of adjusting the minor configuration of the antenna 
construction, so that the mean terminal impedance 
of the antenna would be as nearly 600 ohms as 
possible and that the impedance of the antenna 
would be as uniform over the whole frequency 
The approach to this task was 

By using multiple-wire con- 


range as possible. 
largely empirical. 
struction, the impedance was lowered; increasing 
the spacing between the wires provided an addi- 
tional control of impedance. Further uniformity 
of impedance was obtained by varying the spacing 
between the wires along the length of the antenna. 
The principle of reducing the reflection along the 
electric path by varying the spacing between 
electrically paralleled conductors is well known, 
and the technique of its application is strictly 
empirical. 

The end result of the experimental work is a 
600-ohm multiple-wire delta antenna. Figure 3 
shows the principal dimensions of this antenna. 
For comparison purposes a ‘“‘double-W”’ antenna 
(see footnote 1), developed previously, is illus- 


trated in figure 4 


IV. Performance of the Antenna System 


The performance of the antenna system was 
evaluated by measuring its terminal impedance 
with three different lengths of 600-ohm_ trans- 
mission line; by measuring the radiation patterns 


of a model antenna constructed to scale; and by 
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using the antenna in conjunction with a model C 
automatic ionosphere recorder for making jono. 
sphere records. 

Figure 5 presents the terminal impedance of thy 
multiple-wire delta antenna system with thre 
lengths of 600-ohm transmission line, and, fo 
comparison purposes, the terminal impedance o 
the double-W antenna including the exponentiy 
line transformer. As can be observed, the meg 
impedance of the antenna for all lengths of trans. Hyode! 
mission lines is approximately 600 ohms, and ¢| Up 
variations of impedance are between approximately 
400 and 1,000 ohms. The performance of thy 
doube-W antenna from the standpoint of termina 
impedance is substantially the same; the mea 
impedance is also very nearly 600 ohms, and ¢| 
impedance variations are between approximately 
400 and 900 ohms. 
variations in the case of the double-ll 


The greater number of imped- 
ance 
antenna is caused by the longer transmission-lin 
system of that antenna. The exponential lin 
transformer used in the latter is approximatel) 
288 ft in length. 

Figure 6, a and b, illustrates the expect 
radiation patterns of the multiple-wire dels: 
antenna between the frequencies of 5 and 25 Mi 
as obtained on the model antenna range. In th 
frequency range from 5 to 9 Me a sealing facto 
of 85 was used, and in the frequency range fro! 
9 to 25 Mea scaling factor of 30 was used \t 
Me the results obtained by the two sealing factors 
can be checked against each other. It should } 
noted that the patterns shown are normaliz 
i. e., the gain of the pattern plotter was adjust 
for each frequency to give a full-scale deflectio 
at the point of maximum gain 

Examination of the radiation patterns shows 
that no significant side lobes are developed in tl 
patterns for frequencies below 11 Me. Some 
ing is observed at the frequencies near 14 M 
The maximum of the major side lobe is at 54 
the horizontal. It is believed that this lol» 
caused by radiation from the horrzontal wir 


On the basis of the model measurements, serious 


lobe splitting would not be expected below 2 \ 


Above this frequency the patterns should rapid) 
deteriorate However, even at 25 Me the: 
he expected to be considerable useful radiat 
the vertical direction 

No patterns have been taken at frequ 
corresponding to 12 and 13 Me because no 
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rmins 
mea 


had ti 


mately 


TN pe 
ible -} 
yn-lit 


ul lin 


ter Was immediately available to operate 

ponding scaled frequencies. The patterns 

at frequencies corresponding to 9 Me 

using sealing factor of 85 and 270 Me for 

(0) agree reasonably well. It should be 

hat this is the poimt of least satisfactory 

ration for the model range since at a frequency 

Me the target transmitter is only three and 

wavelengths away from the antenna 
inder test 


pon completion of the scale model tests, a 


| 
ut il 


test of antenna performance was at- 


wmpted by comparing actual ionospheric records 


ile using the various antenna systems. The 
nparison Was principally between the multiple- 
delta and the double-W antennas, both used 
transmission; initially horizontal doublets of 
ous lengths were used for reception. A check 
s also made of the relative merit of the delta 
fenna as compared with a doublet when used 
eption and of the effect of the length of the 
nsmission line on the effectiveness of a multiple- 
delta antenna 
Unfortunately, the reflecting characteristics of 
onosphere may vary in a period of seconds, so 
measurements made by comparing received 
sc amplitudes by switching from one antenna 
nother are subject to error. However, it was 
sidered that a series of as many as ten con- 
records using one antenna could be com- 
d with two other series using another anten- 
me recorded immediately before and the 
immediately afterward. Both receiving 
nas and transmitting antennas were com- 
separately by switching only one of the 
nas at a time. Figures 7 through 11 show 
of the ionosphere records made for these 
risons All of the records for each com- 
were made with the same receiver gain 
inless otherwise noted 
re 7 shows the multiple-w ire delta compared 
double-W’ antenna when used for trans- 
The records were made in the early 
hours so that observations on low fre- 
would be more likely A long receiving 
was also used to favor the low frequencies 
i appearance of these records it is con- 
hat the multiple-wire delta is substantially 
performance to the double-W antenna for 
ssion between 1.5 and about 3.5 Me 


8 compares the same antennas as the 


na for Ionosphere Studies 


previous figure, except that it is intended to com- 
pare the double-W and multiple-wire delta an- 
tennas at higher frequencies. This comparison 
was made with a shorter antenna for receiving 
and during a daytime period when critical frequen- 
cies were high. The conclusion is that these two 
antennas are also practically equivalent at frequen- 
cies between 3.0 and 12.0 Me. 

Figure 9 is the same as figure 8, except that a 
long receiving antenna was used. Both sets of 
records exhibit weaknesses that are primarily a 
characteristic of the receiving antenna. No con- 
clusions were drawn from these records, but further 
tests illustrated by figure 10 were suggested by 
them. 

Figure 10 shows the difference between the 
single-wire delta and a 500-ft doublet for receiv- 
ing. The receiver gain setting was the same for 
each run. It may be noted that much greater 
pickup of both ionospheric signals and noise is 
obtained with the delta. 

Figure 11 is a comparison between the effect of 
using long and short transmission lines to connect 
the transmitter to the multiple-wire delta. For 
receiving, a single-wire delta was erected in a 
plane perpendicular to that of the transmitting 
antenna. The records made with a 200-ft line 
and with a 20-ft line are substantially alike, and 
the slight differences observable are no greater 
than the variations from record to record caused 


by changes in the ionosphere 


V. Conclusions 


On the basis of the preceding discussion the 
following conclusions were reached: 

1. The impedance of the multiple-wire delta 
antenna is reasonably flat in its characteristics 
over the operating frequency range of 1 to 25 Me 
the ratio of maximum to minimum impedance 
being of the order of 2.5. This compares satis- 
factorily with the double-W antenna 

2. The mean impedance of the multiple-wire 
delta antenna is sufficiently close to 600 ohms to 
permit satisfactory operation with a 600-ohm 
transmission line of any moderate length. This 
is also confirmed by the service tests using the 
ionosphere recorde! 

3. The model antenna study of the multiple- 
wire delta indicates that the radiation pattern to 


be expected from this design is such that the 
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direction of maximum intensity of radiation will 
be at the desired zenith angle over all of the 
operating frequency range below 23 Mc. How- 
ever, even at 25 Me. a substantial proportion of 
power is radiated in the vertical direction. 

4. The ionosphere records obtained by the use of 
the automatic recorder indicate that, within the 
limits of the accuracy of this method of compari- 
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Figure 0 Daytime comparison of multiple-wire delta with double-W antenna, each used for transmission 


4 S00-ft horizontal doublet, 25 {t above ground was used as a receiving antenna 
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Mass Spectra of Diborane-d, and Ethane-d, 
By Vernon H. Dibeler, Fred L. Mohler, and Laura Williamson 


The complete mass spectra of B,D, and C,D, and the partial spectrum of B,D,;H have 


been obtained for ionizing voltages of 50 volts and 70 volts 


An analytical procedure is 


given for deriving the contribution of all ions containing a hydrogen atom from the poly- 


isotopic spectrum of a mixture of B,D, and B,D;H. 


A comparison of spectra of B,D. and 


BH, and of C,D, and C,H, indicates that substitution of deuterium for hydrogen has 


only a slight effect on the dissociation probabilities in diborane and ethane 


is the occurrence of a metastable transition 


B.D,*—B.D 2D 


An exception 


in B,D, that does not occur in BoH, viz 


Simple statistical considerations, based upon the assumption of equal probabilities for 


the removal of a deuterium or a protium, give a dissociation pattern that differs by more than 


experimental error from that observed 
I. Introduction 


{ number of theories have been proposed [1]! 
the structure of electron-deficient molecules, 
these, the 
vdregen bridge model and “protonated double 


wticularivy the hydroborons. Of 


md” model have received support from recent 
asurements of the infrared absorption spec- 
im of diborane [2] and the mass spectra of 
borane [3] and pentaborane |4]. In the case of 
wrane, a possible structural correspondence 
th ethane is ruled out and the presence of an 
thylenic linkage is indicated. It is also apparent 
im the mass spectra of diborane and pentaborane 
hat not all hydrogens are bound in the molecule 
equivalent bonds. 
In the present investigation, we have obtained 
mass spectrum of diborane-d, for comparison 
th normal diborane to obtain information on the 
stability of the molecule ion as well as the effect 
iterium substitution on dissociation processes 
orane., 
s convenient at this time also to report the 
spectrum of ethane-d,. This spectrum is 
ired with normal ethane to demonstrate 
fect of complete deuterium substitution in 
Similar data for 


arbons. methane and 


in brackets indie rature references at the end of this 


Mass Spectra of Diborane-d, and Ethane-d, 


Turkevich> 


Friedman, Solomon, and Wrightson [5]. 


methane-d, have been published by 


II. Experimental Details 


Mass spectra were obtained with a type 21-102 
Consolidated mass spectrometer. The tempera- 
ture of the ionization chamber was controlled at 
245° C 
50 and 70 v. 

The B,D, was synthesized according to Finholt, 


lonizing electron energies were nominally 


Bond, and Schlesinger [6]. Approximately 0.3 g 
of 98-percent lithium aluminum deuteride ? [7] 
were dissolved in 25 ml of anhydrous dibutyl 
ether contained in a flask provided with a reflux 
condenser and connections to a vacuum system. 
Solution was accomplished after several hours of 
stirring under an atmosphere of nitrogen at a 
temperature just below the boiling point of the 
ether. The solution was frozen with liquid nitro- 
gen, the flask evacuated, and approximately 1.5 g 
of boron trichloride distilled onto the frozen solu- 
tion. The flask and contents were allowed to 
warm slowly to room temperature, stirring when 
possible. The reaction was moderate, and after 
completion, the volatile material was drawn off 


The deuteride, deuterium gas, and deuterium oxide used in the prepara 
tions of BeD>« and CyD¢ were obtained through the Isotopes Branch of the 


USAE( 
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by a Toepler pump through a trap at 20° C. 
A mass spectrometric analysis showed no hydro- 
borons heavier than diborane, although approxi- 
mately 7 percent of B,D,H was present. The 
BCI, originally obtained from Kahlbaum in 1906 
was quite yellow and probably contained some 
HCl. This may for the 
expected amount of protium in the product. The 


account more than 


rate of decomposition of B,D, to give deuterium 
gas and heavier hydroborons is apparently slow 
even at room temperature, as a sample showed no 
appreciable increase in noncondensibles or heavier 
hydrobarons over a period of several weeks. Be- 
fore admitting a sample to the mass spectrometer, 
however, the precaution was always taken to pump 
off possible noncondensibles at liquid nitrogen 
temperature 

The ethane-d, was kindly 
Wall of this Bureau. It was made by the cata- 
Mass spectrometri: 


supplied by L. A 


lyzed reduction of acet vlene-d 


analysis indicated impurities of approximately 
| percent of deuteromethane, 2 percent of deutero- 


propane, and | percent of deuterobutane, in addi- 


tion to approximately 1 percent of C,D,H. The 


isotopic purity of nearly 99 percent is exceptionally 
good for a compound containing six deuterium 


atoms 


Ill. Results 


The computation of the monoisotopic spectrum 
of B D Is complicated by the presence of B.D H 
The method used, however, is similar to that 


already described for B,H, [3]. By trial approx. 
mation, the abundance ratio, B’®/B"=0.250. wa 
found to give the best fit. This is in agreeme 
with the normal dihora 
(0.251) and (0.250). Usin 


ratio 0.250, one derives the abundance ratio 


ratio obtained for 


pentaborane 


B’B"/Bi=0.500 and Bi°/B'=0.0625 
With these values it is possible to calculate ¢} 
spectrum of ions containing only deuterium ai 


As re lat IV 


few molecules contain protium atoms, the firs 


of ions containing | protium atom. 


spectrum is to a first approximation the B,] 
spectrum, and the second spectrum is the part 
spectrum of B,D;H accounting for all ions +) 
contain protium. There is no basis for dist 
cuishing between ions from B.D, and from B,D. 
with protium removed 

Table 1 shows the B, portion of the B.D 
B.D,H mixture spectra at 70 v ionizing voltag 
and summarizes the calculations of the comp 
monoisotopic spectrum of B,D, and the par 
spectrum of B.D.H Column 1 lists the n 
numbers of the recorded jons; column 2 gives 
mixture spectrum with the ion intensities rela! 
to the maximum peak at m/e=32 taken as 
The monoisotopic spectra are separated by st: 
ing at the heavy mass end and applying 
abundance values calculated above to determ 
the contribution of each B, ion to each va 


‘ 


mie The peak at m/e=34 comes entirely 


B.DoH 
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ion; therefore, the B'’B"D,* ion 
0.17 <0.50=0.08 and the B,’D,* ion 

0.17 <0.0625=0.01. The difference 

n the calculated and observed values at 
‘3 is in the region of experimental uncer- 
since this peak is very small. The B,'D,;* 
100—0.01=99.99. The principal contribu- 
from B,D, but with a small contribution 
B.D,H. The B®B'"D 31) is 
0.50=50.00, and the B2’D;* ion (m/e=30) 
0.0625 measure 
B:'D,H the 
50.00=6.80. From this value we compute 
ontributions of the B’°B"D,H* and B2°D,H 
The B.D, 
the contributions of B{’D 


the 


ion (Mie 


6.25. We now have a 


ion, which is difference 


s as before ion is now obtained 


subtracting and 


B DH from 
,.40 $6.50 


30 peak; I e 56.15 


The B°’B"D, and BD, 


re computed and the cvcle continued until 


obtained 
The 


ndance of the 21 and 20 peaks is calculated to 


contributions to all peaks are 


are two checks on the calculations 
0.82 and 0.10, respectively, in exact agreement 
h the observed abundance 

the 
group for lonizing 
the 
m/e=17 is assumed 
Therefore, the 
250=0.22. By 
0.07 
After 


check 


‘able 2 is a summary of monoisotopic 
tra calculations of the B 
ron energies of 70° \ Starting with 
st mass in that group, 
the B'D 
DD ion (m/e 16) 1s 

ition, the B'D.H 

he B°®D.H* ion (m/e=15) is 


mtributions are 


entirely ion 
O.S4 0 
ion is 0.29 2% 
0.02 
obtained, the single 


Vonoisotopic dissociation patterns of ByDs, BoH 
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point is m/e=10. The observed discrepancy is 
considerably larger than the experimental uncer- 
tainty, requiring a B"/B" ratio of approximately 
0.26 to obtain a zero residual. <A similar effect 
was noted and discussed in a previous paper [4]. 
The B,D, and B,D,;H dissociation patterns are 
obtained by the contributions of all 
ions with the same number of deuterium atoms 
The most 
B,.D,;H 
Other values of ion 
The 
complete spectrum of B,D, and the partial spec- 
trum of B,D,H for 50- and 70-v electrons are given 
The spectrum of B,H, is included for 
arbi- 


summing 


regardless of boron isotopic species. 
abundant B.D, (B.D,;*) 
(B.D,H*) are taken as 100 
abundance are relative to the maximum peak. 


ion of and of 


ip table 3 


comparison Sensitivities (ion current in 
trary units per unit of sample pressure) are given 


for the maximum peaks of B,D, and B,H,. 


caicula- 


y of the 


B.D 


monoitsoto pu 


and BoD.H for 


TABLI 2 Summa spectrum 


tions for the B ‘O-v electrons 


group of} 


B2DsH * 


DH =BD.H 





Omitted from the B,D, spectrum in table 2 is 
a small peak at m/e=12.5 amounting to 0.36 
percent of the maximum at 50 v and 0.40 percent 
at 70 v. This peak is due to the double charged 
ion, B°’B'D,**. The two isotope peaks expected 
coincide with large peaks at m/e=12 and 13. A 
similar doubly charged peak appears in B,H, at 
m/e=11.5, amounting to 0.39 percent of the 
maximum at 50 v and 0.43 percent at 70 v. This 
doubly charged ion of mass 23 was ascribed [3] to 
either BYH** or B’B"H,**. It is very likely 
B’B"H,** corresponding to the B’B"D,** peak 
in the B,D, spectrum. 

Also omitted from the B,D, spectrum are the 
four small peaks listed in columns 2 and 3 of 
table 4. These are the result of dissociation of 
ions after they have traversed the electric field [8]. 
The peaks at 21.6 and 22.6 are isotopic as indi- 
cated by their relative abundance and are due to 
the transitions B’B"D,*—B"B"D,*+2D and 
BUDS —-BYD,+2D. The peaks at m/e=23.6 
and 24.6 are also isotopic and are a result of the 
transitions B“’B"'D,* -B"B"D,* +2D and B,D; 

-BYD,*+2D. The transition BYD;*—-BYD;* 
2D may also contribute to m/e=22.6, but the 
contribution is very small and below the thresh- 
old of observation. The two metastable transi- 
tions observed in diborane are also given in table 4 
These are also isotopic and result from the transi- 
tions B“’B"H,* >B"B"H,*+2H and BYH,*—B, 
H,*+2H. The apparent mass of each metastable 
ion is calculated from the usual relation m, 
(m,)?/m,, where m, and m, are respectively the ion 
masses before and after metastable dissociation. 
The calculated and observed values agree within 


0.1 mass unit. 


TABLE 4 Percentage abundance of metastable ions relative 


to the maximum peak in Bo Dg and BH, 


BoD 


Table 5 gives the dissociation patterns of ethane 
and ethane-d, at 50 and 70 v. The sensitivities of 


the maximum peak in each compound are also 
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given. The peaks at the odd mass numbers 


the C,D, spectrum have been corrected for exrby 


isotope and represent contributions from an esi. 
mated 1.5 percent of C,D;H in the sample. | 
two peaks indicated by the footnote very prob. 


ably include contributions due to the doubl 
C.D,** 
doubly charged ions appear in C,H, at m/e=13; 
and 14.5. The relative abundance of these jops 
is given in the C,H, spectrum. Three metastab| 


charged ions 


and C,D;**. Analogoys 


ions not included in table 5 are observed in t| 


C,D, spectrum 


One 


appearing near mass 24; 


amounting to 0.19 percent of the maximum peal 


at 50 v and 0.18 percent at 70 v is a result of | 
C,D,.*+2D. Another, nea 
m/e=26.5, of 0.19 percent at 50 v and 0.17 perce 
from the transition C.D 


transition C,I 


)," 


at 70 v, comes 


C,D,*4-2D. The third, near m/e=28.4, of 0.08 
percent at 50 v and 0.09 percent at 70 v, comes 
from the transition C.D,*-C,.D,*7+2D. Ana 


ogous transitions 


m/e=24.1, 25.1, 


are 


and 26.1. They are 0.10, 0 


observed for ethane 


and 0.30 percent, respectively, of the maximw 


peak at 50 v 
at 70 v. 


and 


0.12, 0.17, and 0.30. perce: 


TaBLe 5. Monoisotopic dissociation patterns of ©,D 


CH, for 50 and 70 v electrons 


CDs 


® Includes contributior 


f 


rom the 


sensi 


tivity 
doubly charged ions C;D 
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IV. Discussion 


idenced by table 3, the mass spectra of 
e and diborane-d, are quite similar. The 
; in B,D, is considered a good measure of 
}.D,* ion as an impurity on this peak is 
vy in diborane. The lower value of the 
abundance of the B.D," ion (0.18) com- 
to the B,H,* ion (0.34) is perhaps not 
icant. The latter value suggests that a 


of nitrogen impurity may contribute to the 


§ peak of normal diborane. Although the 


ile ion is very unstable, it is evidently not 


mpletely absent from the mass spectrum as 


first supposed [9]. The sensitivities of the maxi- 


im 


reg 


peaks of both compounds agree within 3 
nt, which is the estimated uncertainty, 


monstrating the nearly equal probability of 


ssociating a deuterium from B,D, and a protium 
m BH, In the B, group the relative prob- 


yal 


| 


ies of dissociating successive deuterium 


toms from B,D, and successive protium atoms 


from B,Hg are similar but always somewhat less 


tl 
p 


ires 


1 deuterated compound. In the B, group 
robabilities are roughly equal. This com- 
reasonably well with the similarities pre- 


ously reported [5] and those observed for ethane 


| 
(i ¢ 


\etastable transitions with the loss of mass 2 


thane-d, discussed below r 
) 


by far the most frequent type found in hydro- 


ns [10]. It is not surprising, therefore, to 


d this transition readily occurring in deutero- 


al 


iormal diborane. An interesting difference 


etween B.D, and B.H,, however, is the occur- 


’ 


yf 


of a metastable dissociation of the B,D,;* 
diborane-dy, whereas the analogous transi- 
s absent in normal diborane. As observa- 


of a metastable ion indicates a delay of 
ghly 107° second in the dissociation process, 


evident that substitution of deuterium for 


im in diborane has materially changed the 


life of the transition involved. The other 


etastable dissociations are similar and have 


ul 


\ 


equal probability in diborane and 
ine-d,. 
formula B'B"D,** is unambiguously as- 
to the ion appearing at m/e=12.5 in the 
ine-d, spectrum. The similarity of relative 
lance of this peak and the doubly charged 
pearing in diborane at m/e=11.5 suggests an 
interpretation of the latter that is not 


Spectra of Diborane-d, and Ethane-d , 


possible from the B,H, spectrum alone. From 
the foregoing, it is probable but not certain that 
the doubly charged peak in B,H, has the formula 
B’B"H,** analogous to the doubly charged 
ion in B.D,. 

The part of the B,D;H spectrum that contains 


protium atoms has been included in table 3. 


The method of computation and the small mag- 
nitude of the residuals from which the pattern is 
derived give a relatively large uncertainty for this 
pattern. As expected, the pattern is not like the 
patterns of either B,D, or B,H,. One would 
predict that the abundance of the ions B,D,;* plus 
B,D,H* from B,D,;H should be compared with the 
abundance of the B.D,;* ion from B.D,; the abun- 
dance of B,D,* plus B.D,;H* should be compared 
with the abundance of B.D,* from B.D,, and so 
forth. Ifthe probability of removing a dueterium 
or protium atom is merely a statistical matter, 
then the abundance ratio B.D,H*/B.D;* from 
B.D;H would be 5to1. From table 1 the B.D,H 
peak is 6.8 percent of the B,D;* peak, and on the 
above basis, B.D;H contributes 1.36 percent to the 
B.D,;* peak. As the B,.D;H* ion is immeasurably 
small, there is no check on this. The probability 
of removing a deuterium rather than a protium 
from B,D,H* is 4 to 1, and it might be expected 
that the pattern coefficient of B,D,;H* from 
B,D,H* is 4/5 of the pattern coefficient of B,D, 
from B,D;*. The value thus calculated (37) is 
not in agreement with the observed value of 44 
The probability of removing two deuteriums 
rather than a deuterium and a protium from 
B,.D,H* is 3/4 4/5=3/5, and one would pr dict a 
pattern coefficient of 14 compared with 18 ob- 
served. Similarly, for B,DH* one would predict 
0.4 66.9=27 instead of 29, and for BDH*, 8.8 
instead of 10. For B,H* one would expect 0.9 
instead of zero, and for BH*, 2.4 as observed. 
Analysis of the mass spectra of simple hydro- 
carbons containing one deuterium atom [5, 11] 
indicates that removal of the deuterium and 
protium atoms is not equally probable but in a 
ratio of about 0.5 to 1 for these molecules. It is 
not surprising, therefore, that purely statistical 
considerations fail to account for the B,D,H 
spectrum. The case of diborane, however, is 
more complicated as, according to recent theory 
[1], the six hydrogen atoms are not equivalent, 
but two of the atoms are bound loosely to an 
ethylene-type structure. 
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Insofar as the statistical method roughly ac- 
counts for the spectrum, the following conclusions 
can be drawn. The amount of B,.D;H in the 
mixture is about 8.5 percent. The contribution 
of this to most peaks in the B,D, pattern amounts 
to a few tenths of a percent of the B,D;* peak 
(The B,D, pattern is with B.D,* 
taken as 100, therefore only differential effects 
enter as errors in the pattern The peaks re- 


and BH 


small and near the threshold of observation, and 


normalized 


sulting from the ions B,D;H are very 
it is not unexpected to observe a pattern coefli- 
cient of zero for these peaks 

From table 5, the mass spectra of C,D, and 
C,H, are also quite similar. The sensitivities of 
the most abundant ion in both compounds agree 


is the approximate 


to about 3 percent, which 
uncertainty of the measurements. The similarity 
of the spectra extends to the doubly charged ions 


The amount of 


and the metastable transitions 
C,D;H present is too smail to permit conclusions 


concerning the dissociation of this molecule com- 


pared to ethane-d The abundance of the 
molecule ion in ethane permits an estimate of the 
of C,D,H equality in 


sensitivity of the two ethanes The 


amount assuming only 
35 peak 
indicates approximately 1.5 percent C,D;H in the 
The dissimilarities of C,H, and B.H, 


mixture 


spectra are already noted [3 The same observa- 


tions are applicable to C,D, and B,D,, name'y, ; 

low abundance of the molecule ion and the BD) 
ion in the latter molecule compared to the anal. 
ogous ions in ethane-d,. In addition, a metuasta! 
transition of the B.D,* ion observed in dibo im 


has no analog in ethane-ds. 


V. References 


S. Pitzer, J. Am. Chem. Soc. 67, 1126 
R. E. Bundle, 69, 1327 (1947); and A. B 
66, 747 (1947 

C. Price, J. Chem. Phys. 15, 614 (1947 
Mohler, J. Am. Che 


{2} W 
[3] V. H. Dibeler and F. L 
70, 987 (1948 

{4} \ H Dibeler, F 
R. M 
RP2010 
Turkevich, L. Friedman, E 
Wrightson, J. Am. Chem. Soc. 70, 2638 (1948 

i6] A. E. Finholt, A. C. Bond, Jr., and H. IL. Schles 
J. Am. Chem. Soc. 69, 1201 (1947 

[7] Metal Hydrides, Incorporated, Beverly, Mass 

IS] J. A Hipple, R. | Fox, and I I Condor 
Rev. 69, 347 (1946 

[9] J. A. Hipple, Phys. Rev. 57, 350 (1940). 

{10} E. G. Bloom, F. L. Mohler, J. H 
Wise, J. Research NBS 40, 437 (1948) RPISSS 

fil] M. W. Evans, N. Bauer, and J. Y. Beach, J. | 
Phys. 14, 7OL (1946 


Mohler, L. Williams: 
Research NBS 43, 97 


Reese, J 


Solomon, and | 


Lengel, and 


WaAsHINGtTON, October 21, 1949. 


Journal of Researcs 





y. S. Department of Commerce 
National Bureau of Standards 


Research Paper RP2096 
Volume 44, May 1950 


Part of the Journal of Research of the National Bureau of Standards 





An Improved Apparatus for the Determination of 
Liquidus Temperatures and Rates of Crystal 


Growth in 


Glasses’ 


By Oscar H. Grauer and Edgar H. Hamilton 


The glass is the 


quidus temperature of a 


can coexist between the molten glass and 


the 


highest temperature at which equilibrium 


primary ervstalline phase An improved 


temperature-gradient apparatus is described for the rapid determination of liquidus tem 


peratures and rates of crystal growth in glasses 


ifined i 


ments col 


in previous gradient methods 


in the range from 950° to 1.225°C 


Rates of crystal 
obtained bv this met! 


ogist of data 


I. Introduction 


the 


wortant limitation to the development of new 


vstallization or devitrification is most 


ses. The significance of this fact has been 
nized by many investigators including Little- 
in 1951 1940, 


jhasized the need for further study of devitri- 


and Preston [2] in who 


tion in glass Any such study must take into 


nt two factors: (1) the liquidus temperature 


iXimum temperature at which equilibrium 


ts between the class and the primary crystal- 


- and (2) the rates of crystal crowth at 


hase 


s temperatures. Such data are an important 


t only in the development of new glasses but 


the manufacturing processes of melting, 


and annealing 


of the data that could be used in devitri- 
as phase equilibria 


Hall 


liquidus tempera- 


studies are @xpressed 


compiled by and 


the 


such as thos 


is 

|, which indicate 
or various compositions but which present 
The 


obtaining the liquidus tem- 


rmation as to rates ol crystal growth 
al method for 


re is the quenching method, which was used 


The procedure is simple and rapid 
were determined to 
growth were found to be linear 


od is indicate 


In this method the use of small glass frag 


the cells of a platinum alloy holder largely eliminates the errors inherent 


] iquidus temperatures 
a precision of 4 degrees ( 


The usefulness to the glass tech: 


dl 


nt the Geophysical Laboratory ol the ( arnegie 
of Washington 


In this method a small powdered 


Institution in phase studies of 


silicate systems. 
sample of glass or crystals, encased in a suitable 
container, usually a platinum foil envelope, is 
heated at a given temperature until equilibrium 
is attained and then quenched in air or some suit- 
able liquid. The sample is examined with the 


petrographic microscope for crystals, and the 


process repeated at other temperatures until the 
highest temperature is found at which erystals 
persist 

In 1914 Tammann |4] determined the maximum 
rate of crystal growth of glass rods that had been 


Zschim- 


determined rates of crystal 


heated in a known temperature gradient 
Dietzell [5] 


growth by placing a piece of glass on a platinum 


mer and 


plate in a furnace and maintaining a given tem 
After the 


Wiis Ser 


perature for a specific time interval 


sample was removed from the furnace, it 
tioned and the crystal size measured with a m 


his 


for each temperature at which data wer 


crometel! microscope procedure Was re 


peated 


desired, and from a graph of temperature and 


ervstal size the temperature of maximum devit 
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fication was determined. Tabata [6] used a similar 
procedure to indicate relative rates of crystalliza- 
tion Swift 


modification for determining the rates of growth 


[7] used the same method with a 


at temperatures immediately below the liquidus 
temperature. For these cases, the sample was 
placed over a small hole in a platinum foil, which 
permitted direct observation with a microscope 
after each heat treatment 

Parmelee and Monack [8] set up criteria for 
evaluating devitrification in glasses based on 
prescribed cooling cycles for small experimental 
melts 

A gradient method to determine liquidus tem- 
perature was developed by Silverman [9]. In 
this method a sample of powdered glass was 
placed in a platinum-alloy boat, which was in- 
serted in a furnace in which a known temperature 
gradient was maintained Preston [2] used the 
same method to determine both liquidus tem- 
peratures and rates of crystal growth 

In a previous investigation at this Bureau [10] 
the Silverman method was found to give good 
results except for those glasses that were very 
fluid at temperatures in the region of their liquidus 
temperatures. Such glasses had a tendency to 
flow when the boat was removed from the fur- 
nace, frequently resulting in erroneous values 

In the present investigation an improved ap- 
paratus has been developed that permits a rapid 
determination of liquidus temperatures and rates 


of ery stal vrowt h 
II. Apparatus and Procedure 


The apparatus used for holding the glass sample 
is illustrated in figure 1. It consisted of an in- 
piece of  platinum-—10- 


serted channel-shaped 


percent-rhodium alloy, 5': in. long by ': in. wide 


cells for 


glass were spaced along the length of the holder 


Open containing small fragments of 


The cells were tapered from a 2-mm-diameter 
opening at the top to a l-mm-diameter opening 
at the bottom and extended 1 mm below the 
surface. At one end of the channel the cells were 


spaced in two parallel rows ‘i. in. apart. In 


each of the parallel rows the cells were separated 
by ', in. from center to center, but since one 
row offset the other by ‘\« in., the distance along 
the channel between centers of succeeding cells 


} 


was 'j,in. This interval was found to correspond 
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to a temperature difference of about 5 
in the temperature gradient. At the ot! 
of the apparatus the cells were spaced 
apart between centers in a single row alo 
length of the strip, corresponding to a 20 
temperature difference. In actual operati 


holder was so positioned in the tempi 


gradient that liquidus 


determined within the region of the closely 


temperatures col 


cells of the parallel rows and measurem 
crystal sizes within the cells of the single 
The horizontal tube furnace consisted 
concentric, spirally grooved alundum cores 
sulated from each other and encased in a 
shell packed 


inner core was wound with platinum—20-pe: 


with diatomaceous earth 
rhodium wire along a central zone for a dista: 
6 in. and the outer core with a nickel-chror 
alloy wire along 20 in. of its length 

A fairly linear temperature gradient was 
tained by adjusting the current in both winding 
with small variable transformers. The maxin 
temperature of the furnace was maintained « 
stant to within +2 deg. C by means of a 
mercial-type electronic controller operating 
a platinum—platinum-10-percent-rhodium ther 


As a resu 


stable gradient was maintained in the remain 


couple located in the hot zone 


of the furnace 
Temperatures in the gradient were measu 


with a movable platinum—platinum-10-perce 


rhodium thermocouple in conjunction wit! 


potentiometer The thermocouple was en ased 
a double-bore clay insulating tube except for 
last half inch at the hot junction. 

In making liquidus temperature determinati 
and measurements of rates of crystal growth 
glass, the following procedure was used \ 
piece of glass or of glass-crystal mixture, al 
cm in each dimension, was crushed in a har: 
steel mortar, and the fragments were transfer! 
to a sheet of black paper or black glass to fa: 
the selection of particles of the desire 
Pieces about 2 mm long and 1 mm wid 
loaded into the cells of the platinum holde: 
was then inserted and positioned accurat 
the furnace 
gradient furnace for ihe desired time, usually 
A temp: 
probe at ‘-in. intervals was then made 


The sample holder was left 
to 60 min. for experimental glasses. 
gradient with the movable thermocouple, 2 
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GLASS SAMPLE THERMOCOUPLE 


TUBE 


FIGURE 1 Platinu 
migh A 
vable 


ew; c, Section thr 


was removed immediately (air-quenched) 

om the furnace by means of a long wire hook. 
since the cells confined the glass, ho displace- 
due to flow occurred when the holder was 
moved rapidly the Surface 
on prevented the glass from flowing out of the 


from furnace. 


opening. After air quenching, the holder 
placed on the microscope stage, and the 
nts of the cells were examined with a petro- 

microscope The opening in the bottom 
h cell permitted direct examination of the 
ed sample using transmitted light, thereby 
ng the necessity of breaking away the glass 
The cell that 
s was considered to locate the position in 


the holder. last contained 


mperature gradient corresponding to the 


us Temperatures and Rates of Crystal Growth 


m auoy spec 


Ao 


thorn 


f 





men holder 


view (t showing positior 


ocouplk 


The of the 


cell was obtained by linear interpolation between 


liquidus temperature temperature 
the two nearest points of the temperature probe 

In general, each fragment of glass tended to take 
the shape of the platinum cell when the tempera- 
ture was above the softening point and to assume 
a somewhat lenticular appearance. For photo- 
micrography it was desirable to avoid the optical 
the 
This was effected by placing 


distortion and uneven illumination due to 


curved surfaces. 
cover glasses at both surfaces of the sample and 
interposing a drop or two of a liquid having the 


A 


cover glass with an appropriate liquid sufficed for 


same index of refraction as the glass sample. 


the visual measurement of crystal sizes 


In determining the crystal size, Preston’s 
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method of making measurements of the semimajor 


axis of the erystal or of the radius of spherulitic 
crystals was used. Corrections were made for the 


foreshortening effect of those crystals that did not 


lie in the plane perpendicular to the optical axis of 


the microscope. A microscope, containing a filar 
micrometer eyepiece, was focussed sharply on one 
end of the crystal or of the radius of the spherulite 
The crosshair was brought into coincidence with 
the image, and the readings of the micrometer and 
the mi roscope fine adjustment scale were noted: 
the microscope was refocussed on the other end 
of the crystal or on the radiant point of the 
spherulite; the crosshair was repositioned ; the new 


both The 


horizontal projection of the crystal or of the radius 


si ales were observed. 


readings of 


was obtained from the difference of the micromete! 


scale readings, and the vertical projection from 


the product of the difference of the microscope 
index of refraction of the 
the the 


scale readings and the 


From these data true size of 


glass 
crystal was determined by the geometrical relation- 
ship as indicated in figure 2. An average of four 
different 
The reproducibility of the measure- 


measurements of spherulites was used 
for the data 


ments was about 2 percent 





FIGURI m of spherulitic rad 


( y 4°+( BN 
t amd N 


Measurements by Preston [2] indicated that the 
rate of erystal growth in glass was linear with 
time only during the first stages. A more recent 
7] showed that the rate of 
until the 


crystallization zones from opposite surfaces met 


investigation by Swift 


erystal growth in glass was constant 


In the present investigation the linearity of crystal 


growth was also observed from measurements 


made on various types of glasses 
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Thus, both the rates of crystal groy 


various temperatures and the liquidus temperat 


could be obtained after only one trea 


Or 
I 


However, in those cases where the rate of crys 


growth was very rapid, a short treatment « 
10 min was desirable so that the crystals 
region of maximum devitrification would not 
fere with each other and depart from linea 
their rates of growth. In these few insta 
separate determination was made of the lic 


temperature 


III. Results and Discussion 


1. Liquidus Temperature Determinations 


Liquidus temperature determinations wer 
for five different glasses using the classical qu: 
ing method, the Silverman method, and th 
The results are listed in table 1 ar 
of the table 2 


values of duplicate determinations with thi 


method. 
compositions glasses in 
method show less seatter than those with th 
the 
tained by the quenching method as the stan 


verman method. Considering values 


Ty 


+ 


a! 


for comparison, the average deviation of the va 


obtained by the new method is 3.6 deg C as 


pared to 7.2 deg i 


for the Silverman meth 


The greatest deviation between duplicate det 


minations for the new method was 6 dee C 


pared to 23.4 deg C for the Silverman metho« 


ature determination 


TABLE | 


TABLE 2 Com positions 


BaO Bed 8) 


LayO) Nay 
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mong the glasses that are very fluid at their 
temperatures that the results obtained by 
rman method are erratic, due to the flow 
olten glass in the boat upon removal from 
nce 


new method flow is practically eliminated 


y 
° 


ning the specimens to the small cells of the 
Due to the small mass of the sample and 


8 


nate thermal contact with the platinum 





ll, equilibrium is more quickly attained, 
| cooling of the sample on removal approaches 


2 


enching methed in its rapidity. 


2. Rates of Crystal Growth 


RYSTAL SIZE , SEMI,MAJOR AXIS IN MICRONS 


} gives the values of crystal size and tem- 


c 


for glass E 1105, for heat treatments of 
and 60 min It indicates that the maxi- 


m size for each holding period occurs at about 








same temperature, 895° C. Up to 1 hr the | 
60 


tes of crystal growth (at a given temperature), 
TIME , MINUTES 





FIGURI Effect of duration of treatment at constant tem- 


perature on crystal sir 





Oo—O = 3 min 
x——X * 26 min 
+——+ = 60min 


© 
° 
°o 








| | | 
300 400 500 
EMI-MAJOR AXIS IN MICRONS 





TEMPERATURE ,°C 


@ 
wn 
o 


j 


Effect of temperature on crystal sir 


ments of 13, 26, and 60 min 


are the slopes of the curves in fig. 4, were 
0 be linear. When rates of crystal growth 
ted against temperature, the curve in fig. 5 
ned. In this curve the rate of crystal 
is zero at the liquidus temperature, and 
iximum rate occurs at about 895° C. 








6 and 7 are photomicrographs of devi- F \ 1 1 L 1 





ystals that were grown in glass E 1105 2 4 6 ~ 10 i2 
\a0. 20¢ / -CaO. 10° / - Sil ),, 70% ) at 860° C for RATE OF CRYSTAL GROWTH, yw /MINUTE 


al 26 min. respectively. Fiat RE 5 Rate of ystal growth reraus te mperat ire 
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Ficure 7. Photomicrograph 


of devitrite crystala << 200 
after 26 min treatment at 
SAO) bad focus at center of 


crystals 


Ficure 6. Photon 


of devitrite crysta 


alter 13 min trea 


sao ¢ 





TI 
obta 
ents 
the 
ng 
mick 
some 
brou 
focu 


grow 


Photomicrogr aph 
devitrite crystals 
> with microscope 


/ at ends of crystals 


The photomicrograph shown in figure 7 was 
obtained with the microscope focused on the 
enters of the spherulites. It can be seen that 
the end of the spherulites are out of focus, indicat- 
ng varving degrees of inclination. The photo- 
micrograph shown in figure 8 was obtained when 
some of the ends of the same spherulites had been 
brought into focus; this time the centers are out of 
focus. Values of maximum 


growth for two glasses are given in table 3. 


crystal 
Other 
leterminations of crystal growth as an indicator 


rates of 


Rates of crystal growth for glass forming and 


nonglass-forming materials 


Maximum rate 
of crystal 
growth 


Material 


Glass FE 1105 

Glass E 1222 * (devitr 
fies readily 

Picric acid » (nonglass 860), O00 
forming 

Phosphorus ¢ 60, 000, 000 


glass-forming 


persion, high-index experimental glass. 
jawlensky, Z. physik. Chem. 27, 585 (1898). 
Gernez, Compt. rend. 95, 1278 (1882 


mum 


of devitrification have been made with this ap- 
paratus. From these results, suitable time-tem- 
perature schedules have been derived for the 
slumping of large interferometer plates and other 
large disks of various types of optical glasses 
without the usual devitrification. 


IV. Viscosity and Crystal Growth 


Various attempts [7, 11, 12, 13, 14] have been 
made to find a relationship between the rate of 
Below the 


liquidus temperature, the stable state of a system 


crystal growth and the temperature. 


consists of crystals or a mixture of glass and 
crystals. The glassy phase at ordinary tempera- 
tures is unstable; even at a temperature below 
the softening point of the glass, the viscosity may 
be sufficiently low so that crystals may grow at a 
measurable rate. 

In reviewing the kinetics of crystallization, the 
familiar hump-shaped curve of rate of crystal 
growth and temperature can be considered as the 
result of two opposing forces: a driving force that 
increases as we depart from the liquidus tempera- 
ture and a resisting force related to the viscosity, 
which increases exponentially and plays the domi- 
nant role as the temperature is diminished. The 
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rate of crystal growth has been shown to be in- 
versely proportional to the viscosity: by Frenkel 
[11] on the basis of the classical activation theory, 
by Richards [12] using the concept of thermal 
relaxation time, and by Leontjewa [13], empiri- 
cally in a series of Na,O-SiO, glasses. The latter 


used the equation, 


K nt Ky log n, 


where v is linear speed of crystallization, 7 is the 
viscosity, and A and Ay are constants, Ay being 
very small 

Swift 


Preston, 


{9] used the semiempirical equation of 


where 2? is the rate of crystal growth, 7), is the 
liquidus temperature, 7 is the temperature during 
holding period, 7 is the viscosity, and K is a con- 
stant. A plot of this equation fitted one set of 
his data fairly well but not others 

From a consideration of the difference in vibra- 
tional energy between condensation and melting 
at a plane glass-crystal boundary, Cox and Kirby 
[14] derived an equation for the rate of crystal 


growth of the form, 


where / 


the ions in the glass mobile, /, 


is the vibrational energy needed to make 
the energy that is 


necessary to make the ions in the crystal mobile, 


r the ratio of the surface densities of crystal and 
A and K 


While this equation gave a fair fit 


glass, 7 the absolute temperature, and 
are constants 
for most of their data, there was poor agreement 
at temperatures immediately below the liquidus 
temperature for the two borosilicate-type classes 
that they used 

A comparison of the rates of crystal growth in 
glass-forming and  nonglass-forming materials 
table 3) shows the tremendous differences in their 
rates of growth. As the rate of crystal growth 
increases it becomes increasingly difficult to pro- 


duce glass 


V. Summary 


An improved apparatus for the rapi 
mination of liquidus temperatures and 
crystal growth has been developed. 

The simplified 
apparatus makes it possible for an inexp: 


technique required w 


operator to determine quickly data that form 
required much time and skill 

A rapid means is provided the glass technolog 
for assessing the effects of changes in composit 
upon the devitrification characteristics of glasses 


thanks to William 
Eubank for his aid in determining liquidus t 


The authors express 


peratures by the quenching method, to Charles 
Saylor for aid in preparing the photomicrogra; 
and to Clarence H. Hahner for helpful suggesti: 
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Infrared Spectra of Bromochloromethane, Dibromo- 
methane, Tribromochloromethane, and Tetrabromo- 
methane 


By Earle K. Plyler, W. Harold Smith, and N. Acquista 


The infrared spectra of bromochloromethane, dibromomethane, tribromochloromethane, 


and tetrabromomethane have been measured from 2 to 36 microns 


of other workers it 


have been observed 


2 By the use of the results 


taman spectra it has heen possible to classify all the strong bands that 


Many of the weaker bands were classified as combinations and over- 


tones Oniv a few of the bands of tetrabromomethane were observed, on account of the 


breaking down of the compound in solutions of carbon disulfide and carbon tetrachloride 


The intense bands of dibromomethane and tetrabromomethane were measured in the vapor 


state 


The infrared absorption bands of these compounds had not previously been measured 


over an extended range of wavelengths, and these measurements were undertaken to deter- 


mine the positions of weak bands so that a more complete classification of the spectra of 


these molecules could be made 


I. Introduction 


Other workers have measured the vibrational 
inds of a number of substituted methanes in the 
ired region, and also the frequencies have been 
termined by Raman spectra. From theory it 
been determined that molecules of the type of 
show two active and two inactive funda- 
tals Where interaction occurs the inactive 
encies also appear in the infrared spectrum 
n the substituted atoms are of more than one 
es, the degeneracy is removed and_ nine 
lamentals are present in the spectrum. In a 


ile of the type # % | 


v removed and six fundamentals are found in 


, the degeneracy is only 


spectrum 

four substituted methanes, bromochloro- 
ine, dibromomethane, tribromochloro- 
ine, and carbon tetrabromide, represent the 
hnree ty pes of molecules that vive nine, SIX, and 
etive fundamentals in the infrared spectrum. 
nfrared bands of these compounds had not 
uusly been measured over an extended range 
ivelengths, and these measurements were 


taken to determine the positions of the weak 
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bands so that the combination bands could be 


classified 
II. Experimental Observations 


A Perkin-Elmer spectrometer was used for all 


measurements with lithium fluoride, potassium 
bromide, sodium chloride, and thallium bromide- 
iodide prisms to cover the wavelength range of 2 to 
36 u. In the spectral region of 2 to 24 u, cell 
thicknesses of 0.2, 0.1, and 0.05 mm were used to 
bring out most of the bands For very strong 
bands the compounds were diluted either in carbon 
tetrachloride or in carbon disulfide, or liquid films 
with a thickness of 0.01 mm or less were formed 
In the 


thallium bromide-iodide region between 24 and 


between two potassium bromide windows 


56 w it Was necessary to use cell thicknesses of 1.5 
mm to bring out the weaker absorption bands 
The method of measurement and the reducing of 
data have been described in a previous paper 1] 

Bromochloromethane, dibromomethane, and _ tri- 
bromochloromethane were obtained from the Dow 


Figures in brackets indicate the literature references 


paper 





Chemical Company and tetrabromomethane from not measured, but it is estimated that the 
the Eastman Kodak Company of the order of 0.005 mm. This estimate 

Purification consisted of treatments until an thickness was made on the basis of the per 
examination of the products indicated satisfactory absorption that was observed in a 0.1-m 


agreement with accepted physical constants. To when the compound had been diluted with a trans 


avoid any effect of oxygen, the infrared measure- parent solvent. The reason that the thin ce] y 
ments were made as quickly as possible after employed was to locate the position of the ba 
purification. accurately for the pure liquid in regions of intens 
Figure 1 shows the absorption spectra of the absorption. In figure 1 the portions of t! 

four halogenated methane derivatives in the region sorption curves represented by broken lines ind 
from 2 to 15 w. The cell thicknesses and other cate that small details could not be accurat 

conditions of the measurements are given on the determined on account of the absorption of t! 
figure and in the captions. The intense bands solvents or the absorption of atmospheric bands 
were observed by use of a cell that contained a As tribromochloromethane and tetrabromomet! 
thin layer of the material. The cell was made by ane are solids at room temperature, their spect 
placing a small quantity of the liquid on a plate were determined in solutions of carbon tetrach\ 
of potassium bromide and then pressing another ride and carbon disulfide. There was a tenden 
plate on the top without the use of a shim. The for the solutions of tetrabromomethane to dark: 
thicknesses of the cells made in this manner were on standing, and the absorption spectrum \ 


TRANSMISS/ON, PERCENT 








a 
WAVELENGTH IN MICRONS 
Fiat RE 1. Infrared absorption spectra of bromochloromethane, dibromomethane, tribromoc hloromethane, ar 


tetrabromomethane. 


rribromochloromethane and tetrabromomethane were dissolved in carbon tetrachloride and carbon disulfide The solution concentration of | 
chloromethane in carbon tetrachloride is 3.1 g/ml, and in carbon disulfide, 2.4 g/ml for the observations in 0.2-mm cell. For the insert bands the con 
is 0.14g/ml. Saturated solutions of tetrabromomethane were prepared in each solvent 
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d immediately after their preparation. 
servations were repeated several times and 
e bands were always found. Since it was 
to account for all the observed bands as 


ations, overtones, or fundamentals, it does 


n probable that any of the observed bands 


oduced by products of decomposition or 
on in the solution 
e 2 shows the long wavelength spectra of 
substituted methanes from 14 to 36 uy 
general absorption was observed in the 
a of bromochloromethane and dibromometh- 
om 28 to 36 u, but no definite bands could 
nd in this region 
ochloromethane in the region of 24 to 364, 





For the measurement of 


. 


a saturated solution in methyleyclohexane was 
used. The methylcyclohexane has a high trans- 
mittance in this region, and its spectrum shows 
only one absorption band when a cell 1.5 mm 
thick is used. 

In figures 3 and 4 are shown some of the absorp- 
tion bands of dibromomethane and tetrabromo- 
methane measured in the vapor state. The 
substances were placed in the bottom of a cell 40 
em long, which was not evacuated, and allowed 
to remain in the cell for several hours. On re- 
peating the measurements 6 hours later, it was 
found that the bands did not increase in intensity. 
It is probable that the vapors were at a saturated 


condition in the air of the cell at 25° C. 





WAVELENGTH 


Absorption spectra of bromochlo omethane, dibromomethane, 


2€ 


N MICRONS 


tribromochloromethane, and tetrabromomethane in 


the region from 14 to 36 yu 


n concentration of 


mn 
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n carbon disulfide was used for tribromochloromethane with the 0.1-mm cell 


cel 4 saturated solution in methyleyclohexane was used in the 


A concentration of 0.14 g/ml was used 


gion of 24 to 364 


505 





of W. Bacher and J. Wagner {2}. The oth 
damental bands, as determined by Raman s; 
check within a few wavenumbers with i 
measurements given in table 1, except », 

was estimated by them at 724 cm™' from R 


raABLE 1 Observed frequencies of the infrared 


hromochloromethane dibromomethane tribromoch 





thane, tetrabromomethane, and their assigmnents 
TH IN MICRONS 

Abso plion spec em of dibromomethane ra po 
loIlSp 


turated vapor (45 mm of Hg at t ) THAN DIRROMETHANE 


WAVE NUMBERS IN Cm™' 


CBr, VAPOR 


TRANSMISSION , PERCENT 


40 cm 








14 
WAVELENGTH IN MICRONS 
hands at 12.6 and 14.9 pw fo 


at room tem pe 


III. Discussion 


In table 1 are given the wavelengths and fre- 
quencies of the observed bands. These bands 
have been interpreted as fundamentals, combina- 
tions, and harmonics The numbering of the 
levels is in ace ordance with that adopted for corre- 
lation of a large number of molecules of different 
symmetries, as will be discussed in a forthcoming 
papel by Plvler and Benedict Kight of the 
infrared active fundamentals of bromochloro- 
methane have been observed, and the other, 
occurring at 226 em was outside the range of 
the thallium bromide-iodide prism. In the deter- 
mination of the combination bands for this com- 


pound the Raman value was used from the work 
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and has been observed by us in the 

spectra as a strong band at 852 em™'. 

the fundamental bands of dibromome- 

ve been observed in the infrared spectrum 
at 174 cm 

wed by Dadieu and Kohlrausch and several 

The Raman 


for the remaining eight fundamentals check 


This value has been 
orkers in Raman spectra [3]. 
th these observations, except that »; and 


found at 813 and 1,190 em", 
he infrared as compared to 723 and 1,133 


respectively, 


espectively, in the Raman spectrum as 
The other 


rved bands, except for a few weak ones, have 


d by Dadieu and Kohlrausch. 
assified on the basis of the nine funda- 

\ close similarity between the spectrum of 
bromomethane and of bromochloromethane is 
bserved. The effect of replacing the chlorine in 
omochloromethane by the heavier atom bromine 
o shift corresponding vibrations to longer wave- 
The shift is observed for all frequencies 

ept those involving primarily C—H stretching 
rations. As can be seen from table 1, the 
for the two compounds fall 
This difference in 


es of », and of x 


thin 5 em”! of each other 
venumber is almost within the limit of experi- 
ntal error for this region of the spectrum. 

The spectrum of dibromomethane was also 
termined in the vapor state for some of the 
nse bands, so that a comparison could be made 
th ‘ liquid state. The type of bands at 
could not be determined with 
For the band at 
in the liquid state there is a small 
The other 


5 and 648 em 
rism instrument available 
1,195 em™' in the vapor state 
hifted from 639 to 648 em Two of the 
ng bands located at 810 and 591 em 
finite zero branches when measured in the 
ate. These bands have shifted from 813 
9 em respectively, as measured in the 


state. The vapor band at 745 em is 


em! from the liquid band at 729 em 


v have a zero branch. The comparison of 
tions of the bands in the liquid and in the 
tates shows that there is a small change 
frequencies when measured in the two 
states, 
e six active infrared fundamentals of 
whloromethane, only three bands, », v5, 


vere observed in the infrared region that 
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was covered in these measurements. The other 
three fundamentals fall beyond the range of the 
Calculations 


thallium bromide-iodide 


show that two of these bands, Vr and ,, should 


prism. 


have frequencies at 214 and 215 em respec- 


tively The calculated values are in good agree- 
ment with the wavenumber of a single Raman 
line determined by Lecomte and coworkers [4] 
For the classification of the combination and 
harmonic bands, the Raman values 141 and 2:4 
cm'' have been used. The three fundamentals 
observed in the infrared checked within 5 em 
with the Raman values of Lecomte. The differ- 
ences between the infrared and Raman values 
may be attributed to the effect of the solvent 
carbon disulfide on the tribromochloromethane in 
producing a small shift in the bands, and also to 
the presence at 667 cm”! of the absorption band 
of atmospheric CO,, which made it difficult to 
locate accurately the position of the band at 


675 cm The two absorption bands at 735 and 


‘ 
could not be accounted for on the basis 


9335 em 
of the fundamentals. 

In comparing the fundamentals of tribromo- 
chloromethane and tetrabromomethane, it is 
seen that the substitution of the bromine atom 
for the chlorine atom reduces the number of 
fundamentals from six to four. Only the two 
threefold degenerate Vibrations, vgs and ry, are 
infrared active. The remaining two appear in 
bands, vies, with 


combination. One of these 


frequency 669 em was observed in the rocksalt 


region. This is 4 em! greater than the Raman 


value as determined by Dadieu and Kohlrausch 
[5] This difference can probably be attributed 
to error of observation caused by the presence at 
667 em=' of the absorption band of atmospheric 
Cc ). The other bands in the observed spectrum, 
which are of low intensity, can be interpreted on 
the basis of the four fundamentals. Three of the 
fundamentals listed in table 1 are the Raman 
values determined by Dadieu and Kohlrausch 
The classification of the bands as given in table | 
accounts for most of the observed bands. The 
band observed at 729 em=! in dibromomethane 
was not classified and may have arisen from a small 
amount of bromochloromethane that was present 
as an impurity. There was observed an intense 


band in this region for bromochloromethane 


When the vapor ol dibromomethane Was necias- 


ured, the ratio of intensities of the bands at 735 





were much different than in 
l 


em ' and 813 em 
the liquid, and also a band at 1,233 em™! appeared 
This is near the location of another intense band 
1 


of bromochloromethane at 1,225 ecm”! as measured 


in the liquid state. Another sample of dibromo- 
methane was measured that also contained these 
bands. According to mass spectrograph deter- 
minations it contained about 1 percent of chlorine 
It is probable that a small amount of bromo- 


chloromethane is present in the dibromomethane 


and produces the bands that do not fall in the 


classification scheme. On further purification of 


dibromomethane it was found that the band 


in the region of 13.6 w (729 em disappeared 
when the compound was measured in cells 0.1 
mim thick. 

Without a study of the rotational structure of 
these bands the molecular constants cannot be 
obtained. Unpublished observations made with 
a different instrument show resolution of some 


of the near infrared bands of dibromomethane 


and bromochloromethane into a series 


approximately 1.5 cm”! apart. The 


interval is less for dibromomethane 


bromochloromethane 


The authors thank W. S. Benedict 
Bureau for his discussions on the classific 


the bands of these compounds. 
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Mass Spectrometric Analysis of a Standard Sample of 
Carburetted Water-Gas by Laboratories Cooperating 
With the American Society for Testing Materials 


By Martin Shepherd ' 


The mass spectrometer was used for the analysis of a standard sample of carburetted 


water-gas by laboratories cooperating with Subcommittee VII of Committee D-3 of the 


American Society for Testing Materials. 


The results of the cooperative analysis show the 


reproducibility and, in certain respects, the accuracy of this powerful new apparatus for gas 


analysis 


lytical data were compared with the known values. 


The heating value and the specific gravity of the sample calculated from the ana- 


Some very creditable work is reported, 


but the need for improvement and standardization with respect to the determination of hydro- 


gen, carbon monoxide and nitrogen is evident 


I. Introduction 


report is the fifth of a series of cooperative 
es of standard gas samples * conducted to 
sh basie information for the development of 
ndard methods for the analysis of fuel gases. 
development of these standards is a_ task 


3-VIIL of the ASTM, 


the method of this development has been 


signed to Subcommittee D 


d in four previous reports 


) 


7 laboratories that cooperated in the 
work were widely distributed geographi- 
ut most of them were associated with the 
im industry. There was some representa- 
mm the chemical industry and Government 


all the 


the operating procedures prescribed in the 


ories. In general laboratories fol- 


s furnished with the instruments. Depar- 


\ ty i i ( 


tures from the prescribed course were few and of 
such nature that no significant change was ex- 
pected or observed. The greatest variation was 
in the method of computation. Thus, with respect 
both to apparatus and methods of operation, this 
series of cooperative analyses was conducted with 
a preliminary standardization that was unofficial 
but none the less real; and, while this may or 
may not have been a factor in the accuracy 
achieved, it must have affected the over-all repro- 
ducibility. The need for further standardization, 


particularly in calculating methods, is evident 


II. Standard Carburetted Water-Gas 
Sample ASTM D-3-VII-5 


The preparation of the standard sample of th 
ASTM 


coopera Ve 


carburetted water-gas type, identified as 
D-3-VII-1, 
analysis by volumetric chemical methods, has been 
RPI704 
That account will serve to es- 
ASTM 


which was used in the 


deseribed in detail in’ Research Paper 
see footnote 34 
tablish the history of the present sample 


D-3-VITI-5, 


spectrometet! 


Which was analyzed by the 
since the No. 1 and No 


whi h had by en 


5 samples 
were identical, except for oxvgen 
mostly removed from sample 5. This informa- 
tion was not disclosed when sample 5 was issued 


and accordingly it Was analyzed as a blind sample 


Mass Spectrometric Analysis of Carburetted Water-Gas 509 





In general, it did not find its way into the same 
laboratories that had performed the chemical 
analyses, although, in a few cases, this did happen 

Sample ASTM D-3—VII-5 was issued in the 
same type of eylinder as that used for the sample 
1, and the instructions for transferring it to the 
spectromete! Without contamination were essen- 
tially the same. This information will, accord- 
ingly, not be repeated here 

Chemical and physical methods were thus put 
nto direct competition. The part of the story 
concerned with the analysis by the mass spec- 
trometer is given in this paper. A comparison of 


the two methods will be given in another paper. 


III. Analytical Results 


All of the analytical data submitted have been 
tabulated, together w ith the av erage values derived 
from each laboratory's analyses; but the study of 
these data for a long time would not serve to 
reveal what may be seen at a glance when these 
same data are presented in a series of frequency- 
distribution plots. These plots, which amount to 
actual pictures of the analytical results, will enable 
the reader to arrive at the obvious conclusions. 

In the plots, each circle represents a value de- 
rived from a single determination of the substance 
indicated in the legend. The abscissas are values 
derived from the analyses; and these values are 


plotted equidistant on the ordinates and so indi- 


cate the frequency with which these values occur. 
For example, the frequency -listribution plot for 
carbon dioxide (fig. 1) shows that 4 determinations 
gave the value 4.0 percent, & determinations gave 
1.2 percent, 10 gave 4.3, 18 


t.1 percent, 4 gave 
gave 4.4, and so on The plots are divided into 
two sections, one of which shows the laboratory 
averages, the other all of the determinations 
Thus a laboratory to laboratory comparison may 
be made, but the study of accuracy and repro- 
ducibility may be referred to the plot giving all 
determinations. This system is carried through- 
out the group of plots, and with this in mind, 
they can be studied The analy ses for each com- 
ponent will be considered 

Carbon Dioride (tig. 1 


tions for carbon dioxide are given in figure 1. The 


All of the determina- 


greatest frequency is indicated at 4.5, in agree- 


ment with the arithmetical average (determined 


510 


4 


Fiaure 1 


hon dio 


with four high values dropped), 4.49 +0.14 


age deviation), and with the median, 4.5. Th 
values agree with the most probable valu 
- 0.04, previously determined (RP1704). The la 
oratory average is 4.49 +0.17—or, with one hig 
value dropped, 4.46 +0.14. Accuracy and 
ducibility are thus satisfactory 

Oxygen (fig. 2 Oxygen occurred in this sam| 
in an amount less than 0.05 percent, so tl 
rounded off values indicate its absence in 66 p 
After dropping 


three high values, the arithmetical averag 


cent of the determinations 
comes 0.03 Highest frequency and media 
average (withou 
Most of the di 
minations of this gas are satisfactory 
Carbon Monoride fir. 3). The spectro! 


was not altogether comfortable with respt 


found at 0 Laboratory 


three high rogues) is 0.08 


this determination some analy sts cal 
from the over-populated though most set 
oS peak ; a few others used the less sensit 
peak In one case CO was completely cor 


with N, 


Gaussian pattern 


The plot shows no adherence t 


The arithmetical mean 
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bv the far flung values and disregard for the Gaus- 
sian shape. The mean, computed from values 
between 33.0 and 37.8, inclusive, is 35.1 0.8. 
The median is 35.1 These values agree well 
with the chemical determination of 35.0, even 


though the derivation was from scattered values, 





Improvement is needed, and, indeed, some has 
been achieved since these data were taken Lab- 
oratory averages range from 31.9 to 38.6; their 
average 1s 35.2 1.1] 

Methane (fig. 5 While some difficulty was ex- 
perienced with carbon monoxide and hydrogen in 





this sample (and more difficulty will occur with 
nitrogen), the spectrometer is at home in dealing 
with the hydrocarbon fractions of this gas. The 
change is evident in the methane plot. The 
mean (with three low values dropped) is 8.0 +0.4; 
the median is 8.0; and the highest frequency 
occurs at 8.1 Laboratory averages range from 
6.2 to 9.3; their average is 7.9 +0.4. In general, 
the values are satisfactorily clustered, and a Gaus- 
sian pattern is now evident. 

Ethane (fig. 6). A clean-cut analysis for ethane 
is shown in this plot. The highest frequency 
occurs at 3.3. The mean is 3.3+0.15, and the 


median is 3.3. Laboratory averages range from 2.8 
to 3.5, with an average of 3.3+0.15. The 
answer seems to be 3.3. 

Ethylene (fig.7). Determination of the ethylene 
was not achieved in the decisive manner shown in 
the determination of the ethane In general, 
this is true of the unsaturated hydrocarbons. 
Nevertheless, the determination for ethylene was 





satisfactory. The plot shows a block-like structure 











°o aaah ddd bidder hh hihi bri titiditititiitiiitiiilriiitttt tt tt tt 


of values centered about 13 percent The mean is 


o 


12.9+-0.4, and the median is also 12.9 
PERCENT OXYGEN Other hydrocarbons reported in small amounts 
are shown at the right of the plot for ethylene. 
There are not many, but out of 138 analyses, 
propane was reported 15 times, propylene 18 times, 
from values between 28.4 and 32.3, inelu- n-butane 6 times, butenes 6 times, “C,”’ three 
(0.3 +0.9. which is in reasonable agree- times, C,+ heavier hydrocarbons 12 times, and 

vith the chemical determination of 29.41 acetylene 3 times 

even though its derivation included a wide Nitrogen (fig. 8). The values for nitrogen are 
of values No highest frequency is ap- sufficiently diverse to be statistically inconvenient 
the median is 30.35. The laboratory from the analyst’s viewpoint. The laboratory 
s range from 3.5 (dropping one averages vary from 4.3 to 13.8 (dropping eight 
hen ave rage 5 Need for im- values at 36 . with an average of 6.4+ 0.1 The 
ent is evident mean for all determinations ts 5.9+0.8, if values 
gen (he. 4 The spectrometer Is not al- from 4.2 to 8.2, inclusive, are accepted The 


happy with this determination, as shown pattern is not molded along ( lassi a | lines The 
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5 aaonenene that calculated specific gravity is a direct over-all 

AVERAGES measurement of the accuracy of the analysis. 
Calculated Heating Value (fig. 10). The plot 

again shows no decisive Gaussian pattern, The 


mean of values from 529 to 566 is 550+7, and the 





median is 551. Laboratory average is 547-+8, 


dropping the value 458. The measured value 
is 5498 


IV. The Tabulated Data 





No attempt has been made to set forth all of 
the statistical inferences that can be derived from 
the data given in this report. The whole object 


of the report is to present a general integrated 





picture of the real state of analysis of this kind 
without, for the moment, distracting particulars. 
This is not to say that a statistical study will not 
prove of equal or greater value, but it should be 





presented separately. To tempt someone to do 
this, the data are given in table 1 at the end of 
this paper. 

The table is self-explanatory except for the 
column “Operator-Spectrogram-Calculator Key.” 
The experiment was designed to show the effect 
of these three factors—the operator of the spec- 
trometer, the instrument itself as indicated by 
the spectrogram it produced, and the calculator 
who derived the final answer from the spectro- 


PP +++++++ +4 gram. Each laboratory was asked to produce 
i A. 


four spectrograms of the standard sample, using 











3.0 35 
PERCENT ETHANE two instrument operators and two calculators. 


The complete pattern for each laboratory should 
Figure 6. Frequency-distribution . 


plot for percentage of ethane. have been 


Ol-S! 
; -" OLS! 
an is 5.8, dropping eight values at 36.4 


Ol-S2 
0.06. Certamly this determination must be 02-83 
oved for this type of gas. O2-S3 


_ ’ ’ . . : . (02-84 
lated Specifie Gravity (fig. 9). Specific 


( 
( 
( 
median is nearer the best chemical value of O1-S82-C* 
( 
( 
( 
( 


02-S4 


v calculated from the mass spectrometric 
es show rogues to both right and left of where O means operator, 5 designates spectrogram, 
The mean of all values between 0.626 and and C, the calculator. The O-S-C key of the 
s 0.650+0.007; the median is 0.651. The table is thus explained. The table discloses only 
an pattern is distorted, with maxima sug- seven complete patterns, while others are partially 
at 0.644 and about 0.651. Laboratory complete, or not clearly defined. 
s range from 0.617 to 0.676,with an average Before passing these data along to competent 
i8+0.009. The actual measured value for statistical talent, a few observations may be 
nple is 0.6475.4 It should be remembered briefly made. If the complete patterns for 


ired t Carrol Creitz at this Bureau I he § The he 
ribed in NBS Miscellaneous Publication M177 (1947 his Burean 


ravity was mea 
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laboratories 2, 4,8, 10, 11, 17, and 20 are examined 
in the instance of calculated specific gravity, 
which is a good over-all measure of the success 
of the analysis involving all components deter- 
mined, these facts are evident: 

1) The differences 
using the same spectrogram in the same laboratory 
small. The 


seven laboratories (in 


between two calculators 


were average differences for the 
the order given above), 
expressed in thousandths, are: 1 1/2, 0, 1, 1/4, 
4 1/4, 1/4, 1 1/4; the average of all the differences 
is 1.2 This calculation was based on differences 
between O1-S1-Cl and O1-S1-C2, O1-S2-Cl 
and O1-S2-C2, O2-S3-C1 and O2-S3-—C2, O2-S4 
Cl and O2-S4-C2 

2) The differences between two operators using 
the same instrument in the same laboratory were 
small. The average differences for the seven 
laboratories are: 1 1/2, 1, 1/2, 1, 3 1/2, 1/4, 1 1/4. 
The average of all differences is 1. This caleu- 
lation was based on the average of all Ol measure- 
ments minus the average of all O? measurements 
for each laboratory 

3) The differences between spectrograms were 
larger than between operators and calculators. 
For the seven laboratories, the average differences 
are 1 1/2, 1, 4, 4 1/4, 4 1/4, 1 3/4, 1 1/4. The 
average of all differences is 2.6. (This calculation 
was based on the differences between O1-S1-Cl 
and O1l-S2-Cl1, O1l-S1—C2 and O1-—S2-C2, 02 
S3-Cl1 and O2-S4—-C1, O2-S3-—C2 and O2-S4-—C2 
See O-S—-C key above 

These 


previously in the analysis of a natural gas by the 


observations with those made 


agree 
mass spectrometer (see footnote 3) 

In a paper given December 27, 1949 to the 109th 
annual meeting of the American Statistical Asso- 
ciation, W. J. Youden remarked: “It is a fact of 
experience that a set of measurements made by 
different operators alt different times or in different 
localities is subject to greater variation than a set 
of measurements made by one operator using the 
same apparatus on the same day.’ The data 
presented here offer no exception to this fact of 
Thus, 


experience for example, the mean of all 


values of calculated specific gravity (be we, 
0.626 and 0.677) is 0.650 +0.007; the mean 
L 0.009: b 
reproducibility of each laboratory varies fro 
0.0005 to 0.004, with an average of 0.0017 
times better than the reproducibility ol 


laboratory averages is 0.648 


from the laboratory averages. 

The factors controlling interlaboratory 
ences will probably never be completely evaluat 
One important factor in the present cas 
difference in methods of solving spectrogran 
not enough is known about this whole 
Particular attention is being given this ph 
mass spectrometric analysis in the preparat 
ASTM tentative standards. 


V. Summary 


As often happens, a group of values that 
frequency plot devoid of the classical Gauss 
pattern and are, from the analyst’s viewpoi 
inconveniently various, nevertheless have yield 
a mean value very close to the most probab 
the actual one. This is always a comfort if or 
has an abundance of data from as many sources 
are necessary; but the problem remains to d 
mine the accuracy of a single source in the abser 
of sufficient data from many sources. Thus 
spectrometric determination of hydrogen, ar 
carbon monoxide with nitrogen, has not 
altogether satisfactory. However, the mass sp 
trometer is quite at home in the analysis of 
hydrocarbon fraction of such samples as th 
though some difficulty 
In fact, th 
trometer is able to analyze actual fuel gases of t 


just studied, even 
experienced around mass 28. 
where 20 mixtures ar 


type component 


unusual. Such complicated mixtures hav: 
difficult when analyzed by a combination of 


methods. 


Jean Doyle and Marthada Vaughn k 


checked the calculations of heating valu: 


specific gravity, and this assistance is grat 
acknowledged. 
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Infrared Absorption Spectra of Hexafluoroethane and 
of Chloropentafluoroethane 


By Jose R. Barcelo ' 


In the present paper the infrared absorption spectra of hexafluoroethane and chloropenta- 


2 to 40 microns, and 


rch 


C. \l. Riehards, and H. L 


fluoroethane have been studied The spectra examined ranged from 2 


for this purpose lithium fluoride, sodium chloride 


potassium bromide, and thallium bromide 


iodide (KRSS5 prisms were ised The substances studied were in the gaseous state 


This work was undertaken to compare 


spectra of two molecules of the same tv pe 


of structure but with much different properties of symmetry rhe spectrum of hexafluoro- 
ethane has been previously measured in the infrared to 22 microns, and only small differences 
are noted in the present work from the previo study lables are included, which classify 
the observed spectra of hexafluoroethane d chloropentafluoroethane as fundamental, 


combinatio and overtone bands 


I. Introduction 


This paper presents the results obtained in the 
study of the infrared absorption spectra of two 
fluorocarbons, hexafluoroethane and chloropenta- 
fluoroethane, from 2 uw to 40 u 

In spite of the interest in fluorocarbons, in view 

their characteristics, no data were available on 

nfrared absorption spectrum of chloropenta- 
roethane. The infrared absorptance of hexa- 
oethane has been studied by J. R. Nielsen, 


MeMurry [1], 


Raman shifts for this material in the liquid state 


and the 


been measured by D. H. Rank and E. L, 
Pa Z In the present work, the data obtained 
xafluoroethane are in good agreement in the 
position of the bands with those given by Nielsen, 
1]. The study of the infrared absorption 
im of chloropentafluoroethane in compari- 
vith that of hexafluoroethane, from which it 
s only slightly chemically, is of interest be- 
as a result of the symmetry of the latter, 

| of its fundamental vibrations are degen- 
Because of the asymetry of the former, 

er, all of the eighteen fundamental fre- 


ies of chloropentafluoroethane are active in 


lnirored Spectra of C.F, and C.CIF 


the infrared. Unfortunately, of the eighteen fre- 
quencies expected, seven were not observed, eithe r 
because they were very low in intensity or existed 
in a region of long wavelengths not examined. <A 
majority of the observed bands may be accounted 
for as combinations of the fundamental frequen- 
cies Observed, although some difficulty in inter- 
pretation exists for three of those of lowest wave 


numbet 
II. Experimental Method 


The materials used were of the group of “‘freons,”’ 
supplied by duPont, and the spectra were obtained 
with these substances in the gaseous state. <A gas 
cell of 5 em in length with potassium bromide 
windows was used, except for the longest wave- 
lengths. The initial pressure of the gas in the 
cell was 500 mm of mercury, but as this was too 
high for the observation of some of the absorption 
maxima, it was reduced in each case to the most 
suitable one. The pressure used is indicated on 
the curves. When measuring the long wavelength 
region with a thallium bromide-iodide prism, a 
l-m cell with windows of polystyrene was used 
On account of the lack of strength of the poly- 
stvrene windows, it was not possible to evacuate 
the cell 


to flow in slowly for several minutes 


The cell was filled by allowing the gas 


When the 
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bands were too intense for observing their struc- 
ture, some of the gas was removed by flowing air 
through the cell 


vas in the cell was reduced to a level that vave a 


slowly 


good contrast in the different parts of the band 
By this 


method of reduc ing the quantity of the vas, it Was 


no further removal of the gas was made 


hot possible to determine the partial pressure 
The spectra of these two substances between 2 bu 
and 40 « were recorded with Perkin-Elmer spec- 


12A and 


( hlor ide 


trometers, models 12B, using lithium 
bromide 


With both 


spectrometers scanning and slit drive were auto- 


fluoride, sodium potassium 


and thallium bromide-iodide prisms 
mati Separate motor drives were used, one for 


the wavelength drum and the other to open the 


Fiaut 


slit continuously to achieve approximately con 
stant background energy at different wavelengths 
found 


thr se devices can be 


Pivlet 


Further details of 


in a paper by E. hk R. Stair, and C. J 
Humphreys ’ 

The detector used for the regions observed with 
the lithium fluoride sodium chloride, and potas 
sium bromide prisms was a vacuum thermocouple 
with potassium bromide windows In order to 
detect the weak radiation in the thallium bromide- 
iodide region, a Golay detector was used with a 
window and a lens of thallium bromide-iodide to 
focus the energy on the receive! In both cases 
the signal from the receiver was amplified and 
recorded automatically 

For the recording of the absorption spectra in 
the region bevond liu the tee hnique deseribed by 


Plvler [4. 5] was followed 
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When the amount of 


III. Experimenta! Results 


In figures 1 and 2 are shown the plots of t! 


centage transmission of hexafluoroethan 
chloropentafluoroethane as function of the 
length. Figure 1 shows the short wav: 
region for the two substances, and figurs 
long wavelength region. The dotted parts 
spond to regions in which absorptions of 
vapor and carbon dioxide appear in the 
ground energy, and the spectra in these rm 
have been established with the help of the r 
ing spectrophotometer of Baird Associates 
cancels out the atmospheric bands 

In table 1 are shown the results for hexafl 


ethane, first in brackets the intensities, thi 


Irequencies observed and in the next columt 
P.O, and R branches when 


re presentee 


values for the 
served The intensity scale 
proportional to log ] ] 


show no transmission with a 


numbers that are 
bands that 


quantity of gas are assigned the intensity 10 


the intensity number of the weaker bands 


relative to this band except that (0,0 
bands that have between 2- and 5-percent 


tion and (0.0.0) to bands of 2 pe ent or 


absorption 
The experimental results obtained by obs« 
the infrared absorptior Sper of hexafluoroet 


in the present work show slight variations 


those reported by Nielsen, et al. |] 
The intensities of the absorptlolr bands di 
ol the 


on the cell length and the pressure 


Except for small changes in the transmissiot1 
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WAVELENGTH IN MICRONS 


oroethane and oat 


small bands are listed in the tables of Nielsen 
et al. with values between 750 and 1,000 em 

\{ scanning of the region from 25 to 40 u did not 
reveal any bands of sufficient intensity for detec- 
tion The absence in the 25- to 40-u region of 
any bands of appreciable intensity supports the 
assignment of Nielsen that v, is less than 250 
em ', the value 216 em =! being used 


The band at 518 ecm™' appears to have P, Q 
and R branches when it is observed at a pressur: 
of 100 mm He (see fie. 2 With larger amounts 


ot gas the structure of the band is not observabl 


and the region of greatest absorption is at 522 


cm The experimental results for this region 
nre not shown by Nielsen and it is not possib] 


to determine the part of the band that he used i 


determining the value of 522 em 


The combination band 


\ re ie alor vhen the number of Nielsen, et al. but does not appear in 


s in the path of the radiant energy is the the experimental results. This band 


he intensities of the bands are proportional intensity and, as observed in this work, has a 


} of molecules in the path In transmission of 40 percent at its Maximum whet 


otal number 
ol Nielsen eta cvreatel absorption wis th Tas Is measured with ‘ re | f 500 mn 


in the bands because they used a larger Hg in the 5-em cell 


vy of gas in the absorption path for their In veneral the two set 


ments. This may, in part, account for closelv, and in some bands the 


Terences in experimental results in the et al. show better detail. These slig 


om 10 to 13.5 uw between these two sets in ‘veral parts of the spectrum 


In the pres work only one band Is observations being carried ou 


d in this reg at S94 em while several amounts of hexafluoroethane in th 
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In table 2 the values of the absorption maxima 
found for chloropentafluoroethane are given; the 
order of the wave numbers of the bands is the 


same as in table | As this substance had not 


iW ‘ ambers anc ( ¢ milensilies 
f f silt 


IV. Discussion 
1. Hexafluoroethane 


The substances of this type admit two possible 
configurations, the alternate one with the sym- 
metry model Dsg and the eclipsed one with the 


model DD, Nielsen, et al 1], after a detailed 


discussion, preferred the alternate form, which is 
Hamil- 


Cleveland [6] for the case of hexa- 


in accordance with the hypothesis of D. T 
ton and F. F 
fluoroethane. Their analysis was made from the 
study of the Raman spectra 

(Adopting this hypothesis, and taking into con- 


sideration the Raman data given by Rank and 


Pace [2], the fundamental vibrations will be those 


that are indicated in table 3, in which are shown 
the term, the wave number (cm the activity 
and the group from which they originate 

The values given, as far as the frequencies active 
in the infrared are concerned, differ only from those 


by Nielsen, et al 


frequencies vr and vy, On the othe 


1] by a very few units for the 
hand the 
frequency must be less than 250 em In 
Nielsen’s work it is supposed to be 216 em™', and 
with this value concordant, numbers for the com- 
bination bands can he obtained The region olf 
low frequencies was examined and no absorption 


maximum was found, thus giving the impression 
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been previously measured, a comparison wit 
experimental results is not possible Neil 
table 1 nor in table 2 are given those 
branches not completely resolved. 


ed bands of oroethane 


> 


raBLe 3 Fundamental bands 


listing of hand de signation, activity 


of wmbration 


that any possible absorption band must bh: 
region of an even lower frequency in accor 
with the assumption by Nielsen, et al 
Considering the data obtained here, the p: 
interpretation of the nonfundamental frequé 
as combination tones has been made followi 
1], and the differ 


between the observed and calculated valu 


indications of Nielsen, et al 


of the same order as those of the authors 
tioned, except In two cases where these 
were a little larger. The combination ton: 
sponding to »,+ »;, forecast but not fou 
these authors, is observed with the lithium fl 
prism at 2,534 em These observed and 


lated frequencies are to be found in table 4 
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Observed and calculated wave numbers 


combination bands of herafluoroethane 


2. Chloropentafluoroethane 


presence of the chlorine atom destroys the 
the 


modes of vibration, all active in the infra- 


try, which permits eighteen funda- 


Unfortunately, some frequencies 


) appeal 
edicted to be so low that they are not in the 
which extends to 250 


studied, cm 


which were supposed to be in the region 


est frequencies studied, have not appeared 


v0 seen in the curve. “The region from 25 


u shows only one absorption maximum 


five were supposed to be 
h the 
Irequencies, 
een established To 
hands, the 
they appear were taken into consideration 


sly 


the 18 funda- 
indicated in table 5 


the 


data obtained, 11 of 
whi h ure 
assign funda- 
intensities and the regions in 
them, bringing them in relation with 


ioroethane, “pri frequencies have been 


oethane 


ed Spectra of C.F, and 


introduced where the double degenerates of hexa- 
fluoroethane give way (to two different frequencies 
in the chloropentafluoroethane. In table 5 the 
if found, 


the group from which it originates, are indicated 


term, the value of the frequency and 


Although only a few of the fundamental fre- 
quencies are known, the remaining bands can be 
explained as harmonic or combination tones, with 
the exception of the three with lowest wave num- 


bers. For these, surely the lowest fundamental 
frequencies, which have not been found, must be 
taken into consideration. These 


are given in 


table 6, 
Observed calculated wave numbers of the 


TABLE 6 and 


combination and overtone bands of chloropentafluoroethane 


The 
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acknow ledged 
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Development of Very-High-Frequency Field-Intensity 
Standards’ 


By Frank M. Greene and Max Solow 


on of } development of two verv-high-frequeney field-intensit 
are being used at this Bureau for the calibration of commercial field 
he range 30 to 300 megacveles. These standards are emploved to establis! 


s of field intensity by either of two methods: (a) the standard-antenna method 
circuit voltage at the center of a receiving dipole is measured directl, b 
method in which the current at the center of a transmitting dipole 


The techniques used for determining the antenna current and voltag 


tribution on the antenna is determined theoretically, using Schelkunoff's 
ves the effective length These values are compared with those obtain 


i tests at 100 megacyeles are presented in hi the above two methods 


tercompared using horizontal polarization Their accuracy and limitation 


I. Introduction This necessitated the development of accurate 
field-intensity standards for calibration purposes 
number of years the National Bureau of Such standards are usually employed to determine 
rds has maintained among its other public a value of field intensity in terms of either a known 
s, the calibration of commercial field-inten- voltage or current in an antenna, together with 
ters. These instruments are largely used certain geometrical relationships involved, such as 
neering consultants and broadcast engi- the antenna length, distance of separation, et 
n order to determine a radio station's The standards developed are similar in principl 
efficiency and coverage area, acceptable to those already in use at the lower frequencies 
of which are specified bv the Federal However, it was necessary to modify the tech- 
nications Commission The frequency niques employed for the measurement of antenna 
f this calibration service has in the past current and voltage as well as the over-all tech- 
ed to approximately 30 Me and is main- nique in the use of the standards in order to adapt 
by the Bureau’s Central Radio Propagation them to use at the higher frequencies 
tory In many instances the accuracy with which the 
recently, as a result of the greatly increased value of a physical quantity may be determined 
the very-high-frequency (VHF) band (30 increases with the simplicity of the method used 
Me) for frequeney-modulation and _ tele- to make the determination. Such measurements 
roadcasts, * and other communication pur- usually involve a theoretical calculation, or an 
the Cent: idio Propagation Laboratory experimental observation, or both. In establishing 
ertaken extend the frequency range ol the field-intensity standards discussed herein the 


bration service to include these frequencies methods adopted for current and voltage measure- 


field Intensity Standards 





ment were made as elementary as 
It is believed that this 


enables the achievement of the maximum accuracy 


purposely 
possible for this reason 


possible or practicable at the present time. This 
of necessity results in some limitation of the range 
of measurement and flexibility of use, but the 
premium paid results in increased accuracy. 

The standards or methods used to establish 
accurately known values of VHF field-intensity 
are (a) the standard-antenna method and (b) the 
standard-field method 

It is the purpose of this paper to describe two 
such experimental standards in some detail, the 
techniques involved for measuring the antenna 
voltage or current, as well as the propagation 
tests made to intercompare the field-intensity 
values obtained by both methods. The accuracy 
and limitations of the two methods are also 
discussed 

This work done by the authors represents one 
phase of a program for the development of accurate 
VHF field-intensity calibration standards. Al- 
though these standards were developed to cover 
the entire VHF band, the initial propagation 
tests described herein were made at a frequency 
of 100 Me, using horizontal polarization only; 
and the standard antennas were horizontal half- 
wave self-resonant dipoles in all cases. 

No claims are made as to the originality of the 
individual methods or techniques used, as most 
of these have been reported elsewhere. However, 
the standards discussed here are novel at least 
in their simplicity and in the accuracy obtainable. 
Practical rationalized mks units are used through- 


out this pape! 


II. Standard-Antenna Method 


1. Relation Between Field-Intensity and Induced 
Voltage 


The magnitude of the electric component of 
VHF field-intensity existing at a given point in 
space may be determined in terms of the voltage, 
\’,., induced ina standard receiving dipole immersed 
in the field, together with the antenna geometry 
It will be 


antenna is a horizontal half-wave dipole oriented 


assumed in the following that the 


for maximum response, and that the voltage is 
referred to the center terminals In case the 
antenna is not open-circuited, the voltage meas- 
ured, V,, is of course not the induced voltage 
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but is related to it by the voltage-transf. 


where 7, is the load impedance, in ohms, cor 
to the antenna terminals, and Z, the a 
input impedance in ohms 

The magnitude of the electric compo! 
field-intensity (in volts per meter) is rela 
the induced voltage by 


where /, is the effective length of the ante: 
meters.' 


2. Effective Length of the Antenna 


Assuming a sinusoidal current distribution 
the antenna, the effective length (in meters 


r 
li tan 
T 


/ 


where \ is the wavelength in meters. and 
antenna half-length in meters 
For a half-wave dipole the effective le 


in meters) reduces to 


As is known, the current distribution on 
antenna is not exactly sinusoidal except 
infinitely thin filament. 
of radius a and half-length /, the departure f 
a sinusoidal distribution becomes progressi\ 
greater as the ratio 2//a decreases (i. e 
antenna gets fatter 

An approximate solution for the free-spac 
rent distribution on a cylindrical transn 
antenna ? has been obtained by Schelkunol! 
among others, and is given by the sum of 


and 35 in the appendix. The relative distri! 
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For cylindrical anten 





i Theoretical free-space relative current distribu- 
on a cylindrical half-wave dipole as determined from 
kunoff’s anolysis 

current J(r)/1(0) is shown versus the fractional half-length r// for 


' 


a lj. 


ies of the average characteristic impedance, Ko=120[loge(2 
in the appendix 


ied from the sum of eq 34 and 35 
inction of both the phase length 28l, and 
ratio 2//a, and may be represented by 


I(r 


10 f(B1K, 


current in amperes at a distance r 
meters from the center of the antenna 

current in amperes at the center of the 
antenna (7r=U 

»/ 

A 120 (log. ~ 


a 


| ) ohms, and is the aver- 


age characteristic impedance over 


the half-length / 


2r/X. 


relative current distributions as determined 
he magnitude of the sum of eq 34 and 35 
own in figure 1 for half-wave dipoles for 
values of Ky. Because of the increasing 
xity of the higher order Bessel functions 


infinite series of eq 35, only the first two 


of this series were evaluated. This ac- 
for the fact that the current as shown in 


does not reduce quite to zero as it should 


intensity Standards 


at the end of the antenna, r//=1.0. The contri- 
bution from the neglected terms in this series is 
appreciable only near the ends of the antenna, so 
that this approximation does not appreciably 
affect our results here. 

The effective length of the antenna (in meters) 
from the relative current 


may be determined 


distribution by the relation 
I(r) 


lis é 10 dr 


The magnitude of the sum of eq 34 and 35 may 
be substituted in eq 6 and the integration per- 
formed either mathematically or mechanically. 
We preferred the latter as being less involved. 

In this manner a number of values of effective 
length, /y, were determined for various values of 
the average characteristic impedance, Ao, between 
400 and 1,200 ohms. 
l, over that of an infinitely thin filament, Ay= @, 
is shown versus A, in figure 2 for half-wave 


The percentage increase in 


dipoles.* 

The agreement between calculated and observed 
values of /, for the antenna used will be discussed 
in section LV. 








a 
. HM 


, 


IGURE 2 Percentage increase im effective le ngth of a 


cylindrical half-wave dipole (over that of an infinttely thin 
filament) versus the average characteristic impedance, Ke, 


as determined by substituting the sum of eq 34, and 35 in 6 


The percentage incr 1 /‘w for the case of thin self-resonant dipple 


antenns nearly the same as that for full half-wave 


lipoles for the values isually encountered 
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In order to operate a half-wave dipole at self- 
resonance it Is necessary, as is known, to shorten 
its length somewhat from a full half-wavelength. 

Schelkunoff [2] and King [3] have obtained inde- 
pendent approximate solutions vielding the re- 
self-resonant operation of 
Although obtained by 


quired shortening for 
thin evlindrical antennas 
different 
required shortening as determined from Schel- 


radically methods, the value of the 
kunoff’s solution is in substantial agreement with 
that obtained from King’s second-order solution, 


and is given by 


}/ ; HON) (Qn 97 OS 
A rh k, ’ 


where Sita and Si(2z L.41S815. 


Values of the percentage shortening determined 
from eq 7 are shown versus the average character- 
istic impedance of the antenna Ao, in figure 3. 
3. Antenna Voltmeter 


The standard receiving antenna used was a self- 
resonant (half-wave) dipole shortened approxi- 
mately 4 percent from a full half-wavelength (/ 


100.0 Me 


self-resonant. 


as determined from fig. 3) to make it 
Dural tubing *{. in. in diameter 
was used, giving a ratio 2 //a~600, and an average 
characteristic impedance Ay=120 {log,(2//a)—1] 
~650 ohms. The induced voltage referred to the 
center terminals of the antenna was measured 
directly by means of a relatively high-impedance 
balanced voltmeter connected across the gap at 
the center. In this manner the necessity for a 
separate measurement of the antenna input im- 
pedance was eliminated, which greatly simplified 
the problem. The voltmeter consisted of a modi- 
fied type 1N28 silicon-crystal rectifier built into 
the gap as shown in figure 4. The crystal output 
was filtered by means of a balanced resistance- 
capacitance network, and the d-c output voltage 
was measured directly against a standard cell on 
a special precision slide-wire potentiometer. The 
experimental receiving 


central portion of an 


antenna assembly minus the antenna rods is 


shown in figure 5. The input impedance of the 
resistance-capacitance filter was approximately 
10,000 ohms in shunt High- 
frequency resistors were employed, which were 


with 0.05 wuf. 
available commercially. 
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entade shorten ng ofac ylind ca 
d for se 


in pedanc é, 


Ficure 3 Per 


dipole require -resonant operation 


A ~~ detern 


average characteristi« 


eg 
i 


The effective shunt capacity of the crys 


diode was less than 0.75 uuf. However, sinc 
resulting reactance of some 2,000 ohms was 
quadrature with the radiation resistance of 
antenna, the resulting shunting error from 
source was a small fraction of 1 percent. 

In selecting the crystal rectifier it was necessa 
resistance and 


to consider both its r-f input 








~ 
| 


Cc 


FIGURE 4 Vethod used to mount the ty pe IN 
rystal rectifier in the gap at the center of the 


standard receiving antenna. 


A, Western Electric 1N28 silicon cr 
modified 1N28 as mounted in the antenna g 
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iew of the central portion of an e2 perimental 
d receiving antenna minus the antenna rods showing 


lal rectifier and resistance-capacity filter netwo k. 


Both are a 
tion of the applied r-f-voltage level, decreasing 


ternal) d-c-output resistance. 
th nereasing level, although not necessarily at 
the same rate. 
In order to limit the shunting caused by the 
resence of the crystal voltmeter in determining 
open-circuit antenna voltage, it was necessary 
select a crystal having as high a value of r-f 
sistance as possible. On the other hand, to 
iintain as high a sensitivity as possible in 
easuring the d-c-output voltage of the crystal, 
| Was necessary to select a crystal with a relatively 
v value of d-c-output resistance. 
\ satisfactory compromise was obtained by 
ng a crystal having an r-f resistance between 
10 and 10,000 ohms when measured at an r-f 
of 0.1 v and self-biased. This resulted in a 
ting error of less than 2 percent in determin- 


induced antenna voltage. It was found in 


| that the type 1N28 crystals came closest 
eting these requirements. 

variation of both the r-f resistance and the 
{put resistance versus applied r-f voltage at 


YI \T 
i 


« is shown in figure 6 for the 1N28 crystal 
| when operated self-biased. These crystal 


teristics are not representative of the 1N28 
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crystals in general, as it was possible to find many 
with considerably higher or lower values of r-f 
and d-e resistance than the values shown. The 
variation among other crystals of types 1N21 to 
1N27 was found to be even greater. 

The r-f input resistance of the 1N28 crystal was 
measured by a _ resistance-variation method [4] 


using several high-frequency resistors having 
widely separated values as cross checks in deter- 
mining each value of crystal resistance. The 
crystal was connected in series with a 200-yuf d-c- 
blocking condenser across an antiresonant circuit 
during these measurements. The d-c-output volt- 
age acting as self-bias on the crystal was developed 
across this condenser. 

The d-c-output resistance of the crystal is the 
equivalent internal dynamic source resistance of 
the d-c-output circuit. This is sometimes referred 
to as the d-c, or video, impedance [5] and may be 
defined (in ohms) as 

Ry= ~, 
Ai 
where Ae balancing 


voltage of the slidewire potentiometer and Ai 


a small change in the d-c 
the resulting change in the d-c current in the 
output circuit. If Ae is kept smaller than about 
5 mv, when measuring #,, either a positive or a 
negative value of Ae will result in very nearly the 


same value of R,. 
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FicurRe 6 Variation of the r-f input resistance and dynamic 
d- output resistance as a function of the applied r-f 
oltage level for a selected ty pe 1N28 silicon crystal rectifie r. 
Che measurements were made at 100 Me with full self d-c bias applied to 


These curves are not representative of the type 1N2S crystals in 


d-c output resistance; @, r-f resistance 
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The d-e output resistance of the crystal was 
determined in this manner for a number of r-f 
voltage levels. The crystal was removed from the 
antenna and placed in a coaxial mount fed by a 
25-ohm r-f source. The d-c output voltage was 
measured across a 200- uuf blocking condenser in 
the ground-return circuit 

As the d-e output voltage of the crystal was 
measured on a precision slidewire potentiometer, 
it is necessary to define the sensitivity of this 
instrument. This determined the useful range 
of d-c output voltages that could be measured, and 
hence the resulting range of applied r-f voltages 
that could be determined. 

The voltage sensitivity of the potentiometer (in 


volts) is 


Ae Rai 


smallest change in the balancing voltage 
that can be set on the slidewire dial of 
the potentiometer. 
Ai=smallest unbalance current readable on 
the galvanometer 
R= total series d-c resistance of the measur- 
ing circuit and crystal. 

The potentiometer used in these experiments 
was a commercially available portable type, weigh- 
ing approximately 15 lb. It housed a standard 
cell, dry cells for comparison purposes, and a 
sensitive light-beam galvanometer, together with 
the slidewire dial and standard resisters all in a 
portable carrying case. The instrument had a 
Other 


instruments were available having a range up 


useful measuring range of 0.5 to 161 my. 


to 1.6 Vv. 

The most sensitive portable light-beam galva- 
nometer available was used for measuring the out- 
put voltage of the crystal. Its sensitivity was 
approximately 0.03 wa per scale division (1 mm). 
Assuming that one could detect a change of cur- 
rent, Ai, of 0.1 division, and that the total series 
d-c resistance of the crystal, filter, and potenti- 
ometer was approximately 20,000 ohms, the volt- 
age sensitivity was de= RAi=2 10**310~° 
60*10-° v. Thus, neglecting the other errors in 
the instrument, which 
voltages of the crystal could be determined to an 


were small, d-c output 
accuracy of better than 2 percent from something 
less than 3 mv to over 1 v (the safe upper voltage 


limit of the crystal This corresponded to a 
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range of r-f input voltage to the crystal of 
0.02 to 1.0 v as will be shown. 


The silicon crystal voltmeter was cal 
in terms of a standard voltmeter. On 
used at 100 Me comprised a quarter-way: 
balanced transmission line of accurately 
characteristic impedance terminated in 
viously calibrated VHF vacuum thermo: 

As is known, the magnitude of the trans! 
pedance of a quarter-wave length uniform Joy 
loss transmission line is equal to its characteris 
impedance [6], i. e., 

V, 


7, |= Zo» 


Zr 

where 

\’;=sending-end voltage 

/,=receiving-end current 

7o=characteristic impedance. 
Thus the sending-end voltage may be determin 
in terms of the receiving-end current and thy 
characteristic impedance of the transmission lin 
This relationship is independent of the magnitud 
or phase-angle of the terminating impedance (i 
this case the heater of the VHF vacuum therm 
couple). The line losses involved were found to ly 
small enough so as not to appreciably affect th 
results. 

The frequency error of the thermocouple was 
determined in terms of small bead-type thermistors 
used as current standards. These are known t 
maintain their d-c value 
frequencies well above the VHF band when biased 
The circuit used 


of resistance up 
to low values of resistance. 
to determine the frequency error is shown in figur 
7. The bridge circuit measured the d-c resistan 
of the series combination of the thermistors a! 
thermocouple heater biased with either «d- 

r-f current. A measure of the r-f current wa 
thus provided in terms of a known d-c curren! 
A negligibly small d-c measuring current was us 
in the presence of the r-f current so as not 


appreciably alter the thermistor resistance. T\ 


matched thermistors were used to preser\ 
symmetry, as the frequency correction 0 
thermocouple will in 
different if used in an unbalanced circuit 


general be consid 


The desired value of direct current wa 


established in the thermocouple heater. Kesist 


The VHF thermocouple referred to herein has a separate heat 


ym the couple by means of a small glass bead 
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(D- 166382 
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TC “balun” 
BALANCED 
rf INPUT) 


7 Circuit diagram of the apparatus used to deter- 
the frequency correction of VHF vacuum thermoele- 
when used in a balanced circu arrangement 


L™ anh; C™200uul(f = 100 Mc 


ance R, was then adjusted to balance the bridge, 
t the same time maintaining the value of current 
by readjusting the battery voltage if necessary. 
\fter measuring the d-c voltage output of the 
themocouple, the direct current was reduced in 
alue to 5 percent or less. The r-f voltage of 
lesired frequency was then applied, adjusted in 
value until the bridge was rebalanced, and the 
l-c voltage output of the thermocouple measured. 
The r-f voltage was then removed and the value 
of direct current determined to reproduce this 
The frequency 
correction was obtained from the two values of 
With sufficient sensitivity in the 
bridge and associated controls, it was found 


value of thermocouple output. 
direct current. 


possible to reproduce these frequency corrections 
toa small fraction of a percent. 

Generally speaking, for a given accuracy, a 
current standard is somewhat easier to provide 
at these frequencies than is a voltage standard. 
Fort inately, the relation given by eq 10 enables 

to provide an r-f voltage standard having an 
racy approximating that of the current 
tandard used. The difficulty, however, in this 
ase is in the restricted frequency range over W hich 
tandard may be used. Normally a contin- 
oltage calibration is not required over the 

ney range of the standard antenna. A 
r of individual calibrations spaced at fre- 
Vv intervals of approximately 25 or 50 percent 
the desired range is usually sufficient when 
sponse is broadband. 


Intensity Standards 


Probably one of the most accurate methods 
of determining the characteristic impedance of the 
balanced low-loss transmission line used in this 
voltage standard is from the relation 


where 


characteristic impedance in ohms, 

velocity of propagation in meters per 
second, and 

mutual capacity between conductors in 
farads per meter. 


The type of transmission line used was com- 
mercially available unshielded twin-lead having 
a nominal characteristic impedance of 75 ohms. 
The velocity of propagetion was determined to 
an accuracy of better than 0.25 percent at 100 Me 
by measuring the resonant frequency of a half- 
wavelength section of line short circuited at 
both ends. The length of the short circuit used 
was such as to introduce a negligible error in the 
measurement. The calculated effect of the line 
loss on the velocity of propagation was also found 
to be negligible. The balanced mutual capacity, 
C’, per unit length was determined on a special 
precision Schering bridge at a frequency of 1,000 
cycles to the same order of accuracy. There is 
no reason to suspect that the capacity per unit 
length should change with frequency, at least up to 
Thus 
the characteristic impedance was probably deter- 
The charac- 
teristic impedance of a number of samples 1 m 


frequencies of several hundred megacycles. 
mined to better than 0.5 percent. 


in length cut from the same piece of twin-lead 
varied less than +0.1 percent from the average 
value of the lot. This would seem to indicate 
the exceedingly small departure from uniformity 
existing within a given sample of the line. The 
average value of Z, for the lengths measured was 
83.0 ohms. The variations between samples of 
different production-runs of twin-lead was quite 
large, being as high as 10 to 15 percent, as might 
be expected. The aging of the twin-lead seems to 
be negligible if not exposed to the weather or 
undue stress. This was indicated by having a 
value of the velocity of propagation of a particular 


sample repeat to within 0.1 percent after several 


months. 
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Another type of balanced voltage standard 
used consisted simply of a sensitive thermocouple 
previously calibrated in a manner similar to that 
How- 


ever, in this case a single thermistor was used in 
1.000- 


already described and shown in figure 7. 


shunt with the thermocouple heater. <A 
uuf ceramic disk condenser was used in series with 
heater lead to provide de 


each thermocouple 


isolation and preserve the symmetry. The 
leads were kept as short as physically possible 
(’s In. In Most cases 

The voltage standards described above were 
used at only one voltage level (approximately 
1 v), that 


through the terminating thermocouple 


corresponding to the rated current 
The volt- 
age range was extended downward by means of a 
mutual inductance) attenuator 


better than O.ldb (ap- 


precision piston 

6a] having an accuracy 
on 

1.0% 


total rms harmonic content of the r-f generator 


proximately over a range of 60 db. The 
output was less than | percent, so that any result- 
ing error in the voltage standard from this cause 
was negligible. Pi matching networks were used 
in both the input and output circuits of the stand- 
ard attenuator to decrease the insertion loss. 


blo« k 
equipment and its arrangement for one of the 


Figure 8 is a diagram showing the 
methods used to calibrate the balanced crystal 
The antenna rods were removed dur- 
A broadband balun | 


to transform from the unbalanced output of the 


voltmeter. 


ing calibration 7] was used 


standard attenuator to the balanced circuit of 
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FIGURE 9 Typ cal calibration 
meter, showing the variation in 


appl ed r-f oltage for the 


the antenna and voltmeter. The “ Pi’ mate! 


networks are not shown Two 200-uul 
blocking condensers were used to prevent 
circuiting the d-e output of the antenna crysta 
Typical calibration curves of the antenna 
meter made at 100 Me are shown in figure 9 
usable range of r-f input voltage was from 0.0) 
1.0 v and the probably accuracy about +2 


The solid | 
and curve joining them represent the first o 


cent over most of the range 
calibrations. The circles represent a repeat! 
bration 2 days later at a temperature | 
lower. 

The data of figure 9 show the order of repr 
bility possible with such a crystal voltmeter 


partially controlled temperature conditions.  T! 
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er of the individual points represents ap- 
ately 5 percent of the voltage, so that any 
touching the curve is in agreement to within 
percent, 

vas usually possible to repeat a crystal cali- 
m to within 1 percent if made within a few 
at the same ambient temperature. In no 
vas the calibration relied on for more than a 


naking the field-intensity measurements de- 
din section LV, using this standard-antenna, 
ystal voltmeter was always calibrated in the 
within an heur of the time the tests were 
mpleted. A record of the temperature was kept 
nd usually did not vary more than 2 deg C during 
ny test. Thus all that was required of the crystal 
oltmeter was a good short-time stability. 
\lthough the effect of variations in temperature 
the accompanying hysteresis are known to 
ontribute largely to the instability in a crystal 
voltmeter, no data were accumulated on the exact 
wgnitude of these effects. We have so far avoided 
s by operating under as smal] a temperature 
ration as possible Thus it would probably 
not be feasible to calibrate the erystal in a labora- 
ory and use it for measurements outside at a 
temperature as much as 10 or 15 deg C different 
lhe response of the antenna crystal voltmeter 
as found to be essentially independent of fre- 
quency up to 200 Me with a rising response beyond 
output indication hieh due to series resonance 
vithin the erystal and mount. The series-resonant 
equency for this voltmeter was of the order of 
2,000 Me 


III. Standard-Field Method 
l. Theory 


\ predetermined value of field intensity may be 
stablished at a given point in space in terms of 
irrent distribution in a transmitting antenna, 
effect of the ground, and the geometry in- 
ed 
ln the development of the standards being 
wrted here, only horizontally polarized trans- 
on over plane homogeneous earth having 
values of relative dielectric constant, €,, and 
ictivity, o, was considered. The distance of 
ation used between the transmitting and 
ving antennas was from 2 to 20 wavelengths. 
TI conditions simplify the resulting formulas 
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somewhat and for the accuracy desired are quite 
justified. This range of distances is usually of 
little interest to the propagationist, at these fre- 
quencies, and many of the transmission vagaries 
of concern to him are ignored here. At these 
smaller separations a more exact knowledge of the 
actual angle of incidence and ground reflection 
coefficient is usually needed in order to produce 
the required degree of accuracy in predicting the 
actual value of field intensity existing. The 
geometry of the transmission system is shown in 
figure 10. 

The root-mean-square value of the electric 
component of field intensity produced by a hori- 
zontal transmitting dipole at distances greater 
than about 2\ over plane homogeneous earth ts 
[8, 9, 10] (in volts per meter 


60rl,1 Te77** 
J \ R, 


Direct 


wave 


in which 
l,4=effective length of the antenna in 
meters 
A\=wavelength in meters 
R,=direct-ray path length in meters 
R,=ground-reflected-ray path-length in 
meters 
pe ’*=complex plane-wave _ reflee- 
tion coefficient (horizontal polar- 
ization 
complex surface-wave attenuation 
factor 
rms current in amperes at the center 
of the transmitting antenna 
k=2n/d 
j=yv-! 


The first term of eq 12 represents the field 


intensity that would exist if the transmitting and 


receiving antennas were located in free space 
The remaining terms take into account the pres- 
ence of the earth. The second term represents the 
vround-reflected wave which, when added vee- 
torially to the first term of eq 12, comprises the 
space wave. The third term of eq 12 represents 
the surface wave associated with the r-f currents 
actually flowing in the ground. 

The total field given by eq 12 is usually referred 
to as the ground wave. Propagation via the 


2 
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TRANSMIT TING 
ANTENNA 


[= 


RECEIVING 


ANTE NNA 


marking an 


ionosphere or troposphere, which is subject to 


changing solar or meteorological conditions is 
being ignored here 

The magnitude of the surface wave is usually 
negligibly small for horizontally polarized trans- 
mission at these frequencies [11], and in the re- 
maining formulas used herein we shall neglect its 
presence entirely. 

The complex reflection coefficient of the ground, 
I, may be expressed in terms of the angle y 
tan ~'(A,+A,)/d and a complex dielectric constant 


€), a8 follows, for horizontal polarization [12], 


sin ~W— v€& 
sin ¥+ y« 


cos" v 
cos" ¥ 


where 


€ «(1 j —) €,—]00Xe 


e-—relative dielectric constant of the ground 
referred to free space as unity 
e,€,, Where e, is the permittivity of evac- 
uated free space 
ob 
ag, 10 farad per meter 
rround conductivity, in mhos per meter 
Jai 
wavelength, in meters 


frequency, in cycles per second 


Values of the magnitude of eq 13, p, are shown 
versus tan ~=(h, +h, 
values of the relative dielectric constant, e,, for 


Many 


d in figure 11 for various 


low-loss dielectrics (¢/ew<! types of 
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ground may be treated as low-loss dielectri 
a large portion of the VHF band as far a 
reflecting properties are concerned. This 
ticularly true at the higher frequencies ab 
or 75 Me. 

In this case the phase shift on reflection 
very nearly 180° for horizon:»! polarizati 
that the value of the reflection coefficient is 
approximately by 

'a—~p 


Ignoring the surface wave altogether, the 1 
tude of eq 12 may be written (in volts per 1 
60r/ _1 ] l'¢ R,-R 
~ 


+ 


E » R R. 


Equation 15 may be expanded without f 
approximation, giving (in volts per meter 


6021, 1 
A 


( l 1 p ) 46 ata I R, ~o@ 1/8 
An t/” Re” 2 i 


where p is the magnitude of the reflection 
efficient, T, 
Equation 16 is in a form somewhat more suita 


Ek 


and @ is its phase lag in radiar 


for computation, since no loss of accurac\ 
involved as a result of the taking of small diffi 
ences between large computed quantities. 

In case the approximation given by eq 14 app! 
eq 16 may be written (in volts per meter) 


a 60x/,,1 | l ) 
a ( R, 


4p 
RR, 


The quantity ((1/2,)—(p/R,)? within the 
becomes negligibly small when d>(h,+h/, 


R, >R, >, po >1, and R, R, ~2hhe d. [ 


these conditions, eq 17 becomes (in volts per n 


Ei~ 1202/1 ”” ( 2rhj hz ) 


hd hd 


Equation 18 is obviously not valid 

vicinity of the zeros of the sine term, sin 
first term of eq 17 is not then negligible. — ‘T 
fore use of eq 18 should preferably be restrict 
values of A,ho/Ad<—0.5. If either the wavel 
\, or the distance, d, is increased, eq 18 will 
its last maximum value when h,h./Ad = 0.3229 
after which it will steadily decrease. If the an 
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< 














h, or hy above are varied, this maximum since sin @6~6@ to within the stated 


taceuracy for 
, a 
i when Ayho/d 0.25 6<1/4 radian 
h he Ad] 


+ the field intensity is given 


When attempting to establish values of standard 
ts per meter) to within 1 percent by field Intensity, however, one should rely on eq 16 
BE 24027] ,7 (hh 
~ ( ), 


or 17 until the assumptions applying to eq 18 and 
ad \ 


19 are found to be fully justified 
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2. Description of the Standard Transmitting 
Antenna 


The transmitting antenna used was similar in 
many respects to the standard receiving antenna 
described in section I] \ horizontal half-wave 
dipole was used shortened 4 percent for self- 


100.0 Nik 


diameter was used 


resonant operation at Dural tubing 


In. Im viving a ratio 
2] a~600 and an average characteristic impedance 


K 120 flog, (2la 1}~650 ohms. 


The current at the center of the antenna was 


measured by means of a previously calibrated 
VHF vacuum thermocouple and a section of bal- 
anced twin-lead transmission line one-half wave- 
length long. The antenna and _ thermocouple 
terminated the transmission line at opposite ends, 
and a balanced feed-line from the r-f generator 
was tapped in at the midpoint. Thus the cur- 
rents flowing at the opposite ends of the half-wave 
line were identical in value. The magnitude of 
the antenna current was thus indicated by the 
thermocouple \/2 distant and was entirely inde- 
pendent of the magnitude and phase angle of the 
terminations, namely, the antenna input imped- 
ance, and the impedance of the thermocouple 


heater A 
transform the unbalanced output of the r-f gen- 


broadband balun [7] was used to 
erator to the balanced transmission system. A 
block diagram of the equipment and its arrange- 
ment are shown in figure 12. The thermocouple 
was calibrated in a similar manner to that de- 
scribed earlier for the receiving antenna volt- 
meter, and was found to have a frequency error 
of 0.5 percent at 100 Me 


IV. Propagation Tests 


1. General 


In order to determine the probable accuracy of 
the VHF field-intensity determinations made by 
the standard-antenna and standard-field methods 
actual tests 


described, propagation 


100.0 Me, using 


previously 
were made at a frequen V of 
horizontal polarization only. The purpose was 
lo ac tually intercompare values of field-intensity 
determined with the two standards 

For instance, a predetermined value of field in- 
tensity was established by means of the standard- 
field method, and simultaneously measured using 
the standard (receiving)-antenna method, and the 
instead of relying on 


results compared However 
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Figure 12 


Diagran shou ng the method 
mining the r-f current at the center of the 


dipole 


isolated pairs of such measurements to determir 
the probable accuracy of the two methods, mo 
or less complete propagation runs were mad 
volving some 15 or 20 pairs of determinations 
each, in order to determine the trend with varying 
antenna height or distance of separation. In t) 
manner it was possible to determine the range of 
antenna heights and distances of separation | 
tween the antennas over which the previous as- 
sumptions applying to the earth were valid, ai 
to what extent we were dealing with a unifor 
field. 

Three general types of tests were made as 
lows: 

a Vertical-incidence standing-wavre meas 
ments, in which the standard field beneath 
horizontal transmitting antenna was probed 
tically using the standard (receiving) antenna 

b) Variable-height runs, in which the heights: 
the transmitting and receiving antennas wi 
varied separately for a fixed horizontal dis 
of separation. 

ce) Variable-distance runs, in which the 
zontal distance between the transmitting 
receiving antennas was varied for several 
antenna heights 

In order to conduct these tests it was nec 
to select a special site. The usual require! 
specify that the location be relatively flat an: 


of such reflecting objects as trees, buildings 
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rder to approximate the conditions as- 
to exist in deriving the standard field 
sin section II] 
of approximately 100 acres was made 
at the Beltsville, Md. Airport for these 
Che land had been ( leared and eraded and 
above requirements. Interference from 


was negligible, inasmuch as the airport 
licially closed and used only for emergency 
other requirement, dictated by the lack of 
ty in the standard antenna, made it neces- 
conduct such tests at a location where the 
tensity of frequency-modulation and tele- 
broadcast stations was negligible compared 
standard field being used. Such interfer- 
as further minimized (since the test site 
is sufficiently free of reflecting objects) by prop- 
wienting the twoantennas. This was possible 
case since the interfering stations, some 20 

es distant, were all approximately in the same 


ction from the test site. 


Vertical-Incidence Standing-Wave Measure- 
ments 


This run was originally intended to be used for 
termining the reflection coefficient and hence 
he apparent dielectric constant, and conductivity 
if the ground at the test site. However, it was 
lso found to be useful later in determining the 
fective length of the transmitting and receiving 
tennas, as will be shown. 

For this test the horizontal transmitting dipole 
vas located at the top of an all-wood ladder-mast 


ta fixed height h,=9.227 m. The balanced un- 


shielded 75-ohm r-f transmission line ran at right 
ingles to the antenna and approximately 10 deg 
low the horizontal. The unshielded balanced 
output lead from the thermocouple paralleled 
f line and was spaced approximately 3 ft 
negligible r-f pickup. The r-f power gen- 
and associated equipment were located on 
round approximately 200 ft distant. The 
mitting-antenna current was continuously 
tored by an operator at the generator and 
ld constant to within +0.1 percent. The 
tput of the thermocouple was measured on 
table precision potentiometer. The prob- 
curacy of the antenna-current measurement 
out 1 percent. The frequency of the r-f 


tor was also continuously monitored by 
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means of a calibrated heterodyne frequency meter, 
and held to 100.0 Me 
The standard 


t+ 0.1 percent. 


receiving) antenna (oriented 


parallel to the transmitting antenna), was mounted 
on a movable carriage on the same mast, and was 
readily adjustable in height by means of a contin- 
uous rope passing through nonmetallic pulleys at 
the base and top of the mast. Certain predeter- 
mined heights were marked directly on the side 
of the mast and read by a second operator located 
at some distance off the end of the antenna in a 
null. The actual height of the antenna above 
ground could be thus measured to better than 0.5 
em. The balanced d-c output of the receiving 
antenna crystal-rectifier was carried over rubber- 
covered lamp cord, which passed verti ally down 
the mast to the ground. From this point the d-c 
line went directly to the operator’s position where 
the voltage was measured on the precision slide- 
Tests 


were made for r-f pickup on the d-e line, and in 


wire potentiometer previously mentioned. 


the resistance-capacitance filter at the antenna, 
with the antenna rods removed, but was found to 
be negligible even with the d-c line running the 
entire height of the mast. 

Tests were made for stray radiation from the 
transmittingsystem caused by any slight unbalance 
in the feed line to the antenna or by leakage. The 
transmitting antenna rods were removed and a 

substituted to 
With the normal 
r-f current of 0.1 amp flowing, the resulting field 


65-ohm high-frequency resistor 


terminate the transmission line. 


was too small to be measured for any height of 
the receiving antenna. The vertical component 
of field existing at the receiving antenna was more 
than 40 db below the main field, with the hori- 
zontal transmitting antenna excited. 

Returning to the description of the vertical- 
incidence measurements, the resulting interference 
between the direct and ground-reflected waves, 
set up a standing-wave in the space beneath the 
transmitting antenna, as shown in figure 13. The 
standard (receiving) antenna was moved vertically 
up or down the mast to measure the resulting 
value of field intensity The reflection coefficient 
was accurately determined from the standing- 


wave ratio, as will be shown. Once this was 
known, the expressions for the field intensity for 
the standard-antenna and standard-field methods, 
eq 2 and 15, were equated and solved for the re- 


maining unknown, the effective length of the 
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antennas. This was deemed permissible since 
both antennas were physically identical. 

The field intensity as determined by the stand- 
ard-field method was obtained by substituting 


R, hy, ho, PR, hy, ; h . d 0 in eq 15. 


volts per meter 


viving (in 


. j2kh 

6021, 1 (20 
a h 2 T ho 

If h, >h.. 

eq 20 will have a maximum value if (2kh,+-¢) 

value if (2kh.+¢@ 


the value of field intensity given by 
2nr, 
minimum (2n—1) x, 
= es The magnitude of the 


and a 


Ww here n 


ratio of the maximum to minimum values of field 


intensity is then 

p 

SWR Pema m 

- p 
ho min 
The magnitude, p, of the complex reflection 
coefficient, T, of the ground may be determined 
from the SWR and the heights, h., of the receiving 
antenna corresponding to the maximum and the 
minimum values of field intensity measured plus 

the height of the transmitting antenna, thus 


SWR 
ean ae 
SWR 

hi+h 


The angle of phase lag, o, of the reflection co- 
efficient may be determined from the shift in the 
positions of the maxima or minima from the posi- 
tions corresponding to @=0. For the maxima 


this gives 


where ? | is the number of the maxi- 
mum used to calculate @ counting from the ground 
upward 


In using the minima to determine ¢, the expres- 


SiO) is 


24 


/ and fh. . are the act 


nihimum or 


ial heights, in meters, 
at which the particular maximum 
occurs 

In determining the value of the reflection co- 
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ANTENNA HEIGHT fh 
Ficgure 13 Variation of the observed and ca 
of electric field intensily, E in volts per m eler 


height, ho, in meters, of the standard receiving ant 


snd parallel t 


100.0 Me, J 


smitting antenna was located above 

s fixed height, 4: =9.27 m, d=0 
efficient of the ground at the test site, the hi 
h,, of the various maxima and minima wer: 
fully measured. It was found possible to 
mine these to within +1.0 em. To withu 
accuracy the positions corresponded fo a phas 
The pos 


error in determining @ thus amounts to appr 


shift, @, on reflection of 180 deg 


mately 2 deg. It is interesting to note tha 
calculated phase shift, ¢, for average groun 
15. o@ 5x10 
179 deg at 


of 2 deg in determining ¢ usually has a neg 


mho per meter) at 100 Me is 
nearly normal incidence \ 
effect on the magnitude p of the reflect 
efficient, and on the resulting value of fi 
tensitv calculated 

Although the standing-wave method cat 
used to determine even approximately th 
conductivity at these frequencies, such 
tion is usually not required if the refle« 
eficient can be accurately determined 

An example will be given of the met 
determining the magnitude, p, of th 
coefficient, and the dieleetrie constant. « 
cround from the observed data presented 
3. Three different values of standing-wa) 
were determined from the three maxima a 
three succeeding minima. The minimum 
ring at the ground (/ 0) was ignored b 
of its questionable value 
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value of p corresponding to each SWR was 


2, and the resulting 


etermined from eq 22, 


if the relative dielectric constant ¢, deter- 
from the relation 
p ) 


p 


from eq 13 by substituting 


p and solving for e,, 


was obtained 
and T (c/ew< ] 


sults are given in table 1. 


s determined from standing 


irements 


0. 648 
635 


636 


640 


The values of p and €, measured are determined 
gely by the moisture content of the ground. 
Consequently, their values fluctuate somewhat 
from day to day with changing moisture content. 
Four values of e, obtained in a similar manner over 
i period of a few months are shown in figure 14, 
together with values of daily rainfall. It can be 
seen that the variation in the dielectric constant 
is quite large. Normally, this is not of much con- 
to the propagationist, who is usually inter- 
sted in grazing incidence, but it is apt to affect the 
wcuracy of standard-field measurements at very 
Near 


the surface the moisture content and consequently 


short distances if the incident angle is large. 








witli 


the ground constants themselves may vary with 
depth. Since the skin depth in turn varies with 
frequency and angle of incidence (among other 
the determined for 


these constants may also depend on these param- 


factors (apparent) values 
eters. 

If desired, the reflection coefficient of the ground 
for vertical incidence may be calculated from eq 


I 


13 upon substituting y=7/2, giving 


Values of the magnitude of eq 26 are plotted 
15 


remembering 


dielectrics, 
for this 


for low-loss 


that 


versus ¢, in figure 


a/ew<< | T~ p 
condition, 

Determining the value of the reflection coeffi- 
cient, I, of the ground from the observed data of 
figure 13 actually required only relative values of 
With the value of 
now accurately determined, the only remaining 
20 is the effective length, l,, of 
the 
antennas 


the measured field intensity. 


unknown in eq 
the 
transmitting 
geometrically identical, the assumption is made 
There has 


transmitting antenna. Since standard 


and receiving used are 
that their effective lengths are equal. 
been some indication in the works of King [13] 
that such is not excatly the case, and further that 
the value of ly is also somewhat dependent upon 
the antenna termination. However, it is shown 
that 


are 


for self-resonant antennas these differences 
Shelkunoff 


but 


small does not elaborate on 


these points, our experimental values of 


~ 
c 
DIELECTRIC ¢ 


‘ 


x 
a 








| oo oi 


° 
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FIGURE Chart showing four 


@. Dielectic constant 1 daily rainfall 
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88078950 , 


paiues of 


, of the ground at one particular site, 


25 5S 0 1 20 25 
OCTOBER 
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SEPTEMBER 


dielectric 
100.0 Me. 


the relative 


as affected hy rainfall f 


apparent 





effective length agree well with those calculated 
from his analysis of cylindrical antennas, as will 
be shown. The value of the effective length, ly, 
may be determined by equating eq 2 and 20, and 
solving for ly 

The data shown in figure 13 provide 19 separate 
determinations of /y (ignoring the observed value 
corresponding to h,=0 because of its question- 
The average of these gives a meas- 
0.960 m for the self-resonant 


able value). 
ured value of ly 
antennas used. 

Figure 16 shows the difference 


between each pair of determinations as presented 


percentage 


in figure 13, versus the height, A., of the receiving 
antenna. In this information is included experi- 
mental error, differences caused by changes in the 
receiving antenna current distribution due to the 
presence of the ground, and differences caused by 
the proximity of the transmitting and receiving 
antennas, of which the latter should be small in 
view of the fact that the receiving antenna was 
operated open-circuited. These differences, in 
general, are less than 3 percent, and give a good 
indication of the agreement possible between the 
standard-antenna and standard-field methods at 
least for the case of normal incidence where the 
reflection coefficient of the ground could be 
accurately determined. 

It is also indicated that any varying effect (with 
antenna height) of the presence of the ground on 


NITUDE OF REFLECTION COEFFICIENT 


MAG 








> 
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IGURE 15 Calculated variation of the magnitude, p, of the 


com plez plane-wa e reflection coefficient of the gro ind A 


versus the relative dielectr onstant, é,, for vertical nei 


w/2. 


an 


JIM ELY 
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ANTENNA HEIGHT he , METER 


FicureE 16. Variation of the percentage difference 
the observed and calculated values of field inter 


figure 13, with height, ho, of the receiving antenna, 


the receiving antenna current distribution and 
hence upon its effective length is negligible, a 
least for the present purpose for antenna heights in 
excess of 0.1 wavelength. 

The calculated value of the effective length is 
obtained from the product of the following values 
(a) the classical value of effective length for a 
half-wave dipole given by eq 4; (b) the percentag: 
shortening for self-resonant operation from figur 
3; and (c) the percentage increase for cylindrica 
antennas from figure 2. The resulting value is 
ly d/r 0.96 X 1.0620.97 m for operation at 10 
Me, and agrees to within about 1 percent with th 


above measured value. 
3. Variable Height Runs 


The purpose of the tests described in this sectior 
was to determine the probable measurement error 
existing in the two field-intensity standards under 
discussion. Determinations made by the tw 
standards were directly intercompared, the maxi- 
mum probable error (in the absence of other 


being evaluated as one-hal! 


reference standards) 
the percentage difference between the two me: 
surements. Parallel self-resonal 
dipoles were used at a frequency of 100.0 Me 
Values of field intensity established by means o 


| 


horizontal 


the standard-field method were measured using 


standard (receiving)-antenna method and 
results compared. 

These runs consisted of varying separately 
heights of either the receiving or transm 
antennas for a fixed height of the other ant: 


and for a fixed distance of separation 
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17 View of the various pieces of equipment used in 


7 the measured data of fiqures 13, 14, 18, 19, and 


In the first two runs the receiving antenna was 
ed in height, Ao, from 0.265 to 9.0 m, using the 
mast assembly shown in figure 17. The 
chts, A,, of the transmitting antenna were 1.58 
.05 m, respectively, and the horizontal dis- 
of separation, d, was 30.5 m. The measured 
e of the effective length of the antennas deter- 
ned in the standing-wave test was used, namely, 
~U.90 Mm . 
lhe results are shown in figure 18. In deter- 
ng the field by the standard-antenna method, 
is used (in volts per meter) 


bss 
is the measured value of open-circuit 
value of field determined by the standard- 
ethod was obtained (in volts per meter) 
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<9— 50——s 


a)  WOml Lf l »\ 
e IN Mr rR) 


1 } 
R R, 


0.96 m 

0.100 amp 

3.00 m 
k=2x/d 
R [ (hy h d*| " 
R.=[(hi+hy)? +d] 


Values of the magnitude, p, of the reflection 
coefficient were determined by taking the mag- 
nitude of eq 13, assuming low-loss dielectrics, 
(c/ew< Values of p may also be taken from 


figure 11 


It was not found practical to measure the di- 
electric constant of the ground before each run 
because of the special setup required and the time 
consumed. Previous measurements showed that 
on the average the dielectric constant had a value 
around 15, except after a rain or after a prolonged 


dry spell. This value of e, was therefore used in 


calculating p from eq 13 








Ficure 18 Variation of the magnitude of the electric com- 
ponent of fu ld intensity, |/E n volts per meter as deter- 


standard-field 


methods vs the he ght, h n meters. of the standard receiv- 


mined by the standard-antenna and 


ing antenna 


intenna 
Horizontal 
nant half-w f 0.960 m 


lard-antenna I 


543 





The diameter of the circles representing the 


standard-antenna data in figure 18 is approxi- 


mately 5 percent of the value of field intensity 
Thus the difference is 2.5 percent or less for those 
Over most 


points that touch the solid curves 


of the range of heights, A., the agreement is within 
about 5 percent 

A similar run is shown in figure 19, except that 
the height of the transmitting antenna was varied 
The 
was maintained constant at 
The other 


as those applying to figure 18 


for a fixed height of the receiving antenna 
antenna current, /, 
0.100 amp conditions are the same 
The agreement 
between the two standards is about the same as 


in the two previous runs 
4. Variable-Distance Runs 


In these last tests the field intensity versus 


the horizontal distance of separation, d, was 
determined for fixed heights of the transmitting 
and receiving antennas 


The 


fixed 


maintained at a 
height Ae=3.05 m, and the 


receiving antenna Was 
distance, d, 
varied from 6 to 60 m for each of two heights of 








ov 
meter as deter- 


Standard-Field 


f the transm tling 








FiGurRE 20 


ponent of ela 


the transmitting antenna h/ 1.58 and 3. 
respectively 

Horizontal self-resonant half-wave dipoles 
The field 
mined for various distances of separation 
both the 


methods 


used as before intensity Was 


standard-antenna and = standar 
The remaining conditions are thi 
as those in the previous section. The resul 
shown in figure 20 

The values of field intensity determined by 
the standard-antenna and standard-field mi 
figs. 13, 18, 19, and 20) seem to be in fairly 


agreement over most of the range of a 
heights or distance of separation involved 
reasonable to assume that the effective len 


both antennas remains fixed throughout the 
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The larger disagreements noted are there- 
obably a result of errors in calculating the 
rd field because of the accuracy with which 
tual ground constants are known for the 
; incident angles involved. 


V. Conclusions 


\ description has been given of two standard 
hods of establishing known values of VHF 
ntensity used at the National Bureau of 
andards for the calibration of commercial field- 
nsity meters. Methods and techniques used 
alibrate and verify the accuracy of the stand- 
rds were described. Values of standard field 
nsity can be established in the range from 20 
m to Ilv/m at 100 Me 


increased by 


This range can be 


msiderably varying the trans- 


ting-antenna current 


eans of the 
The propagation tests described in- 


in known amounts by 


standard attenuator previously 
entioned., 
cate that the probable accuracy of the standards 
s about 5 percent under most of the conditions 
ised 
The limitations in the use of the standards lie 
n the somewhat restricted range of field intensity 
vailable, the temperature effect of the crystal 
tifier, and for some uses the lack of selectivity 
the standard receiving antenna. However, for 
i purposes reported here and under the condi- 
ons used these factors presented little difficulty. 
The initial work described in this report was 
nfined to the use of horizontal polarization at a 
frequency of 100 Me 


rify the accuracy of these calibration stand- 


Further tests are in progress 


ds at frequencies up to 500 Me under conditions 
mally used in the field. The additional use of 
eal as well as circular polarization will be 
luded in this program 
Under most conditions normally employed for 
ating VHF field-intensity meters, the value 
antenna constant ® as determined by means 
rizontally polarized transmission would prob- 
be substantially the same if vertical polariza- 
were employed This would be approxi- 
v true provided the antenna were immersed 
iniform field in both instances, and that its 


height above ground was greater than about a 
half wavelength. The antenna constant can be 
shown to be influenced by the effect of the pres- 
ence of the ground on the current distribution 
along the antenna, as well as on its input imped- 
Although, these effects are different for 
horizontal polarization than for vertical, they are 


ance. 


probably relatively small for antenna heights 
above approximately a half wavelength 

Results presented herein indicate that for the 
case of horizontal polarization any possible effect 
of the ground on the antenna current distribution 
probably results in less than about a 2-percent 
change in the antenna constant for antenna 
heights greater than 0.1 wavelength. 

The effect of the presence of the ground on the 
antenna constant as caused by changes in the 
antenna input impedance has been shown to be 
less than about 5 percent for antenna heights 
than 0.65 wavelength above 
ground (€,=15, e=5>107° mho/m) [14]. 


In the calibration of commercial VHF field- 


greater average 


intensity meters the accuracy of calibration as 
certified by this Bureau at the present time would 
probably be from twice to several times the ac- 
depending 


curacy of the standards described, 


upon the merits of the individual meter under test 


A Summary of Schelku- 


VI. Appendix. 
noff's Analysis of the Current Distribu- 
tion on a Thin Cylindrical Antenna 


Schelkunoff [1] 
the free-space current distribution on a hollow cylindrical 


obtained an approximate solution for 


antenna by analogy with a thin biconical antenna, to- 
gether with the theory of the nonuniform transmission 
line 

The antenna region and the space bevond are considered 
as two muitiple transmission lines or waveguides in tan- 
dem, each with its own distinct set of transmission modes 

In the antenna region (that space within a geometrical 
boundary sphere passing through the ends of the antenna 
the field is composed of a principal or dominant mode plus 
a complementary set of an infinite number of higher order 
modes. The principal mode or wave, is that usually as 
sociated with a two-conductor transmission line. The 
complementary waves are generate 1 in order to match the 
field (at the boundary sphere) to the field in the external 
region where no principal mode of transmission exists. 
Because of the nonexistence of physical conductors in the 
external region, the field consists solely of an infinite set 


While these sets 


of higher-order transmission modes 
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of re 


charact 


impedance as a function an expressed ap 


mately as 


The average characteristic impedance over 


~ 


A A 


dr™120( log, 


Ld 


the 


of 
distribution of 


From the theory nonuniform transmissio 


the relative the principal 


K 


current Is 


for the case in which the variation in 


Kuler’s constant 


The 


ntenna 


relative current for a | 


for 


complementary 1c 


and to a second approximation al 


is bv 


given 


Bl Bl) | om 
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e principal current amplitude on the transmission 
4/2 doublets in 


where n=O, 1, 2, 3, 


im is an integer For (2n+1 


we are interested 


) Jom+i1(Bl) and Now+:(8l) are special functions 


to J and N Bessel functions of order n+1/2 


,) 
total relative current on the antenna is given by 
m of eq 34 and 35 


nitely thin filament (A 


For an antenna composed of 
the relative current 
ition reduces to cos Sr, as is usually assumed for 
i CasC 
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Note on the Runge-Kutta Method ' 
By W. E. Milne 


A comparison is made between the 


differential equation y’=f(z,y) and a method of numerical quadrature 


standard 


tunge-Kutta method of solving the 


By examples it is 


shown that the Runge-Kutta method may be unfavorable even for simple functions f. 


One of the most celebrated methods for the 
merical solution of differential equations is the 
originated by Runge? and elaborated by 

Kutta,‘ This 
thod is usually given considerable prominence 


Nystrém,® and others. 
texts where numerical methods are discussed. 
In contrast to step-by-step procedures based 
1 formulas for numerical quadrature the Runge- 
Kutta method (as it is usually called) enjoys 
two conspicuous advantages: 
No special devices are required for starting 
computation 
2) The length of the step can be modified at 
ny time in the course of the computation without 
tional labor. 
On the other hand it is open to two major 
Hyecuions: 
The process does not contain in itself any 
simple means for estimating the error or for de- 
ting computation mistakes. It is true that 
jieberbach ® has found an expression which pro- 
des an upper bound for the error at a given step 
f the Runge-Kutta process (or more accurately, 
However this estimate 


the Kutta process 


lepends on quantities which do not appear 


tly in the computation, and therefore requires 
additional separate calculation. 
Each step requires four substitutions into 


reparation of this paper us sponsored (in part) by the Office 
tesearch 
inge, Ueber die Numerische Aufl 
Ann. 46, 167 (1895 
leun, Zeitschrift Math. Phys. 45, 23 (1900 
itta, Zeitschrift Math. Phys. 46, 435 (190 


Nystrém, Acta Societatis Scientiarum Fenr 


sung von Differentialgleichunger 


Runge-Kutta 


the differential equation. For the case of compli- 
cated equations this may demand an excessive 
amount of labor per step. 

By accident or design it happens that examples 
usually illustrate the 
Runge-Kutta method are such that the method 


chosen in textbooks to 
appears in a very favorable light. 

It is the purpose of this note to exhibit examples 
where the Runge-Kutta method does not make a 
very good showing when compared with a method 
based on numerical quadrature. For the com- 
parison we employ the commonly used form of the 
Runge-Kutta method (actually due to Kutta 
If the differential equation is 


dy 


f(x,y), 
dz ry 


the ‘step from is made by the 


formulas 


where 


h f Bes YU» 


: l 
hf (ZaF 9 


hf Pa T 


hf rant h, 


These formulas yield an approximation of the 
fourth degree in fA, that is, the expansion in powers 
of h of yas, defined by these formulas agrees 
fourth degree with the 
expansion of y, 4, obtained directly from the 


through terms of the 
differential equation 
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This is to be compared with the method? that 


makes the step from s, to r,,, by using 


th V 


then calculating 


to obtain a trial value FP of y 


Y.41 approximately from the differential equation, 
and finally computing a corrected value of y 
by 


h. 
a4 /) y + oy 


> 
» 


This method is also of fourth degree in A in the 
sense that the error terms are proportional to h 


The difference 


serves as a control, since it may be shown that the 
error of y, ignoring rounding errors and errors 
accumulated from 
mately 2¥Y 
tomary to tabulate £ and to correct the trial 
value by estimating /,,,, from the trend of the E 


When this is done and when the 


previous steps Is approxi- 


In actual calculation it is cus- 


already known 
number of places retained in y is kept such that 
ke 29 is not significant, it will not be necessary 


to recompute ¥,4; after y,4; has been corrected. 
Thus after the computation is well under way we 
need only one substitution per step 

Thus this method has two advantages 

1. By internal evidence the computation shows 
how many places in y can be accepted as reliable 
of course accumulation of rounding errors, ete., 
excepted 

2. Except at the beginning, only one substitution 
is required per step 

Two disadvantages may be cited: 

l Special devices are he essary to get started, 

2. Change of the interval A in the course of a 
computation is somewhat troublesome 


Now let us consider a particular example 


dy Sy 


»>y¥=1 when r=0. 
da q 


In each case take A= 0.1 and carry the computation 


to 2 | 


In Kutta’s method, without making sor 


dependent investigation, we have no id 


many places in y to retain. In the com 
below, we arbitrarily retained four decimal 

The results obtained by the two metho 
given in table 1 and are compared with 


obtained from the actual solution y l+a 


TABLE 1 


For this simple differential equation a comp: 
ison of the times required by the two methods 
not very significant. However, since in this « 
ample Kutta requires 40 substitutions whil 
second method requires 15 or so (depending 
exactly how the start is made), the differ 
would obviously be significant for more diffi 
equations. 

It is of some interest to make the compa 
for the equations 

dy ty 
da 
and 
dy 
ag 


with y=1 at r=0, and h=0.1 

For the first of these Kutta is in error at 
by 0.00242, and for the second by 0.0000205 
both of these cases the second method gives 
values 

The foregoing examples serve to show that 
for very innocent looking differential equat 
the Runge-Kutta method may give very 
results. Since the accuracy is in any case diff 
to ascertain, the possibility of such errors occu 
Runge-h 


is a serious indictment of the 


method. 
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A Sampling Method for Determining the Lowest 
Eigenvalue and the Principal Ejigenfunction of 
Schrédinger’s Equation 

By M. D. Donsker and M. Kac 


This is a preliminary report on a sampling method for finding the lowest eigenvalue 
and the corresponding eigenfunction of Schrédinger’s equation The theory underlying 
the procedure is explained and two cases treated numerically Although the initial results 
are encouraging, more experimentation will be needed to fully test the practicality of the 
method 


I. Introduction Carlo methods is to the problem of finding a 


Certain problems leading to complicated partial 
ntegro-differential equations have recently been 
ipproached and some actually solved by utilizing 
sampling 


ious probability techniques and 


ethods Collectively these methods have be- 
ome known as the “‘ Monte Carlo’? method 
lhe problems to which Monte Carlo techniques 
ave been applied seem to be divided into two 
types. Typical of the first type is the problem of 
itrons diffusing in material media in which the 
particles are subjected not only to certain deter- 
nistic influences but to random influences as 
In such a problem, the Monte Carlo ap- 
roach consists in permitting a ‘particle’ to play 


came of chance, the rules of the game being such 


that the actual deterministic and random features 


of the physical process are step by step exactly 


mitated by the game. By considering very large 
numbers of particles, one can answer such ques- 


tions as the distribution of the particles at the end 


of a certain period of time, the number of particles 


scape through a shield of specified thickness, 

One important characteristic of the preced- 

ng approach is that the functional equation de- 

ing the diffusion process is bypassed com- 

ly, the probability model used being derived 
the process itself 

\ more sophisticated application of Monte 


® preparation of this paper was sponsored (in part) by the Office 
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probability model or game whose solution is re- 
lated to the solution of a partial differential 
equation, or, as in the present paper, to determine 
the least eigenvalue of a differential operator by 
means of a sampling process. As an example of 
how the latter problem might be attacked, we 
quote from a paper of Metropolis and Ulam: ’ 
“For example, as suggested by Fermi, the time 


independent Schrédinger equation 
Ad(z,¥,2 r V)o Z,Y,2), 
Reintroduce time 


could be studied as follows 


dependence by considering 
u(r.y,z,t Or, Ze 


then, uw will obey the equation 


This last equation can be interpreted, however, as 
describing the behavior of a system of particles 
each of which performs a random walk, i. e., 
diffuses isotropically and at the same time is 
subject to multiplication, which is determined by 
the value of the point function V. If the solution 
of the latter equation corresponds to a spatial 


N. Metropolis and 8. Ulan The Monte Carlo method, J, Am, Stat 
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mode multiplying exponentially in time, the ex- 
amination of the spatial part will give the desired 
o(7,4,2 corresponding to the lowest ‘eigenvalue’ 
? 

The main purpose of the present paper is to 
present an alternative method for finding the low- 
est eigenvalue and corresponding eigenfunction of 
The chief difference be- 
tween the two approaches is that ours involves 
eliminating entirely the 
multiplicative This alteration in the 
model seems to simplify the numerical aspects of 
the problem, especially if punched card equipment 


Schrédinger’s equation 


only a random walk 


process. 


is to be used. Apart from the possible numerical 
simplification, the method is based on a mathemat- 


ical theory that in itself is of some interest. 


II. Mathematical Theory 


Let X,, X,, X . be independent identically 
distributed random variables each having mean 0 
and standard deviation 1 and let S,=—.X,+.N, 

X,. Under certain general assumptions 
, the most severe of which is that V(r) be 
that the limiting 


on Viz 
nonnegative, it can be shown * 
distribution function o(a,.t) of the random variable 


] — 
1 v(%) ) 
is such that 
e~**"'da(aftjdt ¥(xr)dz, 2 


is the fundamental solution of the 


where V(x 


differential equation 


l ry 


2 dx’ 
subject to the conditions 


y r >() 
y’ (xr) |< M 


LFV 


y’(+0)—yv 0 2. 


The fundamental solution ¥(z) of (3) is expres- 
sible in terms of the normalized eigenfunction 
We have since learned that Dr. G. W. King has made similar calcula- 


lations using a closely related method 
«™M. Kac. On distributions of certain Wiener functionals, Trans 


Math. Soc. 65, 1 t wi 


Am 
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' 
\(¥,(z)} and eigenvalues A, of the one-dimens ops 
Schrédinger eigenvalue problem § 


| d*y 


2 dx’ 


V(x) ¥ (2) Ay >0 { 


as 
‘ so ¥ O) YP, (a ‘ 
az : 
¥@) => Pe | 


Thus from (5 and 2 


we . wa 


Vv(x)dxr ‘g «std a(a,t)dt 


Inverting (6) with respect to s we get, 


f “lo a,f >». Mit (Q)) y¥ t dy r 


0 
. . 


and therefore we obtain the expression for \ 


| . 
; log | e~*d,a(a,t \ 


\, =lim 


If in (7) we neglect all terms in the expansior 
but the first, we obtain 


F } 0 
log é «da at ~log) ¥i(0 Vi q dx, \,f 


Thus, if, by choosing a finite t, we attempt to « ee 
culate \, from (8), we have two sources of erro HL. 
The first, usually a small source of error, is fr by 
the exponentials neglected in the expansion (7 

The second, and more important, is from neg! 


{ F ) 
ing the term (1/1 log vi Q) vi x)dx i 


latter source of error is especially significant, s 

as will be apparent shortly, it is impractical fro! 
other pomts of view to take t very large \ 
this difficulty may be obviated by considerin 


for two distinct values of t, say ¢, and f,, th 








5 From here on all the steps are formal 


they can be justified rigorously 
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Ns ona 


wel ect the exponent ials after the first as before, 
on dividing 


, e~*d,a(a,t;) 
; l 0 

X,~ j log — . 10) 
U “qd a (a,t,) 


lhe Monte Carlo process consists in the calcula- 
f e(a@,t;) and o(a,t,) by a sampling process. 
If instead of o(a,t) we consider the limiting dis- 
ribution o:(a,t) of the random variable 


ie Wee ), (11) 


also satishes (2 and 3), but now the 


ondition y’(+-0 y’ 0 2 


y'(é 2. Therefore, repeating steps 


is replaced by 


é “la (a,t 


» that we are also able to calculate the principal 
gentunction. 
lhe extension of the preceding method to multi- 
mensional Schrédinger equations is immediate. 
s in these cases that the method will probably 


a,t rove to be most useful since, unlike the standard 
variational procedures, the extension to several 
mensions seems to cause comparatively little 
liffieulty. For illustrative purposes we will con- 
Schrédinger’s equation in three dimensions. 
consider three independent 


ar Here we must 


sequences 


Z,, Z2, Zs; - 


ndependent identically distributed random 
bles each having mean O and standard devia- 
Let S..=X,+X;2 - + + +X,, and S,,, 


‘. have the obvious meanings. We consider 

the limiting distributions o(a,t) and o(a, £,'n, ¢, 8) of 
the random variables 

arch 
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l . Ss Sy 8 
Wt k<nt yz yr yi 


and 


l sf. oo , ee: 
DV (§+ 52) oto ¢+>"), 
Th k<nt yn yn yn 


respectively. 
In exactly the same way as in the one dimen- 
sional case, we arrive at 


e~*d,.a(a,t,;) 
A,~ log “. ’ 
i,—t, 
e~*d,.0(a,te) 
and 
; ' ¢ “da (a,t,n,¢,t) 
Wi (E, n, 
¥,(0,0,0 lim 
= , é a] .a(a,t) 
So far, the theory was carried out under the 


assumption that the potential function V was 
nonnegative. In most cases of physical interest 


this is not so. For the hydrogen atom, for 
instance 
const. 


very 


However, the modification is easy, although we 
have not as yet cleared up all the points of math- 
ematical rigor. 

The formula for the 


becomes 


lowest eigenvalue now 


é “lo a,t;) 


€ “doo ats 


and a corresponding modification needs to be 
made in the formula for the principal eigenfune- 
tion. We have not yet tested numerically any 
case with a negative potential function, but we 


hope to be able to report on this in the near future. 


Ill. Numerical Examples and Discussion 


The Monte Carlo procedure used here consists 
in the calculation of the distribution function 
o(a,t) by a sampling process; the principal eigen- 
value is then calculated from (10). For the pur- 
poses of numerical illustration we consider two 
Vi2)=!z|. In both of 


examples, Vi2z)=2? and 
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these cases the eigenvalues are known, and hence 
we have i check on the accuracy of our procedure. 
In order to calculate o(a,t), say when V(r) =2", 
we see from (1) that we must consider the limiting 


distribution as no of 


‘xs 14) 


2 

This means from our point of view that we 
must (a) choose a distribution for the ’s, (b 
choose a sufficiently large n, (c) select an appro- 
priate t, (4) calculate for nt \’s the normalized 
sum (14), (5) 
empirical distribution may be obtained from these 


repeat (4) many times, so that the 
many samples 
Although, 
previously, the distribution 
independent of the distribution of the Y's, the 


under the conditions mentioned 


function ofa.) is 
actual numerical calculation of o(a,.t) is expedited 
by choosing the distribution of the X's to be the 


Bernoulli distribution, i. e 
, , l 
P(X=1 P(X | = 
The sequence of random variables .Y,, X:, 
is then a sequence of l’s and L’s, such 


as might be obtained in coin tossing. This is 
conveniently and rapidly achieved on a calculat- 
ing machine by considering sequences of random 
digits, counting even digits +1 and odd digits —1. 

The value of n to be used must be large enough 
so that the empirical distribution function cal- 
culated is close to the theoretical limiting dis- 
tribution function ofa, Krom (10) we see that 


the two values of ft, f;, and ¢,, to be used in the 
calculation of A must be large enough so that the 
exponential terms neglected are sufficiently small. 
However, since the sample size is nt, the desire 
to make both n and ¢ large must be tempered by 
practical considerations. ‘The number of samples 
to be used must be large enough so that the 
empirical distribution adequately represents o(a,t 

Before discussing these points in more detail we 
consider an actual numerical computation. The 
following data for Viz r? and Viz xr were 
calculated from a certain set of random digits 
on the IBM Electronic Calculating Punch, Type 


604. For both 2? and |\z! n was selected to be 400, 


' 


the RAND Corporat 


* This ‘set of random digits was prepared by 
Santa Monica, Calil 
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%£=5, t.=3.75, 
table 1 is tabulated 


and 100 samples were us 


1H AnH 
l yo S, y l ~ th. Ss 
8000 <> * 8000 —& 
] 1M) l 20H 
B, oo. a 
; 160000 — 160000 


Since each RAND random number card ¢ 

50 random digits, and 2,000 digits are needed 
form one sample, a set of 4,000 RAND cards y 
sufficient for this experiment. It takes 20 m 
utes on the IBM 604 to calculate A, and A 

one sample (similarly for B, and B,), so that 
takes approximately 35 hours to secure the follo 


ing data 


PARLE l Expe rimental data for lowe st eigenvalue 


1 i B I 
l 1. $32 2 738 1. 20 
SOS 7 2u 8. 502 
2. 795 ‘ mn 9. 610 
i 8. 72 17 722 2. 61 
1169 4.7 IZ 6. H02 
1.19 2. 136 us 1. 453 
6.074 12.112 15, 258 9». O58 
‘ 4. 103 242 Sof 7. (04 
’ 751 S. S40 Ww). 440 1s. 14 
1a + YSl 6. 069 6.679 
249 ‘ s44 10. 044 
- 284 i m2 +097 
I 1. 8&8 2 104 1. 50 1. 7x 
4 2. O22 510 1. 638 
1. Ot ns 1. 184 5. 289 
7. 7" 2.712 1 7hy 44. Ue 
l 2s 1. oF wit 4 
Is 1 S44 2.654 1. 523 2. 18. 
r) ‘ 7 O17 Ss aM 
Po ‘ 4 6.380 7. Ole 
4.20 44 7.7. If 
22 } S47 ”) h. 666 ». 24 
= ». 2 6.314 ¥. 148 10. AO) 
ai 1. 962 2. 60S 1850 v2 
2 2 1 S07 4 is2 
»” 3.¢ 7 ol ~ 
2 ). 168 ! v2 ’ 44. #1 
as 62 sUN 11. 528 16.01 
2u 21 4.152 "2 11 
”) Pal 72 1. 136 ‘4 
! 2.71 us 3. OO 4. 02 
2 +. £38 ine 6.114 17. 41 
7. G28 11. M7 5. 184 O01 
‘ 47 su 4. Au ‘ “ 
1. 958 2. 504 1. 457 1 O48 
i 1 ale & S88 11.7% 
7 207 6. 38 + 800 11 
is 1 S83 2. 682 1. 387 1. 905 
a] 4,549 10. 020 &. 308 30, 233 
“ 2. 282 067 2. 258 2. S82 
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Experimental data for 


i i 
4! 4.2 6. 957 f 
42 7. 559 S AO4 17 
45 122 87 
4 1. 435 ti 
i 2. 984 477 2 
4 1. 677 158 l 
‘7 14 + 697 il 
4s 1. GR2 v.14 ) 
49 1. 675 152 l 
"a 1. 580 l 
l 2. te 21 2 
2 4. i SS. 
53 7 } 12. O92 19 
54 bs Tt 
5 4 ‘ 4 
‘ 1. 42 { 
7 Wut l $58 1 
s 1.916 44 l 
; 1. 490 2 
r) ' 82 4 
1 s t 
#2 74 ; 
9g 
i ‘ 
‘ tine x 
wt ‘ x. 7S ‘ 
“7 +. 4) +s t 
Hin | +. 164 l 
u ‘ 
os ® 
2. 4s be | 
O79 a @ 9 
74 s4 “M4 ; 
rf 1s s 
t ] s ] 
41 Su 
s m ‘ 
; 2 ¢ “ 
~ uv 
s! 1 o s 
82 ! s 
s ua 
s4 2 ON4 NTr 
s 74 ] 
st SIS i. 444 l 
aT 4 387 S44 i) 
SS tw 244 ‘ 
sy lil l ; 
w t ~~ » is ! 
’ is nM ) 
2 i 
‘ (is ‘ ‘ ‘ 
4 ’ 4 O04 
4 614 “ 
) 1143 l Is 
17 4.108 s2 
is ‘ 44 7 
au auUA iyi 1 
wn) 6S 4 
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143 
aM) 


SHU 


> BAI 


be | 


by 





lowest eige nvalue 


By 


11, 4s 


14.80 





Con 





Inordertocalculate A, from(10) we need the values 


of e~*d,a(at and e~*d,a(a,t.). Both of these 


integrals were calculated numerically from the 
data in table 1 by adding the exponentials of the 
entries in appropriate columns. For the case 
V(r) = 2, the values of the integrals were obtained 


100 100 
from 1/100 53 e~4" and 1/100 53 e>** and for V(z) 
l l 
100 10) 
z* from 1/100 53 e~*" and 1/100 53 e~*". Inthe 
l l 


case V(r)= 27! the true lowest eigenvalue to two 
places is 0.81, and in the case V(r) =z? it isy2/2 
0.71. 
First Second ALL 10 
imp) I imi 
r 0. 83 0. 79 0.81 
I SU ou io 


An interesting feature of the data is that column 
B,, for example, may be thought of as arising from 
n= 2000, t=1; n=1000, t=2, etc., as wellas from 
n=400, t=5. The larger we take n, the closer 


will be the empirical distribution to ¢(a,f The 


value of e “d,a(a,t) in the ease f=1, and Via 


} 
. 


r? can be calculated exactly and, to three places, 
is 0.678. Calculating this integral from column 
B, with n= 2000, t=1, we get 0.685 

Instead of using Bernoulli distributed variables, 
one might use other distributions. One definite 
advantage of Bernoulli distributed variables is 
that the computation utilizes only the crudest 
properties of the random digits, i. e., whether they 
are even or odd. One possible advantage for cer- 
tain other distributions is that » might not have 
to be taken so large. In particular this should 
be true if we use Gaussian distributed variables 
RAND Gaussian deviates were used in constructing 
table 2. Here we chose f;=3.75, h=5 as before, 
but now n= 100. This means samples of size 500 
instead of 2000 and therefore a total machine 
computation time of 9 hours for 100 samples. 

Envisaging the possibility of calculating the 


second eigenvalue, we considered the quantities 


10 >> (S,+5)? and 10°* 53 (S,+-5)’, 15) 


rl 
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Expe 


mental data 


lo 


second 


eige nval if 


Experimental data for second eigenval 


which correspond to taking &=0.5 in (11 

should not make any difference in the caleu 
of the lowest eigenvalue, and consequently t 
can be utilized in the same way as tabl 

should however be borne in mind that co 
C’, and C, of table 2 represent experimental \ 
of the quantities (18) with the S;,’s being st 
Gaussian deviates. Asterisks on the enti 
table 2 indicate that Sy; or Sip was neg 
Although this information is unnecessary f 
purpose of calculations of the lowest eigen 


it is used in the calculation of the next eigen 


How this can be done is explained brief 
section IV 
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the data in table 2 and again (10) we 


ted the 
value—0.71 


following approximations to ), 


both Via I and Vir r*, all the 
values are known so that in the two illus- 
e examples above it was easy to choose 


priate values for ¢, and t,. The proper 

for n and the appropriate number of samples 
s determined experimentally. In general when 
“rule 


Hav- 


first made a guess atl an appropriate 7 and t, 


genvalues are unknown the following 
imb”’ procedure seems to be feasible. 
a certain number of samples, compute Ay. 
it the computation now for the same n, 
same number of samples and successively 
er t's until the calculated values of 4, become 
If they 


the number of samples 


to the desired number of places 
not become stable 
be increased Ke ping two values of t, for 
h the calculated \’s had the stable value, 
value of A 
large. If it 
stability 


the calculated 
is sufficiently 
until a new 


n and see if 

If not, 7 
change, Increast } 
ars. This stable value is then the appro- 
te approximation to \,. The mere fact that 
lity is observed means the number of samples 


fliciently large 
The Second Eigenvalue 


the principal eigenfunction is even, then it is 
ble to extend the theory of section II in 
that the calculation of the second 
Without going into 


result 


al way 
nvalue becomes feasible. 
just state the 


details we pertinent 


Determination of Eigenvalues by Sampling 


lim Eye 


and 


\ > : 
lim Eye ; " i # 0? 


the mathematical expectations on the right hand 
sides being conditional expectations under the 
0 and S,, 


conditions S,, 0, respectively. 


We then have 


WwW here 
2] 


0) 


log . 19 
/ EK. 7 

Krom the discussion of section I] it should be clear 
how one applies (19) when the data of table 2 are 
available. One must see to it that £ is so chosen 
that ¥.(£) #0, otherwise a higher eigenvalue may 
have been calculated From the data of table 2 


we obtain 


whereas the exact value is 3/2 v2 


2.12. The poor 
agreement could have been expected in view of low 
accuracy in the calculation of \,. 

thank E. C. Yowell, of the 


Standards, for wiring the 


In conclusion, we 
National Bureau of 
boards and for excellent supervision of all the 


punched card work 


Los (AN 


ELES, September 12, 1949 





